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SELECTIVE REFLECTION FROM MERCURY AND 
CADMIUM VAPORS! 


By H. L. Wetsu, J. KASTNER’, AND A. C. LAuRISTON® 


Abstract 


Selective reflection from mercury and cadmium vapors in the neighborhood of 
the resonance lines was investigated quantitatively using incident light of con- 
tinuous spectral distribution. For a given experimental arrangement the lowest 
vapor density (atoms per cubic centimeter) at which selective reflection could 
be detected at Hg 2537A, Cd 2288A, and Cd 3261A was inversely proportional 
to the oscillator strength (f-value) of the absorption line. Contours of the selec- 
tive reflection of the Hg 2537A and Cd 2288A lines were obtained over a large 
range of vapor densities up to 80 X 10'8 atoms per cc. At this density the cad- 
mium reflection extended over several thousand cm.—! compared to about a 
hundred cm.—! for mercury. The general features of the reflection contours can 
be explained by the theory of reflection from an absorbing medium. Fitting a 
theoretical curve to the experimental reflection contour yields values of the 
oscillator strength and the damping constant. For Hg 2537A and Cd 2288: 
the f-values are 0.0268 and 1.40 respectively. in good agreement with those found 
by other methods. The damping constani, ‘y, varies as the square root of the 
number of atoms per cubic centimeter, indicating that the mechanism of selective 
reflection is not the same as that for absorption and emission. An empirical areal 
law for selective reflection confirms the relationship y « ~/N. A shift of the 
resonance frequency of the order of magnitude of the coupling shift calculated 
by Weisskopf was observed for Hg 2537A. Deviations from theory at both high 
and low pressures were observed for Cd 2288A. The low pressure deviation takes 
the form of a line of residual intensity at the resonance frequency, which may be 
due to a different kind of selective reflection. 


Introduction 


Selective reflection of light from a metallic vapor was first observed by 
Wood (20). When a quartz bulb containing mercury vapor at room temper- 
atures was irradiated with the unreversed resonance line 2537A, the vapor 
throughout the whole bulb emitted the resonance line in fluorescence. As the 
vapor pressure was increased the fluorescence was confined more and more to 
the surface of the vapor in the path of the incident light. The surface fluor- 
escence diminished in intensity with further increase in pressure until it was 
no longer observable. However, in this region of vapor pressure (c. 80 mm.), 
Wood observed that the reflection of the resonance line from the quartz- 
mercury vapor surface was enhanced and that its intensity increased steadily 
with increasing vapor pressure. The polarization characteristics of the selec- 
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tive reflection were qualitatively the same as those of light reflected from a 
metallic surface. It was also found that the immediate neighborhood of the 
line was a region of anomalous reflection, with high intensity on the short 
wave-length side and low intensity on the long wave-length side. Wood pointed 
out that this intensity distribution could be explained by the variation of the 
refractive index of the vapor in the neighborhood of a strong absorption line, 
as given by anomalous dispersion theory. 


Resonance fluorescence and selective reflection are of quite different natures. 
Rump (13), using the sharp line from a primary mercury resonance lamp, 
showed that the half-width of the line in selective reflection was. the same as 
that of the incident line, whereas the half-width of the line in fluorescence was 
dependent on the temperature of the vapor. The effect of foreign gases on the 
intensity of the reflection was investigated by Schnettler (14). The reflected 
intensity decreased only 15% when two-thirds of an atmosphere of hydrogen 
or carbon dioxide was introduced into the tube containing the mercury vapor. 
It is well known, however, that only small pressures (c. 1 mm.) of these gases 
are necessary to produce considerable quenching of the resonance fluorescence 
at room temperature. Schnettler therefore concluded that a finite lifetime of 
the excited state cannot be involved in the process of selective reflection. This 
view is confirmed by the recent work of Hansen and Webb (7). When a mercury 
vapor surface was illuminated intermittently with resonance radiation, the 
selective reflection, unlike the resonance fluorescence, did not persist between 
the light pulses. The resolution of the timing method was about one-fortieth 
of the lifetime of the excited state of the free atom. 


It is apparent that selective reflection has been investigated much less exten- 
sively than the related phenomenon of resonance fluorescence, although its 
study might be expected to yield useful information concerning the interaction 
of light and matter. Regular reflection is by its nature a coherent process. 
Classically, partial coherence is introduced into resonance fluorescence when 
the distance between the atoms becomes small compared to the wave-length 
of the light. Theoretical discussions on this basis have been given by Burton 
(2) and Fermi (4). An explanation in terms of the quantum theory is more 
difficult. Weisskopf (16), in the extension of his quantum theory of resonance 
fluorescence to gases at high densities, has deduced that the fluorescence has 
a definite phase relationship with the incident light. He states that this explains 
the selective reflection observed by Wood. The theory, however, is not com- 
plete enough to explain in detail the experimental results. 


In the present experimental investigation attempts were made to observe 
selective reflection from other metallic vapors. For cadmium the reflection was 
found to be easily observable at the resonance line 2288A. The other resonance 
line at 3261 A was also observed in reflection at high vapor pressures (c. 20 atm.). 
The effect could not be observed in thallium because of the difficulty of con- 
taining the vapor in a quartz tube at very high temperatures. 
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Some of the possibilities of a quantitative study of selective reflection were 
explored. In particular the Hg 2537A and the Cd 2288A reflections were 
studied quantitatively over a wide range of vapor pressures. The variation of 
the reflected intensity with frequency could be interpreted in a general way in 
terms of the classical anomalous dispersion equations. There remain however 
certain features of the experimental results which do not appear to be expli- 
cable in terms of existing theories. 


Experimental Procedure 


Since selective reflection is not confined to the resonance frequency alone, 
the selective reflections were detected and studied by using incident light of 
continuous spectral distribution and analyzing the reflected light spectro- 
graphically. The region of anomalous reflection at Hg 2537A is very narrow 
for vapor pressures up to a few atmospheres, so that it was necessary to use 
a spectrograph of high dispersion. At Cd 2288A the region of anomalous 
reflection is many times broader and a low dispersion spectrograph could be 
used. 


A reflection cell with a window of good optical quality was necessary to 
obtain sufficient intensity for the study of the mercury reflection at high dis- 
persion. The cell used was originally a Hilger absorption cell of fused quartz, 
1.4 cm. in internal diameter and 4 cm. long, with contacted plane windows. 
The outer surface of one window was ground and polished at an angle of 1° 
with the inner surface to give an angular separation of the reflections from the 
two surfaces. Mercury was distilled under vacuum into a small side tube, and 
the cell was then sealed off. A nichrome oven with a plane window of fused 
quartz was used to heat the cell, which was supported by a solid quartz rod 
fused to the cell and extending through the wall of the oven to a clamp outside. 
The oven was differentially wound so that the front window of the reflection 
cell could be maintained 20° to 30°C. hotter than the side tube to prevent 
condensation of mercury vapor on the window. Thermocouples were attached 
to the front of the cell and to the side tube containing the liquid mercury. 


A low voltage hydrogen arc of the type described by Allen and Franklin (1) 
furnished the ultraviolet continuum. A condensing lens of quartz with the 
small luminous portion of the hydrogen lamp at its focal point gave a parallel 
beam falling on the quartz—mercury vapor surface at almost normal incidence 
(Fig. la). The reflected beam was focused on the slit of the spectrograph by 
the same lens. A double-prism quartz spectrograph with Littrow mounting 
gave a linear dispersion of 1.7A per mm. in the region of 2537A. With a spectral 
slit width of 1 cm.— the reflection spectrum could be photographed in one hour 
on Ilford Special Lantern plates. 


Several reflection cells were made up for use with cadmium. These were 
constructed by fusing a prismatic window with an angle of 5° to a fused quartz 
tube, 1.2 cm. in internal diameter. Some distortion of the windows was intro- 
duced in the fusion process with a consequent loss in intensity in the reflected 
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image of the lamp. Cadmium was distilled under vacuum into the cell which 
was then sealed off. A differentially wound platinum oven, designed for use 
up to 1400°C., was used to heat the cell. A condensing lens placed between the 
hydrogen arc and the cell gave a converging beam reflected from the cell and 
focused on the slit of the spectrograph (Fig. 1b). A single-prism quartz spectro- 
graph with an f/4.5 camera lens and a linear dispersion of 20A per mm. in the 
region of 2288A was used. With Ilford Q.1. plates exposure times of 30 min. 
were required for a spectral slit width of 11 cm.~}. 
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Fic. 1. Optical arrangements used for selective reflection (a) from mercury vapor, (b) from 
cadmium vapor. A—light source; B—lens; C—reflection cell; D—prismatic window; E—slit of 
spectrograph; F—image reflected from front surface of the window; G—metal. 


The photographic plates were calibrated for intensity measurements by 
photographing the continuous spectrum of the hydrogen are with a stepped 
slit substituted for the usual slit of the spectrograph. Tie calibration spectra 
were photographed with the same exposure times as the .eflection spectra. A 
Leeds and Northrup microphotometer was used to measure the densities of the 
spectrograms. 

The spectrum of a cold mercury are was photographed in juxtaposition with 
each spectrum of the selective reflection of mercury. When the two spectra 
were microphotometered with fiducial marks superposed the position of the 
sharp mercury line at 2534.8A could be transposed to the reflection spectrum. 
It was found from a series of plates taken at the same mercury vapor pressure 
that the calibration reproduced to 0.025A. The cadmium reflection spectra 
were calibrated in a similar manner with a nickel arc comparison spectrum. 


Vapor pressures for mercury were obtained from the formula given by 
Menzies (10), while those for cadmium were taken from the review of Ditch- 
burn and Gilmour (3). Vapor pressures for cadmium have been measured only 
up to 840°C., so that the higher pressures which were obtained by extrapolation 
cannot be expected to be accurate. The number of atoms per cubic centimeter, 
N, was found from the formula, N = p/kT, where T is the temperature in 
degrees Kelvin at the window of the cell. The pressure, », of the vapor in the 
cell was taken to be the same as the vapor pressure of the liquid in the side tube. 
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This procedure is justified according to Weber (15), when the mean free path 
of the atoms is much smaller than the diameter of the connecting tube. 


Selective Reflection from Cadmium Vapor 


Although earlier investigations on selective reflection have been confined to 
the mercury resonance line 2537A (61S o-6*P;), the effect should be observable 
for any strong absorption line, and the intensity of reflection should show a 
relationship to the absorption strength of the line. The second mercury reson- 
ance line at 1849A (6!S)-6'P;) has a much greater absorption strength than 
the intercombination line 2537A, but it lies in a wave-length region incon- 
venient for observation. However, the corresponding resonance lines of cad- 
mium lie at longer wave lengths and can be readily investigated. 


The hydrogen continuum reflected from the cadmium cell was examined for 
cell temperatures up to 1125°C. Selective reflection of the resonance line 2288A 
(5'So-5'P,) was first observed at a vapor pressure of 4 mm. as a faint sharp line 
just discernible on the continuous background of the spectrum. With increas- 
ing vapor pressure the region of anomalous reflection broadened and assumed 
the characteristic appearance of selective reflection, with a less than normal 
intensity on the long wave-length side and a greater than normal intensity on 
the short wave-length side (Fig. 2). At the higher pressures (5 to 20 atm.) the 
region of the anomalous reflection was several hundred Angstroms wide com- 
pared to a maximum width of about 10A for the Hg 2537A reflection at 8 atm. 
The cadmium resonance line at 3261A (5!S\-5°P;) was first detected in the 
reflected spectrum at 1125°C., the highest temperature reached. This temper- 
ature corresponds to a vapor pressure of 15,000 mm. as obtained by extra- 
polation of the vapor pressure equation. With a similar experimental arrange- 
ment the selective reflection of the Hg 2537A line was detected first at a vapor 


pressure of 170 mm. 
TABLE I 
RELATION BETWEEN f-VALUE AND N FOR FIRST APPEARANCE 
OF SELECTIVE REFLECTION 


Resonance Oscillator | No. <= Nf, 
line | strength, f OS 10-17) | (xX 107-17) 
Camek fu | 08 | 06 
Cd 32614 | 1.9 x 10-3 | 1000 ee 
Hg 2537A | 2.68 X 10> | 28.5 © 0.76 

1 


The conditions for first appearance in reflection of the two cadmium lines 
and the mercury line are given in Table I. The f-values or oscillator strengths 
of the absorption lines, obtained from the summary of Mitchell and Zemansky 
(11), are given in the second column. The value of the product, Nf, is found 
to be approximately the same for the two lines for which the data are most 
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reliable, and is not greatly different for the other line. This gives an empirical 
rule which might be used to predict the vapor density of the first appearance 
of selective reflection for any absorption line whose f-value is known. 


Contour of the Mercury Reflection at 2537A 
EXPERIMENTAL RESULTS 


Reduction of photographic densities to intensities over the whole contour of 

. the reflection was carried out for a series of plates at mercury vapor pressures 
of 3, 1,2,4, and 8atm. The intensity contours of the selective reflection plotted 

against the distance in cm.~! from vo, the frequency of the resonance line, are 

shown in Fig. 3. Relative intensities are given as ‘‘% reflected intensity’’ since 
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Selective reflection in the neighborhood of Cd 2288 A for vapor pressures of (a) 34 mm., 


Zatm., (¢) 14.7 atm. Nickel arc comparison spectrum. Wave-length region 2200 to 24004. 
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this quantity is used in the interpretation. The ‘‘% reflected intensities’’ are 
based on a calculated value for the ‘“‘normal reflection”’ outside of the region of 
anomalous reflection, where the refractive index of mercury vapor at the pres- 
sures used here can be taken as unity according to the dispersion formula given 
by Wolfsohn (19). The normal intensity of reflection was calculated from the 
Fresnel reflection formula for normal incidence to be 4.06% at 300°C. and 
4.08% at 500°C., using refractive index measurements for fused quartz by 
Ramachandran (12). 


INTERPRETATION OF THE RESULTS 


The main features of the contour of the selective reflection can be explained 
qualitatively, as Wood pointed out, by the anomalous behavior of the refrac- 
tive index of the vapor at frequencies near the resonance frequency. The most 
obvious approach to the interpretation of the experimental results is therefore 
through the classical theory of reflection from an absorbing medium. The 
refractive index, n, and absorption coefficient, k, are given as functions of the 
frequency, v, by the dispersion equations: 


n—k =n + SS a (1) 
2rmvo (vo << v)? + (y 4r)? 
ida, ye ta (2) 


Sx?mvo (vo — v)? + (y/4r)? 


in which vo is the resonance frequency, mo is the hypothetical refractive index 
if no anomalous dispersion were present, y is the damping constant, f is the 
classical oscillator strength (related to the quantum mechanical transition 
probability), and the other symbols have their usual meanings. 


The reflecting power of an absorbing medium for light at normal incidence 
is given in terms of m and k by the Fresnel formula: 


ae Mere (3) 
(n +n’)? + Rk? 


where n’ is the refractive index of the bounding medium (fused quartz in the 


present case). 


Calculation of the f-value from the Wings of the Contour 

The absorption coefficient, k, decreases rapidly with increasing frequency 
displacement, Av = v — vo, from the resonance frequency. Hence, in the wings 
of the contour the refractive index of the vapor can be calculated from the 
reflected intensity by means of Equation (3) with k = 0. Also, when k is 
negligible the dispersion equations (1) and (2) reduce to the form: 
oe es (4) 


——= 9 


4? — 
2armvy Av 
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in which mp has been put equal to unity. A plot of 1/(m?—1) against Av should 
therefore give a straight line in the region in which k is negligible, and the slope 
of the line can be used to calculate the f-value of the transition. Fig. 4 is an 
example of such a plot for mercury vapor at 8 atm. pressure.’ The fact that the 
best straight line does not pass through the origin indicates a change of vo with 
pressure. The f-values obtained for various vapor pressures by this method 
ranged from 0.023 to 0.027, in fair agreement with values obtained by other 
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methods (11). However, since the calculated values of 1 are very sensitive to 
small errors in the measured reflectivity in the wings of the pattern, this method 
is not the most accurate way of determining the f-value from the reflection 
contour. 


Fitting of a Theoretical Curve to the Complete Reflection Contour 
The explanation of selective reflection on the basis of anomalous dispersion 
can be tested quantitatively by attempting to fit the experimental contour by 
a variation of the parameters f and y in Equations (1), (2), and (3). The pro- 
cedure used in fitting a theoretical curve to the experimental contour was as 
follows: The dispersion equations are written in the form 
Mx 


n?— B= nz —— 


1+ x? 





=A (5) 
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where 


Qnk = ae 
1l+x 


C= oe, and M 
+ 


The solutions for 2 and k are 


2-_VAT+B+A pp 


2Nfet 
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_VA+ B- A 






































n? 
2 2 
| } 
9.0 -— oo 
| Be 
| | | | 
| | 
etn fon - 
Selective Reflection | 
from Mercury Vapour 2 | 
at % atmosphere pressure — #07 [—— ee ane 
¢ Experimental Intensity Se 
Theoretical Intensity {M*0.% | 
sak Theoretical piensa Ub “aa . 3 | 
modified for a slitwidth “|: 3 
of 1.0 cm:! ie & 
' ee 
| e i oe | 
Soar nonce seers 4.06 fo =4.03 = 
zs i] 
| AN? | | 
+___— +— 30 %+4— } | 
| y 
| | | 
—_}—___——_+—— 29-4} and 
| 
a 
——— t— 10+ 
| | Oo 0e=0.em-* 
0, Calibrated 
15.0 10.0 5.0 -O+ 5.0 10.0 15.0 
Ad cm." 





Fic. 5 





101 


(6) 


(7) 
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which can be used to obtain R from Equation (3). It was found convenient to 
obtain the fit by varying the quantities M and 42/y rather than f and y. Since 
the spread of the theoretical pattern along the Av axis is determined by 42/y and 
along the R axis by M, the two quantities can be varied independently to obtain 
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the fit. The experimental contour was plotted against log Av and the theoretical 
contour for a given M against log x on separate sheets of graph paper. A test 
for fit should have been possible by shifting the theoretical and experimental 
curves relative to one another along the abscissa axis, since this corresponds to 
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a variation of 44/7. However, the change in »o with pressure adds a compli- 
cation, the shape of the experimental curve being extremely sensitive to the 
choice of zero of the Av scale. A family of experimental curves was therefore 
plotted corresponding to shifts in the resonance frequency from 0 to 1.5 cm.~! 
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in steps of 0.1 cm.—! from vo, the calibrated frequency of the resonance line. 
For the correct choice of M it was then possible to fit a theoretical contour 
(R against log x) to one of the family of experimental curves (R against log Av). 


Examples of the results obtained by this procedure are shown in Figs. 5, 6, 
and 7. In each case the agreement was improved at the maximum and minimum 
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when the theoretical curve was modified for the finite slit width of the spectro- 
graph. The agreement was excellent for the higher vapor pressures but became 
progressively worse in going to lower pressures. This may be due in part to the 
larger experimental error at lower pressures where the reflection pattern is very 
narrow. However, the deviation is always such that the experimental inten- 
sities are higher than the theoretical intensities in the vicinity of vo. This could 
mean that there is an additional intensity around vo which is not accounted for 
in the interpretation of the reflection contour based on the anomalous disper- 
sion equations. This residual intensity is much greater in the case of cadmium 
at low pressures and will therefore be discussed at greater length later. 


Derived Values of f, y and éyo 


The fitting of a theoretical reflection curve to the experimental contour for 
a given vapor pressure yields values of f, y and dvo= vo— vo’, where vo and vo’ 
are respectively the calibrated and calculated values of the resonance fre- 
quency. The results for five vapor pressures are given in Table II. 


TABLE II 
VALUES OF f, Y, 6 ¥o DERIVED FROM THE REFLECTION CONTOUR OF HG 2537A 


| 
| 











Vapor | No. of aie | Oscillator | ae | Y/NN, a ¢ 
preaeeits a 10-18) strength, f | aie tsi | (X 10-') | om 
O58 | 6.03 | 0.0269 ee “1.09 Md eae, Pee 
1.0 | 11.1 0.0269 | 1.42 | 4.23 | —0.2 
2.0 | 22.2 | 0.0267 | 1.96 bs gles Tae 
40 | 444 | oom | 200 | so | -o8 
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The f-value remains constant with increasing vapor density, and the mean 
value, f = 0.0267 + 0.0003, is in good agreement with values obtained by 
other methods. These have been summarized by Mitchell and Zemansky (11), 
who conclude that the most reliable of the older values is 0.0268 found by 
Garrett (6) from the decay of electrically excited resonance radiation. 


Selective reflection thus provides a new means of measuring the f-values of 
strong absorption lines. The method seems to be comparable in accuracy with 
the best of the older methods and could undoubtedly be improved by more 
accurate intensity measurements. 


The variation of the damping constant, y, with vapor pressure is of special 
interest. If the y measured in this experiment were the same as the damping 
constant due to collision broadening and given by the half width of the absorp- 
tion line, one would expect a direct proportionality between y and N. This 
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relation has been derived theoretically by Weisskopf (18), Margenau (9), 
Furssow and Wlassow (5) and others. However, the y in selective reflection 
varies more closely as~/N as shown by the fifth column in Table II. This result 
is confirmed by an empirical areal law discussed below and also by the results 
on cadmium. Thus the damping present in selective reflection seems to differ 
fundamentally from collision damping in the broadening of spectral lines. 


The shifts in the resonance frequency, 6vo, show a direct proportionality with 
N. However, since the frequency calibrations of the spectrograms were re- 
producible only to about 0.2 cm.~', the values of 6vp can have a large percentage 
error. Weisskopf (17) has shown that the coupling of atoms in the Lorentz- 
Lorenz sense should cause vp to move to lower frequencies at higher densities 
and the displacement of vp should be proportional to N. The frequency shift 
is given by efN 


6vo (in cm.~') = — — 
6rec’mvo 


(9) 


The reality of the effect cannot be checked in absorption experiments because 
the small shift is masked completely. by the resonance broadening of the ab- 
sorption line. Although the shifts calculated from Formula (9) are only a little 
greater than half of those observed, it is considered that 6vo can be interpreted 
as a shift due to coupling. The corresponding change in y due to coupling, as 
derived by Weisskopf, is much too small to be observable here. 


An Empirical Areal Law for Selective Reflection 

It is well known that the integral of the absorption coefficient of a spectral 
line is proportional to N and is independent of broadening processes. A similar 
relationship was found to hold for the contour of the selective reflection. When 
the area between the contour of the reflection and the straight line corres- 
ponding to the normal reflected intensity was plotted against N a straight line 
was obtained (Fig. 8). The existence of such an empirical law implies a certain 
relationship between y and N, since the f-value appears to be independent of N. 
An attempt was therefore made to integrate the metallic reflection equation to 
determine the variation of y with N, but this was not found to be possible. 
However, a numerical evaluation of the integral 


© |, — 1 20R)2 2 
| | (m_— 1.506)? + RP _ 0.0407 | dv 
-o|(n + 1.506)? + R? 


was carried out* for two widely separated vapor pressures. The variation of n 
and k with » was calculated by assuming f = 0.0267 and y = 40 ~/N, these 
values representing an average of the experimental results. The computation 
confirmed the empirical areal law. The observed proportionality of the area 
under the reflection contour to N is therefore consistent with a variation of 
vy as- VN, as indicated in the results obtained in fitting a theoretical curve to 


*The authors are indebted to the Computation Centre, University of Toronto, for assistance in this 
calculation. 
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the reflection contour (Table II). Indeed the areal law can be looked upon as 


the more convincing proof of the relation y « N since it is more direct. 
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Contour of the Cadmium Reflection at 2288A 


Examples of the reduced reflection contours obtained for the Cd 2288A line 
are shown in Fig. 9. The large change in background density over the region 
of the selective reflection made the reduction to ‘‘% reflected intensity’’ more 
difficult at the higher vapor densities. From the shape of the background 
density at low pressures empirical rules were established with which the reduc- 
tion was effected. The ‘‘normal reflected intensity’’ of 4.27%, calculated from 
the refractive index of fused quartz at 2288A assuming a refractive index of 
unity for the vapor, was used for the whole range of the reflection. The error 
thus introduced was probably small compared to the error due to the uncertain 
background density at high vapor pressures. 


Theoretical contours were fitted to several of the experimental curves by the 
method previously outlined. Since the contours for the lower vapor pressures 
were distorted by the low purity of spectrum obtained with the small spectro- 
graph, only the results for pressures higher than 350 mm. were treated quan- 


~ 


titatively. Examples of the results obtained are given in Figs. 10, 11, and 12. 
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In the medium pressure range (1-2 atm.) the experimental curves could be 
accurately reproduced, but at lower and higher pressures only the outer por- 
tions of the pattern could be fitted. 
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TABLE III 
VALUES OF f AND Y DERIVED FROM THE REFLECTION CONTOUR OF Cp 2288A 














a | ee) See constant, | YAN 
atm. (X 10-8) strength, f ¥(X 10-1) | (X 10-*) 
‘el ss Ce 11.0 a 
0.46 3.32 1.38 16.0 8.8 
0.87 | 5.96 | 1.36 21.0 8.6 
2.0 12.4 | 1.38 | 30.0 8.6 
5.4 | 31.1 | 1.36 48.0 8.7 
4.7 | 77.7 | 1.38 | 77.0 8.8: 
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The values of f and y obtained from the fittings are given in Table III. The 
f-value remained constant except at the lowest pressure where the experimental 
error was rather large. The average f-value obtained, 1.40, yields a mean life- 
time of 1.7 10~° sec. for the 5'P; state of cadmium. The value 2.0 107° sec., 
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was obtained by Kuhn (8) from magnetorotation at the edges of the resonance 
line, and also by Zemansky (21) from total absorption of the line in 
cadmium vapor at low pressures. The f-value obtained from selective reflection 
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is in fair agreement with this value and is probably rather more accurate. From 
Column 5 of the table it is seen that the damping constant, ¥, is directly pro- 
portional to~/WN as it was with mercury. The value of y is about 20 times as 
great for Cd 2288A as for Hg 2537A at the same value of N. The large damping 
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constant is chiefly responsible for the great extent of the selective reflection of 
cadmium, as may be seen from Equations (7). 


No shift in vp was required to obtain a fit between theoretical and experi- 
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mental contours as was necessary for mercury. A shift of the size predicted by 
Weisskopf’s formula, Equation (9), would have had little effect on the fitting. 


The area enclosed by the reflection contour and the straight line representing 
normal reflected intensity is plotted against N in Fig. 13. As in the case of 
mercury the area is directly proportional to N over a large range of vapor 
densities. This indicates, as before, the validity of the relationship y « VN. 
At the highest pressures where it was impossible to determine the full extent 
of the selective reflection, the measured areas are too small since they do not 
include the whole pattern. Hence, it is uncertain whether the deviation from 
the straight line at the higher densities is real. 
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For the lower pressures, when the theoretical curve giving the best fit in 
the wings was subtracted from the experimental curve, a residual intensity was 
obtained forming a “‘line’’ approximately symmetrical about vo with a half 
width of about 12 cm. (Fig. 10). At pressures greater than 1 atm. the residual 
intensity is scarcely detectable. Some indication of the effect is present in the 
mercury reflection contours but is of the order of magnitude of the experi- 


mental error. 
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The deviation from the theoretical curve for high vapor pressurés of cadmium 
is of an entirely different nature (Fig. 12). The region of the deviation is much 
greater and neither the maximum nor the minimum is as pronounced as is 
required by the theory. No variation of f and y can lead to a good fitting of the 
experimental contour. With increasing pressure above 8 atm. the maximum 
of the curve begins to fall instead of continuing to rise as the theory predicts. 


The results of the present investigation lead to the conclusion that the 
general features of selective reflection conform to the classical theory of re- 
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flection from an absorbing medium. From the viewpoint of the quantum theory 
selective reflection must be a coherent scattering process rather than an ab- 
sorption and re-emission with a finite lifetime of the excited state. The relation 
y «WN confirms the view that a simple absorption process is not involved 
although no theoretical explanation of the relationship presents itself. For the 
Cd 2288A line, which has much larger f and y values than Hg 25374, devia- 
tions from the anomalous dispersion theory are more marked. The residual 
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intensity observed at lower vapor pressures may be a second type of selective 
reflection arising from resonance fluorescence in a definite phase relationship 
with the incident light as deduced by Weisskopf (16). It should be noted that 
since there are no metastable states near the 5'P state of cadmium the quench- 
ing of the resonance fluorescence with increasing vapor pressure should be much 
less than for Hg 2537A. The reason for the marked deviation of the reflection 
contour from theory at the higher vapor densities.is not immediately apparent. 
A tentative explanation is that the anomalous dispersion equations with a 
unique value of the resonance frequency, vo, are no longer applicable when 
strong coupling exists between the atoms. 
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INTEGRATOR-EXPANDER CIRCUITS FOR THE CONTROL OF 
THE MAXIMUM X-RAY ENERGY FROM THE BETATRON' 


L. Katz, A. G. McNamara,” P. A. Forsytu, R. N. H. HAstam, 
AND H. E. JoHNs 


Abstract 


The electronic circuits known as an integrator-expander were constructed for 
the accurate control of the X-ray energy of the University of Saskatchewan 
betatron. The over-all design and operation of the circuits are discussed. The 
design of the energy bias control, amplifier, and trigger portions of the circuit is 
original, and these are described in some detail. 


Introduction 


The betatron obtained by the University of Saskatchewan in 1948 was the 
22 Mev. Allis-Chalmers machine designed primarily for industrial radio- 
graphy. No accurate means of energy control is incorporated in this make of 
betatron. For nuclear research purposes the maximum X-ray energy must be 
continuously variable over a wide range, and be quite accurately known at all 
settings. 


The so-called integrator method (1, 2, 5) of controlling the maximum X-ray 
energy from the betatron has been developed at the University of Illinois. 
The basic principles used by us are the same as those worked out at Illinois. 
The integrator built at Saskatchewan makes use of circuits which the authors 
feel lead to particularly stable operation. A signal (voltage) is obtained which 
is proportional to the instantaneous energy of the accelerated electrons so that 
expansion of the orbit may be made to coincide with a predetermined electron 
energy and therefore maximum gamma ray energy. 


This method is based on the working principle of the induction electron ac- 
celerator as expressed by one of Maxwell’s equations, 


>> aB 4 
© E.dj = - | fas ; 
Ji s ot 


From this it may be shown that the instantaneous energy is proportional to 
the time integral of the rate of change of flux through the orbit; that is, to the 
total flux linking the orbit. Since the flux linking the orbit is proportional to 
the flux linking the main coils, the voltage across the coils will be proportional 
to the rate of change of flux and, when integrated, will represent the instan- 


taneous energy. 
1 Manuscript received August 19, 1949. 


Contribution from the Department of Physics, University of Saskatchewan, Saskatoon, 


Saskatchewan. : 
2 Holder of a Studentship under the National Research Council of Canada, 1948-49. 
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In the control circuit to be described here, any desired energy may be chosen 
by selecting a value on a calibrated variable bias control. When the signal 
voltage reaches this value, a trigger pulse is developed in the succeeding cir- 
cuits which causes the electron orbit to be expanded with a minimum time 
delay. 

As mentioned above, the over-all method of integration-expansion is that 
developed at the University of Illinois, the integrator stack and expander 
circults being almost identical. The bias, amplifier, and trigger circuits we 
have designed ourselves with the object of obtaining maximum stability, 
dependability, and accuracy. 


The Integrator 
The integrating network (2, 5) (see Fig. 1) consists of R = 10° ohms (twenty 
50,000 ohm, 100-w. resistors) and C = 0.1 uf. The voltage developed across 


BETATRON INTEGRATOR 
TANK CIRCUIT 


AND 
AMPLIFIER 


EXPANDER 
CONTROL 


EXPANDER 
COILS 


Fic. 1. Block diagram of integrator-expander. 


the integrating capacitor is 86 v. at 23.1 Mev. Owing to the large impedance 
ratio of the integrator components, the 180 cycle magnet voltage is phase 
shifted almost 90 degrees, lacking only about 31 minutes of angle, and the 
phase shift is relatively independent of temperature or aging effects in the R 
and C values of the integrator. Actually the phase difference between the 
magnet current and the integrated voltage is even closer to zero degrees 
because of the effective resistance of the magnet coils. If the time constants 
of the coils and the integrator stack were equal, the signal voltage on the inte- 
grating capacitor would be exactly in phase with the magnet current. An esti- 
mation of the effective resistance of the magnet coils results in a calculated 
phase angle of 23 minutes. Thus the net phase advance of the integrator 
voltage ahead of the magnet current is of the order of 8 minutes of angle, which 
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corresponds to a time difference of about 2.1 usec. in the crossing of the zero 
axis. However, it is advantageous to retain this phase advance in order to 
compensate partially for time lags in the circuits that follow. 


The Bias Circuit and Integrator Amplifier 


The first stage of the amplifier (Fig. 2) consists of a cathode follower used 
as an impedance transformer between the integrator stack and the grid of the 
following tube. To prevent grid—cathode voltage drift, the plate supply and 
heaters are stabilized. 


The sum of the integrator signal voltage, the cathode follower grid bias, and 
the voltage on the tap of the Energy Control helipot, appear across the cathode 
resistor of the follower. Whenever this sum exceeds a fixed value of approx- 
imately 105 v., the 9006-diode unclamps and allows the grid of the first ampli- 
fier stage to be driven positive. The rate of rise of the input signal at the instant 
of unclamping is amplified by three stages to produce a pulse rising at least 
10 v. per usec. when the trigger point is near the axis of an 85 v. sine wave 
input. The first amplifier stage is a .cathode-coupled stage for maximum sta- 
bility; the following two stages are single-ended pulse amplifiers with their grid 
biases adjusted to obtain a sharply rising output. All heater voltages are stab- 
ilized and all tube voltages are VR-tube stabilized. 


The 50,000-ohm Energy Control potentiometer is a 10 turn helipot with 
linearity of 0.1%. A stabilized reference voltage of approximately 105 v. is 
applied across this potentiometer. An important feature of the bias circuit is 
that the selected trigger point is independent of the characteristic curve of the 
clamping diode or nonlinearities of the cathode follower input-output char- 
acteristic. At the instant of unclamping, the cathode follower is always at the 
same point of its operating curve (approximately 105 v. across the cathode 
resistor), the diode current reduces to zero, and the current drawn from the 
reference supply returns to its ‘‘zero’”’ value. Hence, assuming that no drifts 
occur in tube emission, etc., the linearity of the bias circuit depends only on 
the accuracy of the control potentiometer and the reference voltage. 


The Trigger Circuit 


The square wave from the integrator amplifier is differentiated and the posi- 
tive pulses used to trigger a cathode-coupled single-shot multivibrator which 
supplies the correct trigger pulse, through a cathode follower driver, to the 
expander. The pulse delivered by the amplifier has a rate of rise sufficient to 
reduce jitter in the multivibrator to about 0.1 ysec. The plate load in the 
output side of the multivibrator is a 13 mh. air-core inductance. When the 
75 ma. current flowing through this inductance is cut off, the inductance reson- 
ates with the stray capacity of the coil, the diode-clipper filament winding, and 
the wiring. The 6AL5 diode-clipper removes the oscillation after the first half- 
cycle. The resulting half sine wave pulse has 350 v. amplitude, 7 usec. duration, 
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and a rate of rise of 250 v. per usec. A cathode follower transforms this into a 
low impedance pulse to drive 40 ft. of coaxial cable and the grid of the 5C22 
hydrogen thyratron. 


The Expander 


The expander circuits are similar to the University of Illinois’ with only a 
few changes of component values. When the 5C22 hydrogen thyratron is trig- 
gered, a 0.2 wf. capacitor which is charged to approximately 3000 v. is dis- 
charged through a 4:1 pulse transformer. The load on the pulse transformer is 
mainly the inductance of the expander coils. The expander coils consist of one 
turn (10 strands of No. 18 wire) on each pole face, at a radius of 15.9 cm. 


The voltage on the 0.2 uf. storage capacitor is obtained from a resonant 
charging choke and a d-c. voltage supply capable of delivering 200 ma. of 
direct current at 1300 v. A variac in the primary of the high voltage trans- 
former controls the output. 


A wave trap, resonant at 180 cycles, blocks induced voltages in the expander 
coils from sending large currents through the secondary of the pulse trans- 
former. The RL-combination (Fig. 2) across the secondary damps out any 
oscillations that may be set up in the wave trap by shock-excitation. High 
frequency oscillations superimposed on the expander pulse are by-passed by 
the RC shunt across the secondary. 


The combination of storage capacitor, pulse transformer, and expander coils 


produces a half-sine-wave current pulse in the expander coils. The pulse has 
a duration of 20 usec. and reaches a peak amplitude of 200 to 350 amp. in 
8 usec., the actual amplitude depending on the power supply variac setting. 


Circuit Performance and Calibration 


Since the integrator amplifier is direct coupled, the jitter in the unclamping 
instant in terms of input voltage was readily measured and found to be between 
0.1 and 0.2 v. Also the ‘‘zero trigger’”’ point, or zero energy point on the control 
dial, may be easily determined. This setting has been found to drift less than 
0.1 v. after an adequate warm-up period. 


The cathode follower driver for the 5C22 has an output impedance of 160 
ohms and delivers a 310 v. pulse with a rise time of 250 v. per usec. to the grid 
of the thyratron. This trigger pulse will reduce the average time jitter in the 
firing of the 5C22 to the order of 0.01 usec. (4). Also, the time delay between 
the triggering of the multivibrator and the 5C22 will be of the order of 0.4 usec. 
Owing to the finite rise time of the expander pulse a further variable delay of 
from 3 to 6 usec. is added, since expansion occurs at different expander currents 
for different energy settings. 


A discussion of some of these errors in energy calibration appeared in a recent 
paper by McElhinney et al. (5). 





118 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


A summary of estimated and measured errors in the linear energy scale of 
our circuits is given in Table I. Corrections which would have to be made for 
an absolute energy calibration (2) are given approximately in Table II. The 
correction for the fixed delay of 2.4 usec. would be at all times less than 0.05 
Mev. Actually this correction is very nearly compensated by the 2.1 usec. 
phase advance of the integrator signal. There is also a variable delay present 
owing to the finite rise time of the expander pulse. In addition, an energy 
change is produced by the flux from the expander pulse. Since both the variable 
time delay and the expander pulse effect are approximately proportional to the 
initial magnetic field at the time of expansion (5), only small nonlinear correc- 
tions to the energy scale should be necessary. McElhinney (5) estimates that 
the nonlinear corrections from the expander pulse effect are of the order of 


0.02 Mev. 
TABLE I 
SUMMARY OF ERRORS IN LINEAR ENERGY SCALE 


Source of error 





Jitter in: 
amplifier 
multivibrator 


5C22 





Total 
Nonlinearity in expander pulse 
correction 





Nonlinearities due to delays in | 
circuits | 0.03 


Error in linear energy scale 0.1 Mev. 


TABLE II 
SUMMARY OF INTEGRATOR EXPANDER CORRECTIONS FOR AN ELECTROMAGNETIC CALIBRATION 





Source of correction Correction 


Correction, 


Fixed delay, Variable delay, 
Mev. 


sec. msec. 
Amplifier 2 
Multivibrator —0 
5C22 0.4 
Expander — 3 to 6 


| 
ia 
| 
| 
| 
| 
| 
| 
| 


+0.05 

| | A+3to+6 |“ +.08 to +.15 
Integrator phase advance —2.1 —0.05 
Expander pulse | “—.1to —.2 


Totals | +9.4 


An absolute electromagnetic calibration (2) of the integrator-expander is 
possible. However, in view of the number of small but somewhat uncertain 
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corrections which must be made, it was thought preferable to assume a linear 
calibration and to use the measured values (3, 5) of the thresholds of several 
of the better-known nuclear reactions to establish points on this curve. As 
indicated by Table I, the maximum deviations from a linear energy scale due 
to jitter and nonlinear corrections should be less than 0.1 Mev. The reactions 
selected for our calibration were C!(yn)C" with threshold of 18.7 Mev., Cu® 
(yn)Cu® at 10.9 Mev., and Be*(yn) Be’ at 1.67 Mev. (5). 
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A MORE EXACT FRESNEL FIELD DIFFRACTION RELATION! 


By G. A. WoontTon 


Abstract 


In general, diffraction relations which were developed for optical purposes 
have been found to be useful in the prediction of electromagnetic fields due to 
radiating apertures but, because many radio-optical measurements must be made 
so close to the aperture that the mathematical approximations in the Fresnel 
field relation often are invalid, new relations, free from these approximations, 
must be developed. ‘The optical, Fraunhofer relation has been found to predict 
the distant field with good accuracy and for this reason the Fresnel field is cal- 
culated from the distant field equation by the Fourier transform method. The 
calculation is found to result in an integral which reduces to Fresnel’s integral 
when the standard optical approximations are made in it. The integral has not 
been evaluated. 


Introduction 


Two types of approximation are present in the calculation of the electro- 
magnetic field in front of a radiating aperture by optical diffraction methods. 
These are approximations in the fundamental physics of the solution and 
approximations in the algebraic and trigonometrical processes by which 
Kirchhoff’s general solution is reduced to the solution of a specific problem. 
The writer, with others (4, 5), has found that for a square aperture as small 
as 10 wave lengths wide, the field can be predicted with good accuracy by 
optical diffraction formulae, provided the point at which the field is measured 
does not lie at an angle of more than 30 or 40 degrees from the normal to the 
aperture, and provided that the radius vector from the center of the aperture to 
that point is greater than three or four times the aperture width. From these 
results it can be inferred that the Kirchhoff solution of the diffraction problem, 
in spite of its inconsistencies, can be used with success to predict the field of 
practical radiators, but that trigonometrical and algebraic approximations 
make the specialized statements of that solution inadequate for use in the 
region close to the aperture. 


The field close to the radiator is particularly interesting in radio-optical 
measurements. Measuring instruments, such as probes and lenses, are usually 
placed in the near field, and it is only by a comparison of the predicted with 
the measured field that the experimenter can estimate the disturbance caused 
by the measuring instruments. Obstructions in the near field modify distant 
patterns; an exact knowledge of the Fresnel field is required to predict the 
effect of the obstruction on the distant field. Apart from these practical con- 
siderations, it will not be possible to press the investigation of the inadequacies 
of the Kirchhoff solution much farther until the near field can be predicted 
from Kirchhoff’s work by equations that are free from extraneous errors. For 
these and related reasons, optical diffraction formulae, free from mathematical 


1 Manuscript received October 28, 1949. 
Contribution from the Department of Physics, McGill University, Montreal, Que. 
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approximations, are needed for the prediction of the Fresnel field of radiating 
systems. The analysis which follows applies only to radiating apertures of the 
linear type, that is, to those with a rectangular aperture and for which the 
aperture distribution of electric field or potential, F(xy), can be expressed as 
the product of separate distributions in the XY planes, Fi(x) Fu(y). The 
analysis leads to an integral statement of the Fresnel field and this integral, 
like Fresnel’s integral (to which it can be reduced by approximation) has not 
been integrated, yet, into a closed form by any of the standard methods. 


The Fresnel Field 


Experimental evidence shows that standard optical formulae can predict 
the distant field of a linear radiator but fail'in the near field, and, for this reason, 
it is proposed to calculate the Fresnel field from the distant field in order to 
take advantage of the relative freedom of the distant field calculation from 
mathematical approximations. 


Define a right-handed rectangular co-ordinate system with its origin at the 
intersection of the diagonals of a rectangular aperture, its X and Y axes per- 
pendicular to the two sides. Let the axes be oriented so that the Z axis projects 
into the space where the field is to be predicted. Let the distribution of poten- 
tial in the aperture be given by the product F,(x) Fii(y). At z = b, along the 
Z axis, set up a subsidiary co-ordinate system (£, 7) with the — and 7 axes 
parallel to the X and Y axes respectively. The aperture distribution F(x) 


F\,(y) will set up a Fresnel field in the plane of this new co-ordinate system 
which will be described as F_(¢) F22(n) (although in fact each of these must be 
a function of } and X as well as — or n). The aperture, because of its electric 
field, will radiate a spherical wave which, by ordinary optical Fraunhofer for- 
mulae, can be expressed to close approximation as: 


eee ay 


o(xWit1 = — 
Ri 


+00 


. W a [rn (x) exp ES — | dx \ 


+a 


x ‘ V Z | Fii(y) exp | i2ny stu dy \ ’ (1) 


where $(x19121) may be interpreted as the electric potential at Pi(xyiz), w 
and X have the significance of the angular velocity and the wave length of the 
radiated field, 7 is the symbol for V — i, Riis the radius vector from the origin 
to P,, 6; is the angle between the Z axis and the projection of the radius vector 
on the XZ plane, and @;; has the same meaning in the YZ plane. [Cos 6 + 1] 
is the inclination factor, and in it @ is defined in the usual manner. 
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Because the aperture field is not unique, it would have been equally correct 
to say that the distant field was caused by the Fresnel distribution F2() F22(n) 
in the £,7 plane. In that case, one would say that a spherical wave was radiated 
from the point (0, 0, b) and an equation similar to (1) would be written in 
which, however, would appear the functions F2,(£) and F22(n) and new polar 
co-ordinates Re, 2, and 422, all appropriate to the new length and angle of the 


radius vector from z = b to P,. 


In order to consider the relation between the two spherical wave fronts, set 
R,i= R2+ b and, because the expression can be separated into independent 
functions of the co-ordinates, restrict the attention to the two circular arcs 
lying in the XZ plane which represent lines along the two wave fronts radiating 
respectively from the origin and from (0, 0, b). These arcs are of radius R; and 
R.. Draw the radius vector R; to some point P; on the circular arc related to 
the wave from the origin. It will cut the arc related to the wave front from 
(0, 0, b) in another point P:. The distance along the radius vector between 
the two wave fronts P, and P; can be expressed as 


5 = b(1 — cos @) (2) 
and this expression can be made as accurate as desired by taking R; and R2 


large enough, compared to 0; the same reasoning leads to the use of common 
angle @ in (2) since for R; and Rez very large 6;= 02.= 0. 


The two wave fronts, because of their separation, 5, along the radius vector, 
are out of phase; the field at P» leads the field at P; by a factor: 


exp| 52%] exp| 37% b(1 — cos 0) | 


oo Jeol -f()-(2E}} 


No inaccuracy will be introduced by taking the modulus of the field to be the 
same at P, as at P, because of the great magnitudes of both R; and R: com- 
pared to b and the consequent insignificance of 6. This information is sufficient 
to fix the relation between the two integral factors: 


+00 


| F,(x) exp E WX sin 4] dx 
—co 


which appeared in Equation (1), and 


+00 
c ‘ 
| PO exp ES ant] dé 


which would have appeared if the equation for the wave radiating from (0, 
0, b) had been written. Rewrite (1) so that it expresses the field along the arc 
on which lies P2, that is, the arc associated with the wave radiating from 
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(0, 0, b). This can be done by advancing the phase of all points P, to the 
phase of points P2, using (3), by writing R: = Re +5 in the phase factor exp 


; 2 ; ’ ; ; ‘ 
| ai(w ~ **) and by interchanging R, with R:2 in the denominator. Equate 
the rewritten Equation (1) to the Fraunhofer equation (similar to 1) that could 


have been written for a wave radiating from (0,0, 5). After cancelling common 
factors, the equation appears: 


saoartee 
: * 9 * si 2 4 , _* 6 
Vy; exp| - 4 (2) —- (22°) f |Jn (x) exp ES sin) a 


+c 
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The two factors Y; and ¥2 are: 


+c 


= a/ dt [ yy exp | jony i Lay: 
n=4 J | Fao) exp | 2m x Ja: 


= 4/1 | a(n) exp | jon Sito 
Y= V x | Fx (n) exp [2mm xX |e 


By restricting this discussion to two circular arcs in the XZ plane, the angles 
61, and 622 have been fixed at zero and further, since Y; and Y2 in general 
describe the intensity of two spherical waves along a pair of arcs in the YZ 
plane, in the special case under consideration they describe the intensity at 
the same point and are, therefore; the same constant; consequently Y¥, and Y2 
can be cancelled from Equation (4). 


Michelson (2) and other authors have pointed out the reciprocal relation 
between the aperture distribution and the distant field that it causes, and have 
justified the calculation of one from the other by the Fourier transform pair 
relation. Following these authorities, one may regard Equation (4) as a rela- 
tion between Fourier transforms: hence define 

sin 6 3 
5 
x (5) 
+o 
Gi(y) = [re exp | — j2ryx]dx 
—co 


+a 


r 


-o 


= exp | - 54 (Bey = (2m) | 





124 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


so that one can rewrite Equation (4) in the form 
G,(y) = Gs(v)Gi(7). (9) 


From the reciprocal property of the transformation it follows that the field 
F.(é) (the prediction of which is the objective of this analysis) is given by: 
+00 


F,(é) = Jeune exp [+ j2ryé] dy 


—-co 
+c 


= foun exp| - b ‘ (fy 2x} exp [+ j2ryé] dy. (10) 


Equation (10) can be transformed to an expression in terms of the aperture 
field, Fi(x), by direct integration or by reference to the tables of Campbell and 
Foster (1) (see pairs’*202 and 867). The transformed integral is 


- [2m f -)2 VP) ae 
qa Slat eE 
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F.(é) = | Fi(x) ( 


In this integral K,(é, x, A, 0). is a modified Bessel function of the second kind 
and first order. It can be written as a Hankel function through the relation 


K,(2) = > nj HH) (ja) = ; nj" H® (7712). (12) 


and in that form can be separated into real and imaginary parts in terms of 
the Bessel functions of the first and second kinds. 


Equation (11) is very like the simple equation for the Fresnel field in terms 
of Fresnel’s integral. The Bessel function is a phase factor of a somewhat more 
obscure kind than the simpler exponential function, while the denominator is 
the distance from an aperture point x to a field point &. The whole expression 
directs the user to find the potential at a point — by adding up, in a particular 
way, elements F;(x)dx in the aperture, having due regard for the phase and 
amplitude changes suffered in transmission from x to &. Had the problem been 
solved for the surface of a sphere, instead of, for the arc of a circle, a 
function F(, 7) would have been predicted at a point in the &, 7 plane as the 
product of two integrals each similar to that in (11) except for the appropriate 
change in the variables. The real and imaginary parts into which it can be 
split are reminiscent of Fresnel’s C and S functions. Because the Bessel func- 
tions, of which the real and imaginary parts are integrals, oscillate, it appears 
that the integral (11) will oscillate in phase and amplitude and that in fact, if 
the integrals were known, a curve, which would bear a relation to Cornu’s 
spiral could be constructed to predict these quantities. 


The standard development of the Fresnel field in terms of the Fresnel’s 
integral makes use of the following approximations: it is assumed that distances 
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between any aperture point x and any field point can be taken as } in the 
calculation of the modulus of the expression but that in the phase term a 
second approximation is necessary. The second approximation is obtained 
from a binomial expansion of the correct expression for the distance between 
x and é: 


r= 1° +(¢-ayt' 0 41 G- (13) 


so that the phase term is written: 


en ~ tee ten ae te ae 
exp| x Jeo bx (é “| (14) 


This development is given in detail by Slater and Frank (3, pp. 180-187). 
Equation (11) can be verified by showing that with the same type of approxi- 
mation, it reduces to the Fresnel’s integral statement of the field. Assume that 
b is large compared to X, so that the argument of the Bessel function must be 
large and the approximation: 

Eis | 2 exp (—2) (15) 
s> 2 VV TZ 


co 
is justified. That is: 
pA ca oan 2 
Kil 9 = (€ — x)? - €X 
. 2xb dl ; 2 \ 
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A wi(—é—x)P?+ Bh} (E — x) + *} 
(16) 
Assume, next, that as in (13), } is large when compared with (¢ — x) so that 
in the exponential function of (16) the argument can be approximated by two 
terms of the binomial expansion and in the denominator (~ — x)? can be 
neglected in comparison to b®. With these further approximations (11) becomes: 
+00 


F,(é) = exp| - 572 | fi ze exp| - i (§ — »| dx. (17) 


—o 


This approximate form is the same as that which is ordinarily obtained from 


the simple theory. The factor 4/2. is {correct in this case since it must be 


multiplied by another similar factor, which is carried by the function giving 
the variation with y, when the complete expression is written. Before Equation 
(1) was written, the expression was approximated in the usual way by removing 
the inclination factor from under the integral sign; this approximation per- 
mitted the inclination factor to disappear by cancellation in Equation (4). It 
should now be brought back and with it the factor 4 by multiplying the right 
i = ; = 1 + cos @ gh 
hand sides of Equations (10), (11), and (17) by ers ) outside integral 
sign, to give equations which are as correct as can be obtained, without intro- 
ducing new complexities. 
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EFFECT OF THE ILLUMINATION OF THE NEGATIVE ON THE 
RESOLVING POWER OF A PROJECTION PRINT 


By P. A. TATE 


Abstract 


Axial measurements have been made of the resolving power of photographic 
projection printing systems. Besides the usual diffuse and specular systems a 
third system was considered in which the aperture of the projection lens is not 
entirely filled by the image of the source. If the relative aperture of this system 
is calculated considering the source image as the entrance pupil of the lens, then 
at small relative apertures this system has greater resolving power than the other 
systems. However, there is no advantage in using the third system at small 
relative apertures. At large to moderate relative apertures there are only very 
slight differences among the three systems. Diffusion screens placed in a con- 
denser enlarging system may cause considerable reduction in resolving power 
if the lens was not previously filled by the image of the source. This effect, and 
also the high contrast that is very frequently present, are thought to be the origin 
of references in the literature to sharper pictures obtainable with condenser 
enlarging systems. 


Introduction 


There are two extreme types of enlarging systems that may be termed 
‘diffuse’ and “‘specular’’. In the diffuse enlarging system the negative is illum- 
inated by a large diffuse source emitting light in all directions. In the specular 
‘system the light is directed exclusively towards the lens through the negative 


by some optical device. An example of each of these types of illumination is 
shown in Fig. 1. A great many enlargers employ systems which are neither 
completely diffuse nor completely specular. 


The properties of diffuse and specular systems differ in several important 
respects. In the diffuse system the lens is always completely flooded with light. 
The effective printing density of the negative is approximately diffuse den- 
sity (5). The amount of light available at the printing plane is relatively low 
owing to the large amount of light which falls outside of the lens and to absorp- 
tion of light in the diffusing screen. 


The lens of a specular enlarging system may not be completely flooded with 
light; e.g., in a condenser enlarger if the source image is smaller than the en- 
trance pupil of the projection lens. Specular density is effective in determining 
the tones of the image in the printing plane (5). Specular systems can be made 
much more efficient than diffuse systems from the standpoint of light available 
on the photographic paper. 


Tuttle (6) has studied the problems arising from the differences in behavior 
of specular and diffuse systems, particularly regarding tone reproduction. He 
has discussed the change in contrast that accompanies a change in the diffusion 
characteristics of an enlarger and compared the effective speeds of specular and 

1 Manuscript received November 1, 1949. 
Contribution from the Division of Physics, National Research Laboratories, Ottawa, 
Canada. Issued as N.R.C. No. 2077. 
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diffuse systems. Information is also given on the effectiveness of diffusion to 
reduce graininess and to minimize contrast of scratch images in the final print. 


It has been mentioned in the literature (2, p. 430) that specular enlarging 
systems give sharper pictures than are obtainable with diffuse systems; it is 
not clear whether the comparisons have been made under like conditions of 
effective relative aperture and print contrast for each type of system. Recently 
this question has become of greater interest with the appearance on the market 
of lamps with very small and uniformly bright luminous areas. Claims have 
been made (1) for these lamps that their use in optical systems can produce 
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Fic. 1. Two types of negative illumination. 


results of great resolution and depth of focus. Any exact information on the 
relative merits of diffuse and different kinds of specular enlarging systems is 
therefore distinctly useful. 


Practical condenser enlargers usually employ tungsten filament lamps in 
clear bulbs as the light source. Most enlargers with large frosted bulbs and 
condensing lenses are closer to diffuse systems than truly specular systems. 
Because of the irregular shape of a clear tungsten source a circular lens aperture 
can never be completely filled by the source image. We see then that enlargers 
actually in use are of three types, or combinations of three types: ordinary 
diffuse, ordinary specular and specular in which the image of the source covers 
only a portion of the projection lens aperture. 


In the last system the resolution and depth of focus depend on the distri- 
bution of light in the projection lens. We shall use the terms “‘relative aperture 
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of the diaphragm’”’ and “‘relative aperture of the source image”’ to describe the 
illumination of the projection lens in discussions of this system. The relative 
aperture of the diaphragm is the normal meaning of the term ‘relative aper- 
ture’, i.e., the ratio of the diameter of entrance pupil of the projection lens 
to the focal length of the lens. The relative aperture of the source image will 
be defined as equal to the relative aperture of the diaphragm when the dia- 
phragm just surrounds the source image. The corresponding terms ‘‘aperture 
of the diaphragm”’ and ‘‘aperture of the source image”’ will also be used. Pro- 
vided that the source is a uniform disk the relative aperture of the source image 
largely determines the illumination in the printing plane. If there were no light 
outside the source image then the relative aperture of the source image would 
be the only quantity needed to estimate the performance of a system of this 
type. However, this condition can never be completely satisfied because there 
is always some light scattered outside the source image from imperfect con- 
denser surfaces and from the negative. Thus, both the relative aperture of the 
diaphragm and the relative aperture of the source image are required to describe 
a condenser enlarging system with an unfilled lens aperture. 


Condenser enlargers are frequently used with a certain amount of diffusion 
introduced into the system. Reports are sometimes seen that this practice 
decreases the resolution of the condenser system. It was thought desirable to 
find the extent of the change of resolution in a number of representative lenses 
when diffusing screens are introduced into a condenser system in which an 


image of the source smaller than the aperture is formed in the projection lens. 
This has been done and is described in Section II of this paper. 


I. Comparison of Enlarging Systems 
EXPERIMENTAL METHODS 


The axial resolving powers of three types of enlarging systems were compared 
under controlled conditions. The comparison was made for one diffuse system 
and two types of specular systems. 

The three systems were: 

(A) A diffusely illuminated system. 

(B) Acondenser system with a uniform diffusely emitting source imaged 

on the diaphragm of the enlarging lens. The source was just large 
enough to cover the lens aperture at its widest opening. Stopping 
down was provided by the projection lens iris. 
A condenser system with a uniform diffuse source imaged in the 
enlarging lens as in (B). The lens was kept at its maximum dia- 
phragm opening and stopping down was accomplished by reducing 
the size of the source and therefore of itsimage. System (C) of a given 
relative aperture of the source image was compared with Systems 
(A) and (B) at an equal relative aperture determined, of course, by 
the diaphragm. 
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The enlarging apparatus was mounted horizontally as shown in the photo- 
graph Fig. 2. A diffusely emitting light source was made from a dense piece 
of opal g'ass, illuminated from behind, and masked off to the right size with a 
sheet of thin brass. The size of the hole in the brass mask was selected to give an 
image of the required size at the projection lens. A gunsight lens of high quality 


Fic. 2. Photograph of enlarging apparatus 


was used as a condenser lens to provide this image. The effect of complete 
diffuse illumination was obtained by putting a piece of dense flashed opal glass 
between the condenser lens and the test object. Adjustment was provided for 
squaring-on the lenses, negative holder, and paper holder. Accurately known 
displacements of the paper holder could be made with a micrometer advancing 
screw for focusing. 


Photographic Details 


In the position of the negative a high contrast, high resolution bright annulus 
type (3) resolving power test object was placed at the center of the field. Suc- 


TABLE I 


Step | Lines per millimeter 


18 8. 
19 9. 
20 se: 
21 B 
22 
23 13. 
24 
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cessive steps on the test object increased in size in the nominal ratio 1: -¥/2. 
The resolving power test object was surrounded by a step tablet of the same 
emulsion for control purposes. The degree of magnification of all enlarging 
systems was kept at 1.8X. Fig. 3 shows a sample enlarged print of the test 
object. Prints were made on Kodabromide F1 paper and were unferrotyped. 


Table I gives the calibration of the resolving power graticule in lines per 
millimeter on the enlarged print as defined in reference (3). 
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Fic. 3. Enlarged print of the resolving power target and step tablet. 


Development was for 90 sec. in D-72 diluted 1:4. Brush development was 
employed with other precautions usual in sensitometry. 


Prints made with all types of enlarging systems had the same contrast. It 
was possible to achieve this because a very fine grain, and hence low scattering, 
emulsion (Kodak 548-G) was used for the test object, and the density difference 
between test object annuli and the background was large enough to take in 
practically the whole scale of the paper. The condition of constant final con- 
trast is essential to valid comparison of enlarging systems. 


Results and Discussion 


In Fig. 4 the curves compare (A) diffuse illumination, (B) condenser illum- 
ination with a flooded aperture, and (C) condenser illumination with a small 
image in the center of the open diaphragm. Systems (A) and (B) having a 
given relative aperture were compared with System (C) having an equal rela- 
tive aperture of the source image. The aperture of the diaphragm in System 
(C) was kept constant at F/4.5 for these curves. All resolution figures are those 
of the focal position determined by maximum axial resolution. No significant 
change of focus with stopping down or diffusion was found for any lens. The 
resolving power is given in steps resolved on the sixth-root-of-two test object 
and is therefore on a logarithmic scale. Relative apertures are also plotted on 
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a logarithmic scale. No significant differences appear among the three systems 
down to an aperture of about F/16 although System (C) seems to have some- 
what lower resolution than the others. Beyond an aperture of the source image 
of F/16 System (C) seems to be considerably superior to the other systems to 
the extent of three steps or more. 
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Fic. 4. Axial resolving powers of three different enlarging systems compared. 


The failure of the resolving power to drop at small relative apertures of the 
source image in System (C) is unexpected because, except for the absence of a 
physical diaphragm around the source image, System (C) is optically similar 
to System (B) in which a decrease in resolving power takes place. The difference 
between (B) and (C) can arise only from the small amount of extra light ad- 
mitted by the open diaphragm of case (C). This light was found to be about 
10 or 15% of the total light of System (B). A series of resolving power deter- 
minations were made with a lens used in System (C) in which the relative 
aperture of the diaphragm was reduced in successive steps, thus progressively 
shutting off the extra light around the source image. 


Fig. 5 shows the result of closing down the diaphragm of a system of Type 
(C) having an aperture of the source image of F/32 to produce finally one of 
Type (B). Resolving power is noticeably reduced at apertures of the dia- 
phragm of F/16 and smaller. The change is everywhere gradual and investi- 
gation showed no sudden changes as the aberrational fringe about the source 
image was cut off. 


An explanation for the behavior of System (C) can be given. It is probably 
true to say that the image in this system consists of an imperfect image resulting 


* 
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from the diffraction effects of a small aperture combined with a fainter, more 
perfect image formed by light from the condenser and graticule surfaces passing 
through the regions of the lens aperture at which the lens develops its greatest 
resolving power. Because the test object was of high contrast, reduction in 
image contrast affected resolution only slightly (4, p. 850) and the final com- 
bined image assumed for System (C) was representative of the good image 
rather than the poor one. A slightly smaller resolving power for System (C) 
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Fic. 5. Variation of resolution as diaphragm in System (C) is closed down to smaller ‘relative 
apertures. Aperture of the source image is F/32 throughout. 
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Fic. 6. Resolution through the focal range for Systems (B) and (C). 


than (A) or (B) at moderate relative apertures of the source image may also 
be explained on the basis of the two types of images. 


Further confirmation of the above view was had from a comparison of the 
curves of resolution through the focal range in Systems (B) and (C) shown in 
Fig. 6. The curves for Systems (B) and (C) at F/32 differ markedly. Depth of 
focus in System (B) is very great and the curve is typical of normal systems 
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of this type at small relative aperture. On the other hand for System (C) the 
curve has a much sharper peak showing that the cone of light affecting reso- 
lution is broader in this case. It corresponds, in fact, more closely to the reso- 
lution curve for System (B) at F/16 which is included in the figure for com- 
parison. 


The depth of focus found for System (C) is not equivalent to that usually 
obtained with systems of such small relative aperture and in extreme out of 
focus positions the resolving power drops below that of System (B) at the same 
relative aperture and focal position. No greater depth of focus therefore may 
be had by using small sources in the System (C) than by using System (B) at 
larger apertures. 


II. Decrease in Resolving Power 
EXPERIMENTAL METHODS 


Experiments were made to determine the extent of the decrease in axial 
resolving power that takes place when diffusing screens are introduced into a 
condenser system in which the projection lens is not completely filled by the 
source. The experimental system reproduced that of a practical enlarger in 
which the projection lens was rather greatly stopped down by the size of the 
source image. 


Apparatus very much the same as that already described was used. The light 


source in this case was either a 3 c.p. or a 32 c.p. clear lamp, depending on the 
time of exposure required. These two lamps had filaments of very nearly the 
same maximum dimensions and gave the same results in similar systems. Each 
source image was adjusted so that its maximum dimension covered a portion 
of the lens corresponding to an aperture of F/16, although the lens was set at 
its largest diaphragm opening. Thus the system was of Type (C) defined in the 
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first part of this paper but differed from the system previously studied in two 
respects. The source was not a uniform diffuse disk as before but a filament of 
irregular shape, and the system was used only at large relative aperture where 
the resolving power effects peculiar to System (C) are not encountered. 


Four different types of diffusive screens were interposed between the con- 
denser and the test object. They were made of glass ground on one and on 
both sides, a medium and a dense flashed opal. Luminance-angle curves as 
seen by a 926 photocell for the four varieties are shown in Fig. 7. To each was 
assigned a degree of diffusion, ‘‘diffusion efficiency” (7), defined as the ratio 
of the area under the luminance-angle curve to the area under the corresponding 
Lambert’s law curve. There are obvious difficulties to any attempt to describe 
the properties of a diffuser by one number, but this is easy to calculate and 
gives a good comparison between diffusers of the same type. 


Photographic details were the same as in the previous set of experiments. 


Results and Discussion 


Curves of resolving power taken through the focus for an enlarging lens used 
with different diffusers are shown in Fig. 8. The curves vary but slightly in 
shape with different degrees of diffusion and except for zero diffusion (specular 
illumination) are very similar. No change in focus (maximum resolution) can 
be noticed in these curves. While the peak of resolution in the diffuse case is 
a little broader than for specular illumination, the latter gives good resolution 
over a larger range, and from this point of view may be said to have greater 


depth of focus. 
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In Fig. 9 the maximum axial resolving power for four lenses is shown plotted 
against the diffusion efficiency of ihe diffusing screen behind the test object. 
The curves all drop rapidly at low values of diffusion and reach quickly a 
constant resolving power. This may be explained by the fact that initially only 
the central portion of the lens is working and resolution drops as light is dif- 
fused over the outer and less well corrected portions of the lens aperture. After 
a considerable amount of light has been diffused over the lens aperture a further 
increase in diffusion does not decrease resolution any more from this cause but 
only contributes to light producing veiling glare. 
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Fic, 9. Axial resolving power vs. diffusion efficiency for an unfilled aperture. 


Conclusions 


Provided that the contrast of the final print is kept constant there is no 
significant difference in resolving power obtainable at large or moderate aper- 
tures for the three different types of enlarging systems studied. This result was 
obtained with a high contrast target, but the result applies with even more 
strength for low contrast detail. It should be noted that it may not always be 
possible to achieve the same print contrast throughout the tone range in dif- 
ferent systems. At very small relative apertures of the source image the System 
(C) is distinctly superior in resolving power to Systems (A) and (B) used at 
equal relative apertures. The value of its superiority is nullified by the fact 
that at least as great resolving power may be had with Systems (A) and (B) 
working at larger relative apertures, and also because no increase in depth of 
focus is achieved with System (C) over that of the same lens used conven- 
tionally at moderate relative aperture. 


Small amounts of diffusion applied to a condenser enlarger system result in 
a rapid decrease in resolving power if the full aperture of the enlarging lens is 
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not already utilized. Probably the extra sharpness achieved with pure con- 
denser systems has been due to the stopping-down effect of an unfilled lens 
aperture and also the increased contrast that may be had with these systems. 
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ISOTOPE SHIFT IN THE RESONANCE LINES OF ZINC! 


By M. F. Crawrorp, W. M. Gray,? F. M. KELLy,? and A. L. SCHAWLOW 


Abstract 


The isotopic structures of the Zn I resonance lines 2139A (4s?'S)—4s4p'P1) 
and 3076A (4s?!S)—4s4p3P1) excited in an atomic beam source have been 
resolved with a Fabry-Perot etalon. Shifts in cm.—' relative to the Zn* com- 
ponent are +0.016 and +0.033 in 2139A, and +0.023 and +-0.046 in 3076A for 
Zn® and Zn®*, respectively. | When the normal mass shiits are subtracted, the 
residual shifts for Zn® are +0.010 (2139A) and +0.030 cm.— (3076A). The 
field effect theory of isotope shift predicts equal shifts and the specific mass effect 
predicts unequal shifts for the 'P and *P levels. Therefore the difference of these 
residuals is a specific mass effect. The difference has the sense predicted by the 
specific mass theory and its value, 0.020 cm.—"', is a large fraction of the observed 
shifts. Thus the specific mass effect must be taken into account before the field 
effect theory can be used to obtain nuclear properties from isotope shifts in 
intermediate elements. 


Introduction 


Two causes of isotope shift in atomic spectra have been recognized, the mass 
effect and the field effect. The mass effect consists of two parts, the normal and 
the specific. The normal effect, which is present in all spectra, is easily calcu- 
lable and decreases, rapidly at first, as the atomic weight increases. The specific 
effect, which was first treated by Hughes and Eckart (14), is present in spectra 
of atoms with more than one electron. It is much more difficult to evaluate 
precisely than the normal effect. 


The contributions of the specific mass effect to isotope shifts have been 
evaluated theoretically for transitions in the spectra of lithium (14), neon (6), 
magnesium (23), boron (19), and helium (1, 8, 12, 17). For about half the 
transitions considered there is fair agreement with experiment, but for the 
other half the agreement is poor. It is probable that the cases of disagreement 
between theory and experiment in the spectra of neon and magnesium are due 
to the use of wave functions that are not sufficiently precise for an accurate 
evaluation of the integrals. He | affords the best test of the theory since there 
is at most one integral for each configuration and it can be accurately eval- 
uated. For the singlet transitions (12), the discrepancy between the observed 
and predicted shifts is 4% on the average and in no case is it greater than 9%. 
Discrepancies as large as 30% are found for the triplet transitions, but these 
discrepancies may be partly related to the fact that the triplet fine structure 
splittings appear to be different for the two isotopes. Thus the He I results 
confirm the Hughes—Eckart theory as a first approximation, but also indicate 
that the omission of the electron spin interaction in the potential function may 
be partially responsible for the residual discrepancies. 
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The field effect (5, 9, 20, 21) has been used to explain isotope shifts in the 
spectra of the heavier elements. Here the normal mass effect is negligible and 
the field effect gives a consistent explanation (11) of the observed shifts. 


No satisfactory explanation has been advanced for isotope shifts in elements 
of intermediate atomic number. Among the lightest elements there is a general 
decrease in the magnitudes of the specific mass shifts with increasing atomic 
number. This has led to the tacit assumption on the part of some authors that 
the specific mass effect is negligible in all but the lightest elements. Isotope 
shifts observed in Cu I and Zn II have been assigned (22) to levels as if they 
were due to the field effect. But the shifts are in the opposite sense to those 
predicted, if the nucleus of the heavier isotope occupies the greater volume. 
This inconsistency raises the question whether the specific mass effect may be 
important in the intermediate elements, as suggested by Bartlett and Gibbons 
(6). Since the number of integrals contributing to the energy level shifts in- 
creases with increasing atomic number, this may well be the case in spite of the 
decreasing factor (1/M,— 1/M2), where M, and M, are the nuclear masses of 
the two isotopes. 


The specific mass effect in any given transition cannot be estimated without 
an involved calculation. But it is possible in certain cases by comparing the 
shifts of a pair of lines to determine the importance of the specific mass effect 
without a detailed calculation. Such a pair is the resonance lines of zinc, 2139A 
(45?'.S,—4s4p'P,) and 3076A (4s?'S)—4s4p%P!). Since these lines have upper 


levels belonging to the same configuration and a common lower level, they 
should have the same field effect shifts, as the field effect shift of either a 4p1/. 
or a 43/2 electron is very small relative to that of a 4s electron. The specific 
mass theory, however, predicts that the shifts of the lines should differ by an 
amount depending only on C(4p, 4s), one of the integrals contributing to the 
shifts of the upper levels (6). 

Experimental 


The resonance lines of Zn I were excited in an atomic beam source and re- 
solved with a Fabry-Perot interferometer. The crystalline quartz optical flats 
of the interferometer were coated with aluminum by rapid evaporation (10). 
The interferometer was placed between the collimator lens and the prism of a 
fast spectrograph in which an off-axis paraboloidal mirror is used as the camera 
objective. The source and spectrograph will be described in detail in forth- 
coming publications. 


The intercombination resonance line (3076A) was photographed on Ilford 
Ordinary and Eastman 103-0 plates. For the singlet resonance line (2139A) 
Ilford QI plates were used. 

Results 
The line 3076A was resolved without difficulty by a 40 mm. etalon and a 


55 mm. etalon. This line is not strongly absorbed in the atomic beam. The 
density of the beam was therefore made as high as possible without spoiling 





140 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


the collimation by collisions. No interference was experienced from the second 
member of the sharp triplet series (3072A) since the intensity of this triplet in 
the atomic beam source is negligible compared to that of 3076A. 


The structure of 3076A consists of three components, the strongest having 
the lowest frequency and the weakest the highest frequency. The relative 
intensities of the three components were measured by calibrating the plates 
with stepped-slit continuous spectra of a hydrogen lamp. The density con- 
tours were obtained with a Leeds and Northrup microphotometer. The peak 
intensity of each component was determined and corrected by subtracting the 
intensity contributions of the other two components. These corrections were 
small for the 55 mm. pattern in which the components were well separated. 
The measured intensities were in good agreement with the relative abundances 
(18, 15) of Zn*, Zn®, and Zn**. The separations of the components of 3076A 
relative to the Zn™ component are +0.023 +0.001 for Zn® and +0.046 +0.001 
cm. for Zn®8, 

On the more strongly exposed 40 mm. patternsa weak fringe is observed and is 
attributed to Zn® which hasa nuclear moment. The Zn® line should have three 
components with the separations of the hyperfine structure states of 4s4p *P). 
These separations have been determined by Lyshede and Rasmussen (16) 
and confirmed by Arroe (2) from the structures of other lines of Zn I. The 
observed fringe is due to overlapping of adjacent orders of the two outer com- 
ponents of Zn” 


The line 2139A was difficult to resolve because of the strong self-absorption 
in the atomic beam. Broadening of the components and even partial reversal 
were observed when the density of the beam was high. The half-widths of the 
components of this line were made sufficiently small to permit resolution by 
reducing the density of the highly collimated beam. Under this condition 
exposures of three to four hours were necessary, even with thinner aluminum 
films in the etalon. Exposures with a 3 mm. spacer showed that the instru- 
mental width of the fringes corresponded to an effective reflection coefficient 
of 0.6 at 2139A. With these films the structure of the line was just resolved 
using a 73 mm. spacer. 


The wave-number separations of the components of 2139A relative to the 
Zn* component are+0.016 + 0.002 for Zn® and + 0.033 + 0.002 cm.“ forZn®. 
These values include corrections of +0.002 and +0.003, respectively, for the 
overlap of the components. 


Each of the lines 2139A and 3076A has three equally spaced components 
due to Zn*, Zn®, and Zn®*. The structure observed for 3076A disagrees in 
magnitude and in sense with that estimated by Billeter (7) from a study of the 
absorption by zinc vapor in a magnetic field. Billeter’s method was indirect 
and did not give an accurate determination of the structure. He did not deter- 
mine the sense of the shift, but assigned it by analogy to the shifts observed in 
the 6s?1.S)—6s6p*P, transition in Hg I. Auslander (3, 4), in an experiment 
similar to Billeter’s, was unable to detect any structure in 2139A. 
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Discussion 


In the following the separations between the Zn* and Zn** components alone 
will be discussed. Everything said applies to the separations between the Zn“ 
and Zn® components except that they are half as great. Shifts are considered 
positive when the component due to the heavier isotope has the higher wave 


number. 
TABLE I 


ISOTOPE SHIFTS IN THE ZN I RESONANCE LINES 


Observed shift Normal mass shift | _ Residual shift 


2139A +0. 033 +0. 023 | +0.010 cm.- 
30764 +0. 046 +0. 016 | = +0.080 cm. 


} 


In Table I are shown the observed isotope shifts together with the normal 
mass shifts and the residual shifts found by subtracting the latter from the 
former. Since the field effect theory predicts equal shifts for the two lines with 
the sense opposite to that of the observed residual shifts, it cannot be the sole 
cause of these shifts. On the specific mass theory the difference, 0.020+0.003 
cm.~!, between the residual shifts corresponds to a single integral, C(4p, 4s) in 
the notation of Bartlett and Gibbons (6). The observed difference has the sense 
given by the theory and must therefore be regarded as a specific mass effect. 
As it is comparable in magnitude to the observed shifts, it follows that the 
specific mass effect makes an important contribution to the observed shifts and 
cannot be neglected in considering isotope shifts for elements of intermediate 
atomic numbers. 


A rough calculation of the field effect shift between Zn** and Zn® in the 
resonance lines leads to the value —0.03 cm.—! on the assumption that the 
charge is uniformly distributed in the nucleus and that the radius of a nucleus 
is proportional to the cube root of the mass. This value is probably an upper 
limit since the isotope shifts in heavier elements are smaller than that predicted 
by the field effect theory for this model of the nucleus (11). Since the field 
effect also is of the same order of magnitude as the observed shifts, but in the 
opposite sense, the observed residual shifts must be regarded as the resultants 
of both the field effect and the specific mass effect. 


The average of the residual shifts of the resonance lines is 0.020 cm. If the 
field effect is about —0.03 cm.~', the average specific mass shift must be about 
+0.05 cm.—' Since the integral C(4p, 4s) cancels out in the average of the 
residual shifts of the two resonance lines, +0.05 is an approximate measure of 
the specific mass effect that corresponds to the following sum of integrals (6), 


C(4p, 1s) + C(4p, 2s)+ C(4p, 3s) + 10C(3d, 4p) — 3C(2p, 4s) — 3C(3p, 4s). 


The experimental proof of the importance of the specific mass effect in Zn I 
has a bearing on the interpretation of the isotope shifts observed for other 
intermediate elements and the change of the sense of the isotope shift with 
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increasing atomic number. Schiiler and Westmeyver (22) have compared the 
shifts in corresponding lines of Cu I, Zn II, Cd II, and Hg II arising from the 
transitions between the levels of the configurations d‘s? and d'°p. In all these 
spectra the observed shifts have been arbitrarily assigned to the d°s? configura- 
tion by analogy with Hg II where the field effect alone will qualitatively explain 
the shifts. This assignment of shifts leads to the result that the electrons of the 
d°s* configuration are more firmly bound in the heavier isotopes of copper and 
zinc. If this result is interpreted on the basis of the field effect alone, the 
nucleus of the heavier isotope has the smaller volume. This conclusion is now 
questionable in the light of the relatively large specific mass effect established 
in Zn I. In Cd II the sense of the shift in the d*s? configuration is that pre- 
dicted by the field effect theory and is opposite to that in Cu I and Zn IT. 


It is possible to explain the sense of the shift in Cu I and Zn II. The specific 
mass shifts for both the d'p and d°s? configurations result in the energy levels 
of the heavier isotope lying higher. This is the same sense as the field effect 
shift so that here the two effects are additive. The d'°p configuration contains 
no s electron and will have negligible field effect shift. The d°s? configuration, 
with two s electrons, has a larger field effect shift. Thus if one assumes that 
the specific mass shift for the d'°p configuration is larger than the sum of the 
specific mass and the field effect shifts of the d°s* configuration the sense of the 
observed line shifts is explained. 


Since the field effect shift increases with increasing atomic number and the 


specific mass effect probably decreases owing to the factor (1/4,— 1/Ms2), the 
sense of the shift in these transitions may be expected to change with increasing 
nuclear charge. The changed sense is observed :n Cd II. However, the data 
available are not sufficient to verify this suggestion. To establish such an 
explanation the shifts in the levels relative to a datum of zero shift for ioni- 
zation must be determined uniquely by measuring the shifts in a large number 
of lines. 


The discovery of a relatively large specific mass effect in the Zn 1 inter- 
combination resonance line has a bearing on Mrozowski’s conclusion (18) that 
there is evidence of a ‘“‘pure nuclear’, or field effect isotope shift in the band 
spectrum of ZnH. The specific mass effect is not negligible as assumed by 
Mrozowski and possibly could account for the small residual shift (0.007 cm.—!) 
attributed by him to a pure nuclear effect. 


The magnitudes of field effect shifts give information about the nucleus (11), 
but, as shown by the results of this investigation, these cannot generally be 
determined for the intermediate elements without numerical calculations of the 
specific mass shifts. The nature of the integrals that must be evaluated is such 
that precise values can only be obtained when very accurate wave functions 
are used. Until accurate integrals are available the prospect of getting signi- 
ficant information about nuclei of intermediate mass from the magnitudes 
of isotope shifts is not good. 





CRAWFORD ET AL.: ISOTOPE SHIFT -IN THE RESONANCE LINES OF ZINC 


References 


. ANDREW, A. and CARTER, W. W. Phys. Rev. 74: 838. 1948. 
2. ARROE, O. H. Phys. Rev. 74: 1263. 1948. 
3. AUSLANDER, J. Helv. Phys. Acta, 11: 562. 1938. 
. AUSLANDER, J. Helv. Phys. Acta, 12: 168. 1939. 
5. BARTLETT, J. H. Nature, 128: 408. 1931. 
. BartLett, J. H. and Gippons, J. J. Phys. Rev. 44: 538. 1930. 
. BILLETER, W. Helv. Phys. Acta, 7: 524. 1934. 
3. BRADLEY, L. C. and Kuun, H. Nature, 162: 412. 1948. 
9. Breit, G. Phys. Rev. 42: 348. 1932. 
. CRAwForD, M. F., Gray, W. M., ScHawtow, A. L., and KELty, F. M. J. Optical Soc. 
Am. 39: 888. 1949. 
. CRAWFORD, M. F. and ScHawtow, A. L. Phys. Rev. 76: 1310. 1949. 
. Frep, M., Tomkins, F. S., and Bropy, J. K. Phys. Rev. 75: 1772. 1949. 
. Hess, D. C., INGHRAM, M. G., and HaypDEN, R. J. Phys. Rev. 74: 1531. 1948. 
4. HuGues, D. S. and Eckart, C. Phys. Rev. 36: 694. 1930. 
5. LELAND, W. T. and Nrer, A. O. Phys. Rev. 73: 1206. 1948. 
. LYSHEDE, J. M. and RasMussEN, E. Z. Physik. 104: 434. 1937. 
. MANNING, T. P. Phys. Rev. 76: 173. 1949. 
3. MrozowskI, S. Phys. Rev. 58: 597. 1940. 
. OpECHOWSKI, W. and DEVRIEs, D. A. Physica, 6:913. 1939. 
. Racan, G. Nature, 129: 723. 1932. 
. ROSENTHAL, J. E. and Breit, G. Phys. Rev. 41: 459. 1932. 
. ScHULER, H. and WeEsTMEYER, H. Z. Physik, 82: 685. 19338. 
. Vinti, J. P. Phys. Rev. 56: 1120. 1939. 





FORBIDDEN TRANSITIONS IN DIATOMIC MOLECULES! 
I. THE QUADRUPOLE ROTATION-VIBRATION SPECTRUM OF H: 


By G. HERZBERG 


Abstract 


Four lines of the 2-0 band and four lines of the 3-0 band of the quadrupole 
rotation—vibration spectriim of Hz have been observed and measured using effective 
absorbing paths of 10 to 50 km. atm. From this spectrum improved values of 
the rotational and vibrational constants of H¢ in its electronic ground state have 
been obtained. 


A. Introduction 


The possible transitions between the various energy levels of diatomic mole- 
cules are determined by well known selection rules.* These selection rules 
depend on the type of radiation considered (electric dipole, magnetic dipole, 
electric quadrupole radiation, etc.) and on the coupling conditions in the 
molecule. By far the strongest transitions are those corresponding to an alter- 
nating electric dipole moment of the molecule, and usually Hund’s cases (a) 
or (b) are good approximations to the coupling conditions. Transitions vio- 
lating the selection rules holding under these conditions are designated as for- 
bidden transitions. 


Three types of forbidden transitions may be distinguished: 


(1) Those that violate a selection rule that holds only in a certain approxi- 
mation or for certain coupling conditions; for example, the rule that states of 
different multiplicity do not combine with one another (AS = 0); 


(2) Those that violate a selection rule that holds rigorously for electric 
dipole radiation; they may occur weakly as magnetic dipole or electric quadrupole 
radiation; 

(3) Those that violate a rigorous selection rule but occur on account of the 
presence of external fields or at elevated pressures (enforced dipole radiation). 


Several examples of each of these three types of forbidden transitions have 
been found.* Some of these occur in absorption in the earth’s atmosphere, 
others in emission in the spectrum of the night sky; still others are suspected 
to play a role in the spectrum of the night sky and in the spectra of planets. 
A further study of these forbidden transitions may therefore be expected to 
shed more light on the problem of the upper atmosphere of the earth and of 
planetary atmospheres. In addition such studies will lead to a better under- 
standing of the operation of the selection rules, of the electronic structure of 
the particular molecules, and to a more precise determination of certain impor- 

1 Manuscript received November 22, 1949. 
Contribution from the Division of Physics, National Research Laboratories, Ottawa, 


Canada. Issued as N.R.C. No. 2086. 
* See, for example, Herzberg (9). 
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tant molecular constants. On the basis of such considerations a series of inves- 
tigations of forbidden transitions in diatomic molecules has been initiated. 


The present first paper in this series deals with the quadrupole rotation- 
vibration spectrum of H2.* A second paper will deal with the forbidden near- 
ultraviolet bands of Oz. Further work on other molecules is being prepared. 


B. Theoretical Considerations 


It is well known that the occurrence of ordinary infrared vibration spectra 
of diatomic molecules depends on the presence of a permanent dipole moment. 
For homonuclear molecules the dipole moment is zero and therefore no dipole 
vibration spectrum can occur. This conclusion may be verified in the energy 
level diagram of Fig. 1 by noting that the transitions with AJ = + 1 would 


5 -a 


-a 


2 +S 
| = 
0 +s 


Fic. 1. Energy level diagram for the quadrupole rotation-vibration spectrum of a homonuclear 
molecule. 


violate the rule that symmetric rotational levels do not combine with anti- 
symmetric ones while’transitions with AJ = 0 would violate the parity rule 
(positive rotational levels combine only with negative ones). Since no other 
AJ values are allowed for dipole radiation it follows that the dipole vibration 


*For preliminary notes on this work see Herzberg (7, 8). 
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spectrum is rigorously forbidden. However, as was pointed out several years 
ago (6), since homonuclear molecules may have a quadrupole moment, a quad- 
rupole vibration spectrum may occur. 


The selection rules for quadrupole radiation are 
4J=0,+1,+2 (1) 
with the restriction that J’+ J” must be at least 2, and the rule 
t+, - SH, +o, (2) 


where —|— stands for ‘does not combine with”. On the basis of these selection 
rules and the additional general rule that symmetric do not combine with 
antisymmetric rotational levels it is immediately seen from Fig. 1 that only 
the three branches with AJ = + 2,0, —2canoccur. According to the inter- 
national nomenclature these branches are designated as S, Q, and O branches. 
The first lines in these branches are those with J = 0, 1, and 2 respectively. 


In the case of the H. molecule James and Coolidge (10) have made more 
detailed calculations of the variation of the quadrupole moment with inter- 
nuclear distance and have predicted the relative and absolute intensities of the 
1-0, 2-0, and 3-0 infrared vibration bands. According to them the integrated 
absorption coefficients of the first lines of the S branches of these bands are 
in the ratios 8.1 X 107°, 7.3 X 107°, and 1.5 X 107° respectively to the first 
line of the R branch of the (allowed) fundamental of HCl. According to 
Herzberg (6) this means that the minimum path lengths required for the Hy. 
quadrupole bands are 2.5, 2.7, and 13.0 km. atm. respectively, assuming a 
temperature of 150°K. (at room temperature the minimum path lengths would 
be slightly higher). It may be noted that the intensity decreases in the series 
1-0, 2-0, 3-0 much more slowly than for ordinary dipole vibration spectra, for 
example that of HCI. This is due partly to the fact that the absorption co- 
efficient for quadrupole radiation is proportional to »* rather than » as for dipole 
radiation and partly to the unusually large mechanical and electrical anharmo- 
nicity of the H, molecule. 


James and Coolidge have also given theoretical expressions for the line 
strengths of the individual lines. They are (disregarding the intensity alter- 
nation due to the nuclear spin): 


(J + 1)(J + 2) (3) 
(27 + 3) 


2)(J + 1)(2J + 1) 

3(2J — 1)(2J + 3) 
9 - I=). 

(2J — 1) 


for AJ = + 2(S branch): 


ton As = 0 (Q branch): (4) 


for AJ = — 2 (O branch): (5) 


It is easily seen on the basis of these formulae that the intensity within each 
branch follows a course similar to that in ordinary infrared bands. However, 
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in the case of H2 because of the small moment of inertia the intensity maxima 
are expected at the first or second line of each branch. 


With the long optical path recently obtained by multiple traversals (2) it 
appeared possible to surpass the minimum absorbing paths required for ob- 
taining the quadrupole rotation—vibration spectrum of He. Therefore an 
attempt was made to photograph this spectrum. Some more general interest 
attaches to this problem since up to now no case of quadrupole radiation had 
been found in molecules and no case of quadrupole absorption for any atomic 
system; all previously established casés of quadrupole radiation were observed 
in emission (e.g., the auroral green line). 


C. Observed Spectrum 


The large absorption tube described in more detail in a previous publication 
(2) is 22 m. long and has a diameter of 25 cm. Using White’s method the light 
from the filament of a small lamp can be made to traverse the tube up to 240 
times, depending on the quality of the mirrors and the spectral region to be 
photographed. Thus an actual path length of 5300 m. can be obtained. By 
increasing the pressure to 10 atm. an effective path length of 53 km. atm. can 
be reached which is well above the minimum path length predicted for the 3-0 
band of the quadrupole spectrum of Hoe. 


The absorption tube was therefore filled with H», to a pressure of about 9 
atm. and the spectrum in the region of the 3-0 band (8000 to 8600 A), for about 
240 traversals, was taken in the second order of the 21 ft. grating spectrograph 
at Yerkes Observatory. Absorption lines were indeed found near four of the 
predicted wave lengths which were expected to be the strongest. However, 
since there is a fairly strong H,O band in the same spectral region and since the 
Hz» lines do not form an obvious band but are some 200 A apart from one 
another it is necessary to compare the spectrum obtained in He with one in air, 
both containing water vapor. This comparison is illustrated in Fig. 2a for three 
of the lines of the 3-0 Hz band. The reality of the H, lines is clearly demon- 
strated by this comparison. 


In the region of the 2-0 band (11300 to 12400A) the plate sensitivity is so 
low that only 100 or even fewer traversals could be used. Fortunately the 2-0 
band is stronger than the 3-0 band (see above) and there is no difficulty in 
observing it even though again a strong H2O band occurs in this region. Fig. 25 
shows three of the four observed lines obtained in the first order of the 21 ft. 
grating. 


Table | lists the wave lengths and wave numbers of the eight lines that were 
found. All but S(O) and S(2) of the 2-0 band and S(2) of the 3-0 band were 
measured on two or more plates. Second and third order iron lines were used 
as ‘standards, and the vacuum correction was calculated using the index of 
refraction data of Barrell and Sears (1). 








148 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


It is noteworthy that the Hg, lines are appreciably sharper than the accom- 
panying H,0O lines and in fact as sharp as the slit width and resolution of the 
spectrograph will allow, in spite of the comparatively high pressure of 5 to 10 
atm. No systematic investigation of the pressure dependence of the H2 absorp- 
tion was made but a qualitative comparison of exposures made at different 
pressures indicates a variation with the first power rather than the square of 
the pressure. 

TABLE I 
QUADRUPOLE ROTATION-VIBRATION LINES OF He 




















Designation air | Pvae (observed) | o-C1 O-C2 
Q(1) 12379 .90,4 8075 .39s | -+2.598 +0.014 
2-0 | S(0) 11892.49, 8406 . 36, +2.427 -—-0.028 
| S(1) 11618.975 8604.25, +2.343 4+-0.018 
\ SQ) 11379 .22, 8785.54 +2.166 --0.003 
{ Q() 
3-0 4 S(O) 8272.686 12084.65; +3.766 +0.022 
{ S(1) 8150.680 12265 .54s +3.534 —0. 
S(2) 8046. 438 12424. 44; +3.293 +0.001 
A cm! cm.—! ema 





8497.491 11764.94, +3.908 —0.012 





D. Rotational and Vibrational Constants of H; in its Ground State 


The observed wave numbers of the Hz quadrupole lines deviate slightly but 
systematically from those calculated from the data derived by Jeppesen (11), 
Urey, Teal, and MacWood (14, 13), from the ultraviolet spectrum and the 
Raman spectrum. These deviations are given in the fourth column (0-c;) of 
Table J. The reason for these deviations lies partly in the much lower wave 
number accuracy in the ultraviolet compared to the photographic infrared 
(0.1A corresponds to 6.9 cm.—! at 1200A and 0.069 cm.— at 12000A) and partly 
in the lack of reliable standards in the ultraviolet which may lead to rather 
large systematic errors. 


The vibrational and rotational energy of a diatomic molecule (in cm. units) 
is given by 


T(v, J)= G(v)+ Fi(J), (6) 
where 

G(v) = w.(v + 3) — wx.(v + $+ wy.(0 + P+... (7) 

Fy(J) = ByJ(J +1) — DLT +:1)2°+ AJIT +13 4+... (8) 

B,= B.— a.(v + 3)+ ve(v + $)?+... (9) 

Dy= D.+ Bev + $)+..- (10) 


with the usual meanings for the various symbols. 


Unfortunately the eight quadrupole lines measured are not sufficient to 
determine all the vibrational and rotational constants independently. Nor is 
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Fic. 2. Spectrograms of several lines of the quadrupole spectrum of H2. 
(a) the 3-0 band (second order of 21 ft. grating, equivalent path 45 km. atm., exposure time 
60 min.) 
(b) the 2-0 band (first order of 21 ft. grating, equivalent path 12 km. atm., exposure time 
10 to 40 hr.) 
In each case the bottom strip shows the spectrum obtained with hydrogen, the top strip that obtained 


with air, both containing some water vapor. 


i 
i 
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it possible to determine the rotational constants by the usual method of 
combination differences. An estimate of the consistency of the data can be 
obtained by evaluating the differences of corresponding lines in the two bands 
(ie. S,-3(1) — S,=2(1), etc.). These differences equal the differences of corres- 
ponding rotational levels (of equal J) in the two vibrational levels v = 3 and 
v = 2 and according to Equations (6) and (8) are given by 


T(3, J) —T(2, J) =AGoy, — (Be — Bs) (J + 1) + (D2—D3) J{J + 1)?-... (11) 


where AG2:; = G(3)— G(2) and where J refers here to the upper state of the 
line considered. In Fig. 3 these differences {after addition of 2.7/(J + 1) in 
order to take out most of the slope] are plotted against J(J + 1) as circles. 
The slope of the line thus obtained near J = 0 added to 2.7 gives B, — B; while 


Ty(J)-Tylv) 
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Fic. 3. Energy differences of corresponding rotational levels with v = 3 and v = 2 as obtained 
from the quadrupole spectrum (0) and from the ultraviolet Lyman bands (x) and Werner bands 
(*) of hydrogen. 


the deviation from a straight line is due to the quartic term. Using Urey 
and Teal’s D, values and a B, — B; value corresponding to the final B, values 
(see below) the full-line curve in Fig. 3 is obtained which represents the ob- 
served values (circles) very satisfactorily. 


The same differences T(3, J) — T(2, J) can also be obtained from the ultra- 
violet Lyman and Werner bands measured by Jeppesen (11). Using the P and 
R branches of the 1-2 and 1-3 Lyman bands and the P, Q, and R branches of 
the 0-2 and 0-3 Werner bands and plotting against the J value of the lower 
state the points represented by small crosses and dots respectively in Fig. 3 
are obtained. It is seen that these points show considerable scatter. Moreover 
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the points derived from the Werner bands deviate systematically from those 
derived from the Lyman bands by an amount of the same order as the o0-c, 
values in Table I. 


The method adopted to improve the vibrational and rotational constants of 
H; on the basis of the measurements of the quadrupole spectrum was as follows: 
The individual lines were expressed in terms of the band origins v,-9 and v3_ 
and the rotational constants B., B;, De, D; and H,»= H;= H,, for example 

Sya2(1) = ve-o+ 12B,— 144D.— 2B)+ 4D + 1720H.. 
Using this and seven similar equations for the other observed quadrupole lines 
and substituting the values of Do, D2, D3, H, from Urey and Teal (14) a least 
squares solution for v2_9, v3—o, Bo, Bo, Bz was carried out. The values obtained 
for these constants and the values of w,, w.%-, @-Ve, B., a, and y, derived from 
them are listed in the last column of Table II. The preceding column of Table 

















TABLE II 
ROTATIONAL AND VIBRATIONAL CONSTANTS OF He IN ITS GROUND STATE 
From the ultraviolet spectrum From the quadrupole spectrum 
Bo 59.331 cm— 59.318 cm—! 
Bz 53.506 ee 53.480 , 
B; 50.673 = 50.635 5 
B. 60.848 #5 60.809 + 
Qe 3.0670 - 2.993 
Ye 0.06840* “ 0.024, ms 
De 0.04648“ ) 
Be —0.00134 * ( Taken over 
Al, 5.18 X 10-5 
Vi-o 4160.27 “ J 
Vo 8084.44 ve 8087 .05o re 
V3—0 11778.34 11782 .31o ry 
We 4405 . 30 si 4395 .24 * 
WeXe | 125.325 Es 117.99; af 
Wee +1.9473¢ “ +0.29; - 
i 0.74142 K 10-8 cm. 0.74166 X 10—-* cm. 
ke 5.7600; 105 dynes/cm. | 5.73379 X 10° dynes/cm. 





*Based on a four term formula with —0.0065(v + 4)8 as the last term. 
{From the Raman spectrum [Teal and MacWood (13)]. 
Based on a four term formula with ®.ze= —O0.11265. 


II contains the corresponding constants of Urey, Teal, and MacWood (13, 14). 
As a check, using the new set of constants the positions of the individual 
quadrupole lines were calculated. The deviations from the observed values 
are listed under o-c. in Table I. Considering that only five constants were 
determined from the eight lines, the agreement of the observed and calculated 
values within + 0.02 cm.—! must be considered significant and indicates that 
the new constants represent a genuine improvement. An even better agreement 
might have been expected if it were not for the fact that all the lines except 
S(1) are rather weak and difficult to measure. 


The constants D,, 8,, and H, used here were obtained by Jeppesen (11), and 
Urey and Teal (14) from their values of B,, a,, w,, w.«. by means of theoretical 








HERZBERG: FORBIDDEN TRANSITIONS IN DIATOMIC MOLECULES. I. 151 


formulae and not by direct determination from the combination differences.* 
When the new constants of Table II are substituted in the same formulae for 
D., 8., and H, the following values are obtained 
D,= 0.04656, 8.= — 0.00180, H.= 5.016 x 10 

of which the first two deviate rather more than expected from the original set 
of constants. Using these constants for a second least squares solution for 
¥2—0, ¥3—o, Bo, Bz, Bs yields values that give a much poorer representation of 
the quadrupole lines (deviations + 0.05 instead of + 0.02) and this would be 
still worse in the next higher approximation. Therefore the D., 8., and, H, 
values of Urey and Teal (14) are considered preferable. As shown by Dunham 
(4) a breakdown of the theoretical formulae for D,, 8,, and H, may be expected 
particularly for a molecule like H»2 for which B,?/w,? is especially large. More- 
over it must be remembered that the B, and w, values given here are based on 
only three B, and AG,,; values with low v respectively which indicate a 
positive curvature of the B, and AG curves while for higher v values a negative 
curvature is certainly present. 


While it is considered that the new rotational and vibrational constants are 
superior for a representation of the first four vibrational levels, they are defin- 
itely not applicable to the higher vibrational levels. It must also be emphasized 
that the impossibility of a direct determination of the D, values introduces 
an appreciable uncertainty into the B, values. The AG values are much less 
affected by this uncertainty. 


The old and new values for the force constant k, and the internuclear dis- 
tance r, [using modern conversion factors, see Herzberg (9)] are included in 


Table II. 
E. Discussion 


Shortly after the preliminary announcement of the above results, Crawford, 
Welsh, and Locke (3, 15) reported the observation of the fundamentals (1-0 
transitions) of Ns, Os, and H2 in absorption at high pressure. This absorption 
is clearly due to induced dipole radiation since it varies with the square of the 
pressure. In agreement with this interpretation the line-width is of the order 
of 250 cm.—'. In contrast the line-width of the He lines here observed is as 
small as the resolving power of the spectrograph used permits, that is, smaller 
than 0.05 cm.~! even at 10 atm. pressure. Such a small line-width is expected 
for quadrupole radiation. This together with the fact that the intensity of 
absorption is close to the predicted value leaves no doubt that the spectrum 
here reported is a genuine quadrupole rotation—vibration spectrum. t 


Since the absorbing path used here is considerably larger than that used by 
Crawford, Welsh, and Locke at the same pressure (10 atm.) one might have 
expected to find the induced dipole spectrum in addition to the quadrupole 


*Fig. 3 shows that the scatter of the ultraviolet measurements is much too large for such a direct 


determination. 

tThe original suggestion of Crawford, Welsh, and Locke (3) that the spectrum here considered 
is an induced dipole spectrum was withdrawn in their second note (15). The similar suggestion by 
Mizushima (12) has also been withdrawn (private communication). 
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spectrum. No evidence of broad absorption bands underlying the observed 

. lines was found although this is difficult to ascertain on account of the varying 
plate sensitivity in the photographic infrared. The absence of the overtones 
of the induced dipole spectrum is, however, not surprising if it is remembered 
that some of the reasons for the slow decrease of intensity in the series of 
quadrupole bands do not apply to dipole bands. 


The first discussion of the H, quadrupole spectrum (6) was made with a view 
to detecting the large amounts of Hz assumed to be present in the atmospheres 
of the major planets. The knowledge now obtained of the correct wave lengths 
and of the intensity of absorption should help in solving the problem of Hy in 
planetary atmospheres. Unfortunately planetary spectra thus far available 
have a much lower resolution than the laboratory spectra used here and do not 
show the quadrupole lines. Since the Hg lines are very sharp the amount of He 
used here can be considered as an upper limit to the amount of H: in planetary 
atmospheres only if the same resolution is employed. 


It would be interesting to ascertain whether similar quadrupole rotation— 
vibration spectra are observable for O2, Ns, and other diatomic molecules. 
However it would be difficult to predict the intensities of these spectra. At 
any rate the decrease of intensity in the series of overtones would be expected 
to be faster and since the fundamentals would be practically unobservable 
because of overlapping by the dipole transitions of HO and COs it seems likely 
that much longer absorbing paths are required than for He. In the case of No, 
Goldberg (5) has looked for the 2-0 band in the solar spectrum taken at very 
low sun with an equivalent path of about 50 km. atm. but did not find any 


absorption. 
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A SATURATED CORE RECORDING MAGNETOMETER! 


By D.C. RosE AnD J. N. BLoom 


Abstract 


A recording magnetometer operating on the saturated core inductor principle 
has been constructed. An inverse feed-back system supplies a current to neut- 
ralize the field being measured. The neutralization is such that the sensitive 
element is maintained within 2 or 3 gammas of zero field. The feed-back current 
is measured by a recording milliameter. Using a feed-back coil of suitable 
dimensions, the record gives a direct measure of the field. 


Introduction 


In this instrument the saturated core method of measuring small magnetic 
fields has been adopted to make a recording instrument of high accuracy and 
stability. The instrument is absolute, in that the field is calculated from the 
current in a coil of known dimensions. Though the instrument was devel- 
oped to record variations in the earth’s magnetic field, it should be useful for 
any problem where fields from a few gammas to a few hundred gammas are 


to be measured. 


The saturated core method was applied as early as 1931 for measuring fields, 
and was used in a magnetic compass (3). It was developed by Vacquer, Simon, 
and Hull (5) as an airborne instrument. During the recent war the principle 
was used extensively in a magnetic detector for various applications, for ex- 
ample, the U.S. Navy Magnetic Airborne Detector. 


Principle of Operation 


The sensitive element of the instrument consists of a strip of high permea- 
bility alloy such as permalloy or u-metal. The strip used in this instrument is 
0.006 in. thick, about 23 in. long, and 1/8 in. wide. A wide variation in these 
dimensions makes little difference in the sensitivity provided the coils are 
suitably designed. This strip of permalloy is excited by an alternating field of 
sufficient strength to be well over saturation on the peak of each half cycle. 
The excitation is supplied by a coil about equal in length to the strip. Fig. 1 
shows a sketch of the sensitive element of the instrument. The inner coil is 
the exciting coil, the outer being for inverse feedback. 


If the permalloy strip is excited by an alternating voltage of pure wave form, 
the resulting current wave will be distorted by the shape of the B-H curve. If 
there is no field other than the exciting field, the distorted wave will be sym- 
metrical on the + and — sides of the zero axis, so will contain only odd har- 
monics. If the instrument is in an external field (the earth’s field or the field 
to be measured) the exciting voltage will saturate the core more easily in one 


1 Manuscript received November 9, 1949. 
Contribution from the Division of Physics, National Research Laboratories, Ottawa, Canada. 
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direction than in the other. The wave form then will not be symmetrical about 

the zero axis and will contain even harmonics the strength of which is a function 
soe 

AO i 


of the external field to be measured. 
YA TILL. ! i | \ 


" EXGITING WINDING FEED BACK WINDING 


MULL pect! y 


Fic. 1. Schematic diagram of sensitive element. 


Fic. 2. The B-H curve, idealized with no hysteresis. 


This can be shown analytically. If one considers the ideal case where the 
hysteresis loop has no area, and the curve has corresponding equal values in 
the first and third quadrants, the magnetization curve (Fig. 2) can be repre- 
sented by an equation in the form of a power series 


B’ = a,H’+ a;H"+ asH” + 


Applying a small external field (so small that it is on the linear part of the B/H 
curve) is equivalent to shifting the axis an amount 


B’= B+ uH, 
H'= H + H., where H, is the small external field. 


H is applied as a sinusoidal magnetizing force Ho sin wt, so the equation 
becomes B+ uH,= a,(Ho sin wt + H,) 
+ a3(Ho sin wt + H,)' 
+ a;(Ho sin wt + H,)'+..., ete. 
This equation may be reduced to the form 
B = —uH,— Pot+ Pisin wt + P2 cos 2wt + Ps; sin 3wt + 
where the P’s are power series in Hp and H,. The P’s with even subscripts 


have no terms which do not depend on H,. The problem then is to select an 
even harmonic or group of harmonics, and a measure of the amplitude of these 
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gives a measure of H,. Experience shows that when the second harmonic is 
chosen the intensity of this harmonic is roughly proportional to H, over a 


considerable range. 


In our instrument the field is measured by amplifying the second harmonic, 
rectifying it to a direct current, and feeding it back through the outer coil 
(Fig. 1). The outer coil is of known dimensions and is sufficiently long to give 
a reasonably uniform field over the length of the permalloy strip. The feed 
back neutralizes the external field H, and, if the feed-back constant is suffi- 
ciently high to leave the residual field negligibly small, the field is measured 
by the current and dimensions of the feed-back coil. There are two well known 
ways of selecting appropriate even harmonics. In the first only filters are used. 
In the second two strips of u-metal are used instead of one. These are placed 
side by side, each with an exciting coil, connected in opposition. A pickup coil 
on each connected in series, or one coil surrounding both, will pick up the 
harmonics not balanced out by the opposition connection of the primaries. 
Both methods are practical, but in the instrument described here the method 
wherein the second harmonic is separated by high quality filters was chosen. 
This has the advantage that the sensitive head is geometrically simple, and 
only one winding (apart from the feed-back coil) is required. It has the dis- 
advantage that special filters have to be built.* The U.S. Navy M.A.D. 
equipment (4) uses this single core method, and the design used here is a 
modification of theirs. 


The Oscillator, Filters, Amplifier, and Phase-Sensitive Detector 


The exciting frequency chosen was approximately 1000 c.p.s. (Actually the 
filters as constructed are peaked about 990 cycles and the oscillator is tuned 
to this.) As the same coil is used for excitation and for pickup the exciting 
voltage must have no second harmonic and must be very stable. Fig. 3 shows 
the circuit diagram of the oscillator and filters. The oscillator is a well known 
type wherein the feed-back circuit contains a resistance capacity network, which 
fixes the frequency. After amplification the exciting voltage is fed to the coil 
through the 1000 cycle filter. The second harmonic is selected by connecting an 
appropriate filter across the exciting coil in parallel with the 1000 cycle filter. The 
values of the filter components are as indicated. These filters were manufac- 
tured to our specifications by the United Transformer Company. The 2000 
cycle filter has an impedance transformation with an actual voltage step up 
of about six. Both filters were designed so that there would be a correct im- 
pedance match and the exciting coil was wound to give the correct impedance 
(1000 ohms at 1000 cycles) with the permalloy strip in place. Care was also 
taken_to match the loads on the other sides of the filters. The 1000 cycle filter 
attenuates 2000 cycles by about 80 db. The oscillator produces only a fraction 


*Another successful recording instrument using the two core method of separating the required 
harmonics has been developed by P. Serson of the Dominion Observatory, Ottawa. We wish to thank 
Mr. Serson for showing us the details and results of the use of his instrument. Double-core instru- 
ments had also been made in the N.R.C. laboratory by Ralph Bailey (1) for geomagnetic prospecting. 
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of 1% of second harmonic. The sensitivity of the system varies with different 
coils and y-metal strips. An average value at the output of the 2000 cycle filter 
is about 25 uv. per gamma of H,. With the circuit as used any spurious second 
harmonic that might get through from the oscillator would be about 1077 v. so 
the sensitivity is not limited by unwanted 2000 cycle from the exciting voltage. 


TO FREQ 


© TO AMPLIFIER 


Fic. 3. Oscillator and filters. 


C\— 63? pf. 
Ca, C;, C:—0.1 uf. 


R,—0.25 meg. 
R.—1000 ohm, variable 


R;—0.25 meg. 

R;—56K 

R;—0.47 meg. 

R;—389K 

R;—0.4 meg., variable 
Rs, Ro—1.5K 

Rio 2 5.6K 

1000 cycle filter 

Co, Cio—1164.1 wu. 
Cn, Cre, Cis —14788 Mu. 


C;—100 wyf., variable 
C;—620 wf. 
Cy, Co— 10 ef. 
Coo —20 pf. 
T,\—Primary, 10,000 ohms 
— Secondary, (a) 20,000 ohm CT 
(b) 1000 ohm 


2000 cycle filter 
Cu—1994 pf. 
Cis—126. 2 pups. 


Li—1.591 h. 
ly -1.498 h. 
L3—1.591 h. 


Cie—98474 pp. 
Cy IIS Myf. 
Cis—871 up. 

CG g—l {81 Myf. 
Ly—1.591 h. 
Ls—6.366 h. 
Le—0.0638 h. 
L-—8.988 h. 


The amplifier circuit is shown in Fig. 4. The 2000 cycle filter attenuates 
1000 cycles by 80 db. relative to the 2000, and an additional 30 db. rejection 
of the 1000 cycle signal is achieved by a tuned feed-back network in the second 
stage of the amplifier. This network, described by Valley (6, pp. 398-408), has 
a high inverse feedback except at the desired frequency. In spite of the atten- 
uation a certain amount of 1000 cycle gets through. This, with the usual noise, 
fixes the limit of sensitivity of the instrument. 


There is an optimum value of the exciting voltage, a typical curve being 
shown in Fig. 5. The voltage used is about 3.8. The 2000 cycle filter and 
amplifier reduce this, relative to the 2000 cycle signal, by 80 + 30 = 110 db. 
or a factor of about 3.2 & 10° in voltage. There will be, therefore, about 12 pv. 
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of 1000 cycle unwanted, to be compared with the wanted signal of 25 uv. per 
gamma. This 1000 cycle residual can be seen when the sensitive head is turned 
to null field, the output of the amplifier being displayed on a cathode ray 


oscilloscope. 


The circuit of the phase sensitive detector is shown in Fig. 4. Sensitivity in 
direction of the d-c. output to the phase of the 2000 cycle signal is necessary 


FF 





Fic. 4. 


R,:—22K 
R,—150K 
R;—5.6K 
Ry—0.9 meg. 
R;—2.2 meg. 
R,—120K 

R;, Rs—0.12 meg. 
Rg 56K 
Rio—2.2 meg. 
Ry, —0.5 meg. 
Ri—2.2K 
Riz—0.12 meg. 
Ria —56K 
Ri;—d.6K. 


Amplifier, detector, and filter. 


Rie, Riz, Ris —lK 

Ci —4 bf. 

Cy—0.01 pf. 

C;, Cy —886 Mu. 

C;—1372 pf. 

Co—4 mf. 

C;—0.25 pf. 

Cs, Co—0.01 pf. 

Cio, Cu—4 uf. 

Cp, Cis —10 uf. 

Ly —10 h. 

T,:\—Primary, 10,000 ohm 
—Secondary, 45000 ohm each 


Bc 
1 


OUTPU 


2 


to make the sense of the feed-back current vary as the field passes through null. 
The circuit for this is a modification of that used by Cosens (2). A reference 
voltage of the same frequency and of constant relative phase to the signal 
voltage is required. This is obtained from the 1000 cycle oscillator by a simple 


nm 


RELATIVE SENSITIVITY 


Fic. 5. 


- os  & 5 6 
EXCITING FIELD 


Relative sensitivity of sensitive element at different exciting voltages. 
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circuit (not shown). It is a conventional frequency doubler using a push-pull 
rectifier to give 2000-cycle pulses. These are amplified, the wave form being 
purified considerably by a parallel-tuned circuit of high Q. A very pure wave 
form is not essential here. For maximum sensitivity the reference voltage must 
be in phase (or 180° out of phase, depending on the direction of H,) with the 
signal voltage. This is accomplished by a fine adjustment of the oscillator 
frequency to take advantage of phase shifts in the filters and amplifiers. 


The over-all gain of the amplifier from the 2000 cycle filter output to the 
output of the phase sensitive detector is of the order of 60 db. A simple low 
pass filter (critical frequency about 35 c.p.s.) is included after the detector to 
give a smooth d.c. If this is connected directly to a meter such as an Esterline 
Angus Recorder of 1 ma. full scale, a sensitivity of the order of 100 gammas for 
the full scale of the meter is ordinarily achieved. The instrument thus far, if 
calibrated frequently, would make a complete measuring instrument in itself, 
but it would not be stable in sensitivity over long periods. For convenience 
in design the oscillator, filters, amplifier, phase sensitive detector, low pass 
output filter (and power supply for these), are included in one chassis. The 
feed-back amplifier and servo system are mounted separately with their power 
supply in another chassis. This design was followed so that the instrument 
could be used independently without feedback, and also because it was found 
advisable to use a separate power supply for the feed-back section of the 
instrument. 


The Feed-back System 


In principle, the rectifier output is fed back through a coil surrounding the 
sensitive element (the outer winding in Fig. 1). It is fed back in such a direction 
as to neutralize the field almost completely, leaving the sensitive element in a 
very weak field, but still sufficient to supply the feed-back signal. If the resi- 
dual field is small enough to be neglected the current in the feed-back solenoid 
gives a measure of the field. The design of such a system has definite limit- 
ations, owing to the requirements of available recording meters and inherent 
instabilities of inverse feed-back systems. 


If the filtered output of the phase sensitive detector were fed back directly 
it would make a complete instrument. However, with the circuit as shown, the 
output at this stage would make the feed-back ratio too small to meet the 
objective of making the residual error signal negligible. 


A balanced d-c. amplifier was considered but discarded on the grounds that 
it adds an objectionable source of zero drift. The only balanced component 
at present which might cause zero drift is the phase sensitive rectifier. We 
found no way of eliminating all circuits which might cause drift, but, with the 
circuit used, this has been reduced to a negligible value. 

To make recording of the feed-back current readily possible the current 
available to cover the selected range in field intensity should be 1 ma. This 
will operate an Esterline Angus Recorder. Leeds and Northrup and Brown 
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make more sensitive recorders and these have been used with the instrument 
but it was considered desirable to design it for use with either type. One milli- 
ampere was therefore chosen as the’ feed-back current to neutralize a field of 
1000 gammas. Toreduce toa minimum, errors due to the residual field, the feed- 
back ratio N should be as high as possible. At least 100 was considered a good 
objective. WN is defined by the equation 


1 
H=He(1+4), 


where H is the field to be measured and HH, is the neutralizing field produced 
by the current in the feed back coil. 

An attempt was made to amplify the 2000 cycle signal by about another 
10 db. and connect the output of the phase sensitive detector through an 
appropriate filter directly to the feed-back coil. For this plan the filter design 
is a difficult problem. It must block the 2000 cycle current completely; it must 
have a reasonably low impedance to pass the 1 ma. Yet the whole system must 
be stable. A consideration of the stability requirement and some experimenting 
led to the conclusion that a high value-of N could not be obtained in this way 
without using capacities that were much too high to be practical, and specially 
constructed inductances. One instrument was built and operated successfully 
for several months but the best value of N achieved under practical conditions 
was about 30. 

The Servo System 


A servo feed-back system as shown in Fig. 6 was found successful. The prin- 
ciple of operation is well known. The d-c. voltage from the filtered output of 


75 VOLTS 
AC 
TO FEEDBACK 
COIL 


FROM AMPLIFIER 


6.3v ac— 0000) 


Fic. 6. Servo amplifier and feed back circutt. 


Ri -A6K C, Cy —0.25 pf. 

R.—22K C;—10 pf. 

R;—150K Cy—4 uf. 

Rs—0.5 meg. C;—0.01 pf. 

R;—0.38K Cc—10 pf. 

R;—22K C;—4 bf. 

R;, Rs—100 ohms Cs—0.01 pf. 

M,—servo motor T,— UTC A-11 

M2—dam ping generator T2—Primary, 3500 ohms 

M;—recording milliameter Secondary, 900 ohms 
B—1.5 v. dry cells. 
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the sensitive detector is broken up into pulses by a vibrating reed chopper or 
converter (a Brown converter is used). The pulses are at 60 cycles and after 
amplification drive a two-phase servo motor. The servo motor drives a poten- 
tiometer which controls current to the feed-back coil from a battery circuit as 
shown. A Diehl Manufacturing Company 60 cycle a-c. control motor was used, 
geared to a General Radio potentiometer type 214A with a gear ratio of 25 
to 1. The potentiometer has a resistance of 500 ohms and the winding is fine 
enough to give steps of less than 2 gammas. Potentiometers of this type have 
been operated for several months with no trouble. They have the advantage 
that they are cheap and easy to replace. A precision potentiometer is not 
required since the feed-back current is measured independently. The ampli- 
fication, when all available gain is used, is such that the servo motor will start 
in either direction with a signal of considerably less than 1 gamma. 


The servo motor used with the gear train and potentiometer has inertia and 
torque characteristics which make it unstable unless there is auxiliary damping. 
Damping is achieved easily by another small two-phase servo motor (mounted 
on the same shaft as the driving motor) connected to act as an induction 
generator. It is excited from the 60 cycle power supply and the output fed 
back into a winding of the converter transformer. The damping can be con- 
trolled by the exciting voltage. The damping generator is a Pioneer Bendix 
type CK-2 designed for 400 cycles but works very well at reduced voltages at 
60 cycles. The design of these servo parts is not critical; several mock up 
systems were made using assorted types of motors obtained from war surplus. 
The over-all sensitivity of the instrument is such that the servo mechanism 
will balance it so that the residual field at the sensitive element is about 2 
gammas. When using this servo system the definition of the feed-back ratio 
in the previous section does not apply, but the residual error is always small. 
Used with an Esterline Angus recording meter set for 1000 gammas full scale 
the instrument is more accurate than the meter. With a Brown or Leeds and 
Northrup recorder, precision up to the limit of the magnetometer can be 
obtained. 


Stability 


With the high inverse feed-back ratio used, the gain of the amplifier can 
change considerably without appreciable change in the accuracy of recording. 
Apart from component failures the only changes which can affect the output 
seriously are a change in frequency of the oscillator and a change in the balance 
of the phase sensitive detector. The oscillator circuit was chosen because of its 
stability in frequency and purity of wave form. Its power supply is regulated 
by VR tubes, and experience shows that its frequency is very stable. 


Five complete instruments have been made and operated continuously for 
varying periods up to several months. The instruments have been operated 
successfully without adjustment for periods up to two months but for con- 
tinuous use a regular tune-up as described below is advisable every two or three 
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weeks. Once familiarity with the instrument is achieved the adjustments take 
only a few minutes. 
Adjustments 


In tuning up the instrument the feed-back control potentiometer (R2 in 
Fig. 3) on the oscillator section is reduced until the circuit will just oscillate 
steadily. The gain control (R; in Fig. 3) feeding the exciting coil amplifier is 
then adjusted until the voltage across the exciting coil is 3.8 v. With a signal 
corresponding to about 50 gammas and the servo amplifier disconnected, the 
frequency of the oscillator should be adjusted until there is a maximum signal 
in the 2000 cycle amplifier. The frequency is adjusted by the small variable 
condenser (C; in Fig. 3). This tunes the circuit to the 2000 cycle filter system. 
The 1000 cycle filter is relatively broader. Therefore the frequency is correct 
for both if correct for the 2000 cycle filter. Phasing of the reference voltage 
and the signal to the phase sensitive detector is the next adjustment. The 
horizontal and vertical plates of a cathode ray oscilloscope are connected to 
the two inputs of the detector. Test point terminals are brought out for this. 
The gains should be adjusted to give a good pattern, then the frequency of the 
main oscillator readjusted slightly until the pattern is as nearly as possible a 
straight line. This adjustment may put the oscillator slightly off the peak fre- 
quency of the filters, but the gain is sufficiently high that this is not important. 
The next adjustment is to balance the phase sensitive detector. Connect the 
servo amplifier and turn the gain control (Ri in Fig. 4) ef the 2000 cycle 
amplifier to zero. Adjust the balancing potentiometer (Ris in Fig. 4) until the 
servo motor remains still in any position. All these adjustments except phasing 
the detector can be made with nothing more than a high resistance a-c. volt- 
meter. An indicating voltmeter built into the instrument is connected across 
the filtered output of the phase sensitive detector. If the instrument is used 
without the feed-back system this meter indicates the field measured. If the 
servo feed-back system is in operation this meter should always read so nearly 
zero that the current is not observable as different from zero. With practice, 
phasing may be accomplished by using this meter without a cathode ray oscillo- 
scope, but a 3 in. oscilloscope is very useful in making the adjustments quickly. 


When the above adjustments are complete, the gain of the 2000 cycle 
amplifier and the servo amplifier should both be turned up to give a good 
torque on the motor without instability. Too high a gain will cause the system 
to oscillate. The perfection of adjustment can be tested at any time by feeling 
the torque on the servo motor gears when they are given a small displacement 
from the balanced position. 


Design of Sensitive Element 


The handling and mounting of the permalloy strip is very important in 
making successful magnetometers. This is particularly so if wide changes in 
sensitivity are to be avoided. The metal used was Allegheny yu-metal supplied 
in sheets 6 mils thick. The strips were cut and heat treated in the laboratory. 
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The heat treatment must be carried out accurately according to instructions 
of the suppliers. After heat treatment the strips must be handled very care- 
fully so that no bending or strains are introduced. Otherwise the desired mag- 
netic properties are lost. 


The strip is mounted in the center of a lucite bobbin on which the exciting 
coil is to be wound. The bobbin is made accurately in two halves (see Fig. 1), 
being split longitudinally through its axis. Two slots are milled in the two 
halves slightly wider than the strip and making a total thickness of 
slot of 8 mils, 2 mils greater than the strip thickness. To assemble, the 
strip is fixed in the slot in one half of the bobbin with a tiny spot of soft wax 
at its center. The two halves are then pressed together and the coil wound. 





Fic. 7. Complete magnetometer showing oscillator and amplifier box, servo box and recorder, 
with sensitive element in foreground. 


The wax holds the strip at its center, the ends being free to expand and contract 
in the slot, but not to move laterally, by more than the 2 mils clearance. The 
feed-back coil is wound on a lucite cylinder about 9 in. long. The cylinder has 
an axial hole just large enough to hoid the bobbin which is placed at the center 
of the coil and stuck in with a little coil dope or wax. 


This design gives very satisfactory results. Tests were made for temperature 
variation of sensitivity, the head being placed in a box where the temperature 
was changed between 0°F. and 60°F. When the detector is connected directly 
to the meter without feed-back, a change in sensitivity of 0.3% per Fahrenheit 
degree was noted but with the feed-back connected no temperature effect could 
be observed. 


Fig. 7 shows a photograph of a complete instrument. The coil here is de- 
signed for recording variations in declination of the earth’s magnetic field. 
To do this it is set in the usual position, horizontally at right angles to the 
steady horizontal component. The recorded variations will then be H sin D. 
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where D is variation in declination angle. To use the instrument to record 
variations in horizontal or vertical force some additional system to neutralize 
the steady component must be included. Fig. 8 is a sample record of a mag- 


netic storm. 
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Fic. 8. Records of a magnetic disturbance taken with Leeds and Northrup Speedomax (top) and 
Esterline Angus milliameter recorders. Time marks are G.M.T. (The Leeds and Northrup 
recorder paper drive had been altered. Therefore the time marks do not fit the paper scale.) Full 


scale in each case represents 1000 gammas. 
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FORMATION OF RADIOACTIVE SURFACE FILMS ON MINERALS! 
By T. G. CuurcH 


Abstract 


With the requisite chemical treatment, specific minerals will assume a surface 
radioactivity by undergoing reactions with radioactive ions in solution. It is 
essential that the mineral have a suitably tarnished or coated surface, for in 
general no exchange will occur ona fresh surface. It is suggested that the radio- 
activity so induced could be used in conjunction with a mechanical device to 
separate and concentrate the valuable minerals in an ore, or to control hydro- 
metallurgical plant processes. 


Introduction 


This paper describes experiments performed to determine whether mineral 
surfaces can be made radioactive by undergoing exchange reactions with radio- 
active ions in solution. Such radioactive coatings, if selective, would serve as 
a means of differentiating between minerals in an ore, and enable them to be 
separated by a suitable mechanical ore sorting device actuated by a Geiger 
counter. Thus the metallurgist would have at his disposal an additional means 


of beneficiating certain ores. 


F. Paneth (1) has pointed out that natural minerals offer great resistivity 
to exchange reactions as compared with laboratory-prepared compounds of the 


same chemical formula. His work clearly indicated that (except in rare cases) 
no exchange could be expected on the fresh surface of a primary mineral. There 
remained, however, the possibility that tarnished or weathered mineral surfaces 
might react more readily. Therefore, in the present work attempts were made 
to coat the mineral with an adherent film of a substance which would exchange 
readily with active ions in dilute acid solutions. The experiments were of a 
preliminary nature only, to test the feasibility of the method. It is anticipated 
that every ore will present a somewhat different problem. 


Experimental 

General 
Synthetic ores were used; the constituent minerals were clean, well crystal- 
lized pieces that offered a surface for counting of about 15 mm.? Fragments 
of quartz and feldspar representative of gangue minerals were carried through 
the same treatments given the ore mineral. A few pieces of each mineral were 
placed in a beaker containing a conditioning agent designed to produce a surface 
coating. After the conditioning period the pieces were washed and transferred 
to the active solution for the required time. Then, after washing and drying, 
the activity of each piece was measured with a 110 volt a.c. operated Geiger 
1 Manuscript received in original form August 15, 1949, and, as revised, December 1, 1949. 


Contribution from the Research Division, Atomic Energy Project, National Research 
Council, Chalk River, Ont. Issued as N.R.C. No. 2079. 
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Probe Monitor fitted with an end-on counter having a window of 2.7 mgm. per 
cm.” thickness. The samples were placed 6 mm. below the window. Activities 
were recorded as counts per min. per cm.? of mineral surface. 


Galena (PbS) 

Galena oxidizes fairly readily in moist air to give a surface coating of lead 
sulphate. For this reason S* in the form of SO,” in acid solution was used in 
attempts to activate the surface by exchange reaction. To produce the desired 
coating the freshly broken mineral was given a pretreatment with dilute nitric 
acid before transferring to the active solution. The lead sulphate coating was 
found to be firm and adherent. 


No exchange was observed to occur on the fresh surface of the galena. As 
oxidation of the surface progressed, exchange took place more and more readily. 
Even when only a slight tarnishing was apparent, the surface assumed a quite 
high activity after a few minutes in the active solution. Representative tests 


are recorded in Table I. 


TABLE I 


EXCHANGE OF S*O,g IN SOLUTION WITH GALENA SURFACE 












Time i Activity, 

Conditioning | Appearance 2S ee | counts/min./em. 

treatment of coating acu solution . Seer 

: has (1 we./ce. in 0.01% H2SO,)|} A 
| Galena Gangue 

1% HNO; Bright, no tarnish 3 min. to 20 hr. j Nil Nil 

quick rinse 

1% HNO; | Slight tarnish 10 min. 4,200 | 200 

for 1 hour | | 

1% HNO; Visible coating of 10 min. | 9,800 | 350 

for 20 hr. | PbSO, 

1% HNO; | Heavy coating of 20 hr. 31,500 | 2,100 

for 20 hr. | PbSO, 


In these tests the active solution contained added (inactive) sulphuric acid, 
thus lowering the specific activity of the S*. When only the tracer SO,~ was 
present a much higher activity on the galena surface was obtained, but in this 
case some S*O,> was absorbed on the quartz and feldspar, so that the ratio of 
the activities, galena/gangue, was very little higher than when inactive SO,> 
was present. Representative tests with carrier-free SO, are shown in Table II. 


Pyrite (FeS) 

Some tests were carried out on pyrite under the same conditions as for galena. 
The oxidation of the surface in this case produces a coating of ferric oxide 
instead of sulphate, as any sulphate formed goes into solution. Consequently, 
the mineral surface could not be made radioactive by using S*O,~. 


A patite (Ca F..3Ca 3P'5Os) 
It was found that no exchange occurs between the surface of the apatite and 
P® as phosphate in acid solution. Attempts to produce a coating on the mineral 
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TABLE II 
=XCHANGE OF S*Og IN SOLUTION WITH GALENA SURFACE 





| ers 
peli Activity, 
\ ppearance Time in, counts/min./cm.? 
a eoatine active solution —| oe 
as (2 we./cc. in 0.1% HNOs)) _ | 
Galena | Gangue 


Conditioning 
treatment 


1% HNO; Slight tarnish 10 min. 
for 1 hour 


3,150 | SiO.-2800 
| feldspar- 

4200 
1% HNO; 3,000 
for 20 hr. 
1% HNO; 
for 20 hr. 


Visible coating of 10 min. 
PbSO, 
Heavy coating of 
PbSO, 


15,400 | 


20 hr. 140,000 | 3,000 


—— ee ead sali ancl enim aiadionsivoee - eh dita 


that would exchange with P®O,> in acid failed. This was to be expected since 
calcium phosphate is soluble in dilute acid. It was found, however, that dilute 
sulphuric acid produces a relatively insoluble coating of calcium sulphate on 
the apatite, and that this coating exchanges very readily with S®O,> in solu- 
tion. In this way the mineral surface could be given a high activity. Some tests 
under rather extreme conditions are shown in Table III. 


TABLE III 
EXCHANGE OF S*Oyg IN SOLUTION WITH APATITE SURFACE 


Time in Activity, 


Conditioning 
treatment 


Warm conc. HeSO, 
for 2 hr. several 


Appearance 
of coating 


Very thick coating 
of CaSO, 


active solution 


(5 we./ec. in 0.02% H2SO,4) 


15 min. 


counts/min./cm.? 


Apatite | Gangue 


2,100 


water washes. 


a | 126,000 


10% H.SO, for 20 
hr., then several 
water washes. 


260,000 


Cobalt Minerals 


Some work has been done on cobalt minerals. Tests were carried out on 
smaltite (CoAss), pentlandite (Fe, NiS) and niccolite (NiAs), using Co®’** 
as the active ion. Smaltite generally contains some nickel diarsenide and 
niccolite has many impurities, including cobalt. Owing to their different chem- 
ical compositions and physical characteristics, it is difficult to coat the three 
minerals uniformly. However, digestion with boiling 3N sodium hydroxide 
solution produces a surface film of hydrated oxides. Then, after standing in 
slightly acid active solution with a strength of 1 uc. per cc. for some hours, the 
tarnished minerals become radioactive. Some of this induced activity seems 
to be due to surface absorption and some to ion exchange. Probably close pH 
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control and the use of suitable amounts of carrier would reduce straight absorp- 
tion to a minimum. But, regardless of the mechanism of these effects, the 
activity measured is a definite characteristic of each mineral and the differ- 
ences in the activities of each are quite marked. 


Conclusion 


We may say that certain mineral surfaces can be made radioactive by 
exchange with active ions in solution, provided that the surface is coated with 
a slightly soluble compound produced by weathering or by chemical treatment. 
The number of radioisotopes available should make the method generally 
applicable. The experimental data which has been presented is rather crude 
and scanty but it is hoped that sufficient work has been done to indicate the 
general method to those who may be interested. 
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ELECTRON DIFFUSION IN A SPHERICAL CAVITY! 


By A. D. MACDONALD AND SANBORN C. BROWN 


Abstract 


The diffusion equation for electrons in a nonuniform field is solved and the 
breakdown condition derived. The breakdown condition is expressed in such 
a manner that an effective characteristic diffusion length Ae is determined; the 
meaning of A, expresses the equivalent characteristic diffusion length for unifonm 
electric fields. From the experimental breakdown fields, in which the electric 
field is uniform, A, is determined and used to predict theoretical breakdown curves 
for the nonuniform cases. Theory and experiment are compared, the agreement 
verifying the correctness of the approach. 


Electrical breakdown of a gas at microwave frequencies has been considered 
theoretically and experimentally for the case in which the electric field is 
uniform (2, 5, 6). It is sometimes important to consider cases in which the 
electric field is not uniform in order that the shape of the volume from which 
diffusion takes place may be varied. The case of a cylindrical cavity of arbi- 
trary length has been solved by Herlin and Brown (3), and the procedure of 
this paper is similar to theirs. The diffusion equation for electrons in a non- 
uniform field is solved and the breakdown condition is expressed in such a 
manner that an effective characteristic length A, is determined; the meaning 
of A, expresses the equivalent characteristic diffusion length for uniform elec- 
tric fields. From the experimental breakdown fields measured using uniform 
fields A, is determined and used in the theories developed in References (5) 
and (6) to predict theoretical breakdown curves for the nonuniform case— 
these curves are compared with experiment. 


I. Diffusion Theory 


The differential equation to be solved for the characteristic diffusion length 
of electrons in a gas is 


A*y + SE*y = 0, (1) 


where y = Dn, the electron diffusion current density potential; D is the electron 
diffusion coefficient; 1 is the electron concentration; E is the r.m.s. value of 
the electric field, and ¢ is defined by the equation 

¢ = »/DE’, 
v being the ionization rate per electron. 


The boundary condition on Equation (1) is that the electron concentration 
and thus y go to zero within a mean free path of the walls of the metal cavity. 


1 Manuscript received November 18, 1949. 
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Command, and O.N.R. 
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¢ may be expressed as a function of E/p and pid, where p is the pressure 
and ) is the free space wave length of the applied electric field. Integration of 
Equation (1) is simplified by the use of the approximation employed by Herlin 


and Brown (3, 4). E\e-? e\? (E\8-? 
=7,(£ od ae 2 
oridia: (2) (£) (=) @) 


where { is the value of the ionization coefficient at the maximum field point; 
k is introduced for mathematical convenience and has the units of reciprocal 
length. The quantity 8-2 is obtained as the slope of the ¢ vs. E/p plot on a 
logarithmic scale. The slope required is that of ¢ for which pd is kept constant. 
The approximation of Equation (2) is very good where the ionization is high. 
It is inaccurate only where there is little ionization and, since the regions in 
which there is little ionization contribute only slightly to the determination of 
breakdown fields, the procedure leads to correct results. 


The electric field in the lowest electric mode in a spherical cavity may be 
given by (10, p. 562). 


- 


-éor 


E,= E,(—* ) cos 6 j:(2.75r/a), (3) 


E,= Ea( 2) sin 6 : [ (2.75 r/a)ji(2.75 r/a)] , (4) 


7dr dr 
E,= 90, 


where 7, 6 and ¢ are the spherical co-ordinates, a is the radius of the sphere, 
and j; is the first-order spherical Bessel function. 


It is seen that the electric field depends on both r and @, the introduction of 
which makes Equation (1) inseparable. However, in breakdown, we are inter- 
ested only in the energy transfer from the field to the electrons so that we need 
take into account only the magnitude of the field at a given point. Near the 
center of the cavity, where the field, and therefore the ionization, is high, the 
magnitude of the electric field is approximately spherically symmetric. This 
fact is illustrated in Fig. 1, which shows the maximum variation of fields with 
6, as a function of r. If we assume that the electric field may be expressed as 
the average of these values over the whole of the cavity, we may write 

E = E,[(1 —(r/a)?]. (5) 
Equation (5) is also plotted in Fig. 1, where it is seen to be a good approxi- 
mation to the average electric field, except near the boundaries where it does 
not matter. 


The assumption that the electric field is independent of @ and ¢ leads one 
to the independence of y on these variables. Therefore, we may write Equation 
(1) with the assumption of Equation (2). 


1 d (nt) (EY 
— —[rr—) +h ( — = 0, 
r? dr : dr Eo ¥ 
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and introducing the value of E from Equation (5), we have 


14 (nd) + [1 — (r/a)P y =0. (6) 


r? dr dr 


We expand the term in r/a by the binomial theorem and drop powers of (r/a) 
greater than 2. This makes an appreciable error only near the boundaries, 
where again the accuracy of the method is unimportant. Then 


ad Py 2 dy 


r, 


ss +e *r’)y =0, (7) 
where p2= B/a®. 1/y is the radius at which the ionization goes to zero under 
the il assumptions. Beyond 1/y these assumptions lead to a negative ¢ 
which is not physically correct, so we set ¢ = 0 for r > 1/u. 


1.0 ————E,(8=0) 
—-——Eg (6=7/2) 
seo — AVERAGE 
ANALYTIC APPROXIMATION 


08 


E/E, 


06 


04 


0,2 





S 0 02 04 O6 08 1.0 


v/a 


Fic. 1. Electric field in a spherical cavity showing the maximum variation and the analytic 
approximation to the average field as a function of radius. 


For mathematical convenience, we transform to a dimensionless independent 
variable and let kur? = x; the equation (7) becomes 


ay, : 3 dy “( 1 =.) k 
— +-—-— +—(— —-—Jy=0 x<—. (8) 
dx? 2xdx 4ku =A? ¥ be 
We transform the dependent variable by letting 
ya=e"*s, 
and then 
d’g +28( 3 1) ag 
= a er 9 
dx? - 2x 0 (9) 
where 
3 k 
a=>—-—-—. 
4 du 
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Equation (9) is the equation for the confluent hypergeometric function in 
parameters 3/2 and a (Tables of this function are available (7, 8).) The second 
solution is not allowed by the boundary condition and therefore 


y= rua = x), (10) 


where we designate by ¥ that part of y for which 7 is less than 1/y or x < R/p. 
When x > k/u, ¢ is zero and the differential equation (1) becomes 


4 (gat) 9 
dx dx 


n= cfi-(2)"] on 
x 


where x)= kya? and is determined by the condition that ~: be zero on the 
boundary; C is an arbitrary constant, We must match the solutions of Equa- 
tions (10) and (11) at the point where ry = 1/y and therefore 


of which-the solution is 


vii _ V2 
v1 fe 
which gives us 
2/3aM(a +1; 5/2; 31) _ 


M(a; 3/2; x1) 


where 


Equation (12) relates a, the radius of the cavity, 6 determined from the slope 
of the ¢ curve and k which is inversely proportional to the characteristic dif- 
fusion length. The equation may be written in the form 


| 2B ade + 5/2; y) 7 x (13) 
M(a; 3/2; 9) x 

where a = = z. y= = and x = 6”, Equation (13) is an equation in 
which the left hand side is a function of ka/B'/? only and the right hand side is 
a function of g'”? only. Therefore it is a simple matter to find ka as a function 
of-8. For the case of a uniform field, the characteristic diffusion length in a 
sphere is a/x and inspection of Equation (1) indicates that for a uniform field, 

= r/a. k now may be considered as a measure of the effective radius of the 
discharge for diffusion and ka/w then is A/A,, where A is the characteristic 
diffusion length as determined by the geometry of the container, A, is the 
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effective characteristic diffusion length. Equation (13) is solved numerically 


and plotted in Fig. 2. 











Fic. 2. Ratio of effective diffusion length to characteristic diffusion length. 


II. Experiment 


Breakdown fields have been measured for hydrogen in a spherical cavity 
operating in the lowest electric mode. The details of experimental method are 
similar to those previously reported (5). The microwave apparatus is shown 
in Fig. 3. Microwave power with a free space wave length of approximately 
10 cm. generated by a c-w magnetron is coupled to a microwave resonant cavity 
through coaxial transmission lines. A known fraction of the power delivered 
is measured by a bolometer. The power absorbed by the cavity is combined 
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Fic. 3. Block diagram of experimental microwave apparatus. 


with the cavity Q and the known field configuration to determine the electric 
field by standard methods (1, 9). The cavities in which breakdown takes place 
are made of oxygen-free high conductivity copper and connected through 
Kovar to an all-glass vacuum system. The vacuum system would hold at a 
pressure of better than 10~* mm. of mercury for a period of about two hours 
with the pumps turned off. A single series of breakdown measurements takes 


about this time. 
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III. Calculation of Effective Diffusion Length 


Fig. 4 is a ¢ vs. E/p plot for hydrogen. Constant pd lines are plotted, the 
data being combined from breakdown measurements of cavities with different 
A’s. Those having the two smallest values of A are taken from previously 
reported data (6). . This plot illustrates the use of the theory developed in this 
paper for computing an effective A. Both the corrected and uncorrected data 


are presented. 
yee | | ". 5 


-—-—-p = CONSTANT 





AT 


‘ SPHERICAL CAVIT 
CORRECTED 


SPHERICAL CAVITY 


NUNGORRECTED 


Fic. 4. Plot illustrating the actual calculation of effective A's. Both the corrected and uncor- 
rected £ curves for the spherical cavities are shown. 


nial se aac ad sul : 


VOLTS 


E/D (ca=mm.H ) 


The procedure used in finding A, begins with a plot of ¢ vs. E/p, using the 
geometrical value which is a first approximation. If the geometrical character- 
istic diffusion length is equal to the actual diffusion length, we shall find A./A 
from Fig. 2 to be equal to 1. This is not generally true, so we take the value 
of 8 for the slope of a constant pd curve for the pressure and wave length in 
which we are interested, find A,/A and recalculate ¢ to obtain a second approx- 
imation. This process is convergent and is continued until successive calcu- 
lations agree. 


The sphere used in the experiment had a radius of 4.69 cm. and a geometrical 
A of 1.49cm. The effective A determined from the theory varies from 0.43 cm. 
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for a low E/p to 0.61 cm. for a high E/p. (The effective radius of the discharge 
depends on the ionization coefficient.) 
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Fic. 5. Comparison of experimental electric field and values predicted from theory. 


In Fig. 5 are plotted the theoretical values of E as a function of p, using the 
calculated values of A,. These are calculated from the theory developed in two 
previous papers (5, 6), the value of A, being inserted for A in this theory. 
Experimental verifications of these theories for the uniform field case was 
presented in the papers referred to. Experimental breakdown fields are com- 
pared with theory on the same figure. The good agreement illustrates the 
validity of the method. 
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TABLES OF THE FUNCTION e~**” M(a; 7; 2)! 


By A. D. MacDona.Lpb? 


Abstract 


Six-figure tables of the function u(z) = e-22/yM(a :¥ : 2) have been computed 
for the following values of the parameters, a = 0.001, 0.01, 0.1 (0.1) 1.0, 0.25, 
0.75;* y = 0.5 (0.5) 2.0, and z = 0 (0.1) 1.0 (0.5) 8. They may be used in con- 
junction with tables of the exponential function for accurate values of the con- 
fluent hypergeometric function. 


Tables of the confluent hypergeometric function multiplied by an expon- 
ential have been computed. The product e~*’’M(a; y; 2) which is herein 
designated yu(z) has appeared in analytical solutions of the Boltzmann transport 
equation for electron distribution functions and is useful in predicting electrical 
breakdown in gases (2). There are two principal reasons for computing this 
particular function. In some cases, the breakdown equation may be written 
(2) simply by setting u(z)= 2.f In addition, the confluent hypergeometric 
function has an exponential type behavior in many regions, so that interpolation 
in the tables is difficult. The function tabulated in the present work is slowly 
varying, and good accuracy is obtainable by simple interpolation, requiring 
only first and second differences. Because extensive tables of the exponential 
function are readily available, tables of u(z) may be used in place of tables of 
the confluent hypergeometric function. 


The tables of u(z) are more extensive than those of the confluent hyper- 
geometric function previously published along with some of the more important 
properties of the function (1). 


The tables have been computed by the Joint Computing Group of the 
Massachusetts Institute of Technology. 
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TABLE I 
e—*8/0-5 M4(a; 0.5; 2) 
a 
a 0.001 0.01 | 0.05 0.10 0.20 0.25 0.30 
0.1 | 1.00001 1.00007 | 1.00031 1.00054 | 1.00080 | 1.00083 | 1.00080 
0.2 | 1.00"3 1.00028 | 1.00125 | 1.00220 | 1.00324 1.00334 | 1.00317 
0.3 | 1.00006 1.00064 | 1.00288 | 1.00504 | 1.00733 | 1.00751 | 1.00710 
0.4 | 1.00012 1.00116 | 1.00524 | 1.00913 | 1.01311 | 1.01337 | 1.01257 
0.5 | 1.00019 1.00187 | 1.00839 | 1.01451 | 1.02063 | 1.02093 | 1.01956 
0.6 | 1.00028 | 1.00276 1.01237 | 1.02128 1.02993 1.03019 | 1.02807 
0.7 | 1.00040 | 1.00387 1.01724 1.02951 | 1.04106 1.04119 | 1.03809 
0.8 1.00053 1.00521 1.02308 1.03929 | 1.05406 1.05395 | 1.04962 
0.9 | 1.00070 | 1.00678 | 1.02995 | 1.05070 | 1.06900 | 1.06848 | 1.06266 
1.0 | 1.00089 | 1.00863 | 1.03793 | 1.06385 | 1.08594 | 1.08483 | 1.07720 
1.5 | 1.00235 | 1.02273 | 1.09762 | 1.15969 | 1.20322 | 1.19516 | 1.17285 
2.0 | 1.00498 | 1.04792 | 1.20093 1.31913 | 1.38328 1.35788 | 1.30833 
2.5 | 1.00941 | 1.09004 | 1.36819 | 1.56702 | 1.64142 | 1.58179 | 1.48726 
3.0 | 1.01663 | 1.15814 1.63001 1.93955 | 1.99867 1.87907 | 1.71513 
3.5 | 1.02819 | 1.26640 | 2.03268 | 2.48924 | 2.48357 | 2.26596 | 1.99940 
4.0 | 1.04654 | 1.43698 | 2.64635 | 3.29250 | 3.13459 | 2.76367 | 2.34976 
4.5 | 1.07556 | 1.70474 | 3.57747 | 4.46043 | 4.00323 | 3.39966 | 2.77842 
5.0 | 1.12142 | 2.12453 | 4.98793 | 6.15478 | 5.15835 | 4.20915 | 3.30056 
5.5 | 1.19891 | 2.78295 | 7.12434 | 8.61128 | 6.69183 | 5.23723 | 3.93487 
6.0 | 1.30864 | 3.81690 | 10.3629 | 12.1738 8.72622 | 6.54143 | 4.70421 
6.5 1.49057 5.44328 | 15.2785 17.3444 | 11.4249 8.19507 | 5.63650 
7.0 | 1.77959 | 8.00619 | 22.7505 24.8567 | 15.0057 10.2916 6.76574 
7.5 | 2.28963 | 12.0524 | 34.1258 | 35.7842 | 19.7592 | 12.9499 | 8.13326 
8.0 | 2.97325 | 18.4519 | 51.4698 | 51.6979 | 26.0728 | 16.3214 | 9.78933 
TABLE II 
€*8/1-0 4 (a; 1.0; z 
ws | | | 

N | 0.001 0.01 0.05 | 0.10 0.20 | 0.25 | 0.30 
0.1 | 1.00000 1.00003 | 1.00012 | 1.00023 1.00041 | 1.00047 | 1.00053 
0.2 | 1.00001 | 1.00010 | 1.00049 | 1.00093 | 1.00164); 1.00192 1.00214 
0.3 | 1.00002 | 1.00024 | 1.00114 1.00214 | 1.00375 | 1.00437 | 1.00486 
0.4 | 1.00004 | 1.00043 | 1.00206 | 1.00388 | 1.00677 | 1.00786 | 1.00873 
0.5 | 1.00007 | 1.00069 | 1.00330 | 1.00617 | 1.01073} 1.01244 1.01378 
0.6 | 1.00010 1.00102 1.00485 | 1.00907 | 1.01569 1.01814 | 1.02005 
0.7 | 1.00014 | 1.00142 1.00675 1.01259 | 1.02169 1.02503 1.02759 
0.8 | 1.00019 | 1.00191 1.00902 | 1.01678 | 1.02879 | 1.03314 | 1.03645 
0.9 | 1.00025 | 1.00247 | 1.01168 | 1.02169 | 1.03703 | 1.04253 | 1.04668 
1.0 | 1.00032 | 1.00313 1.01476 | 1.02735 | 1.04648 | 1.05326 | 1.05832 
1.5 | 1.00082 1.00805 | 1.03758 1.06881 | 1.11426; 1.12941 | 1.14009 
2.0 | 1.00168 | 1.01652 | 1.07638 | 1.13819 1.22405 | 1.25077 | 1.26826 
2.5 | 1.00307 1.03014 | 1.13799 1.24660 | 1.39005 | 1.43121 | 1.45566 
3.0 | 1.00525 | 1.05132 1.23253 | 1.41022 | 1.63236 | 1.69016 | 1.72002 
3.5 | 1.00857 | 1.08366 | 1.37505 | 1.65277 | 1.97949 | 2.05476 | 2.08578 
4.0 | 1.01363 | 1.13265 | 1.58800 | 2.00906 | 2.47191 | 2.56287 | 2.58644 
4.5 | 1.02128 | 1.20663 | 1.90513 | 2.53040 | 3.16713 | 3.26733 | 3.26796 
5.0 | 1.03288 1.318385 | 2.37725 3.29268 | 4.14707 | 4.24196 | 4.19330 
5.5 | 1.05049 1.48737 3.08119 | 4.40844 | 5.52852 | 5.59000 5.44880 
6.0 | 1.07729 | 1.74387 | 4.13351 | 6.04506 | 7.47838 | 7.45600 | 7.15292 
6.5 | 1.11825 | 2.138446 | 5.71152 | 8.45209 | 10.2356 | 10.04266 | 9.46838 
7.0 | 1.18108 | 2.73145 | 8.08571 | 12.0025 | 14.1426 | 13.6346 | 12.6189 
7.5 1.27781 3.64726 | 11.6698 | 17.2550 | 19.6914 | 18.6319 | 16.9123 

| 1.42726 | 5.05713 | 17.0979 27.5895 | 25.5979 | 22.7727 
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TABLE II 
‘1.0 V4 (a; 1.0; 2) 


| | 
0.40 | 0.50 0.60 0.70 | 0.75 | 0.80 


| 
| 


.00060 | 1.00063 | 1 
.00242 | 1.00250 | 1 
.00548 | 1.00563 | 1 
-00980 | 1.01003 | 1 
.01540 | 1.01569 | 1 
02231 | 1.02263 | 1 
.03057 | 1.03086 | 1 
.04020 | 1.04040 | 1 
.05124 | 1.05127 | 1 
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46374 | 1.43047 | 1 
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.04912 | 1.92526 | 1 
.49095 | 2.27959 | 2 
2 
2 
3 
3 
4 
5 
6 
7 


.00039 
00156 
.00346 | 
00608 
00939 
01337 | 
01799 
02324 
02910 
03555 
07623 
12998 | 
19599 
27387 
.36363 
.46555 
58014 
70814 
85048 
00827 
18280 
37554 
58815 
82248 


00060 | 1.00052 
00238 | 1.00206 
00534 | 1.00461 
00945 | 1.00813 
01473 | 1.01260 
02115 | 1.01802 
02872 | 1.02436 
03743 | 1.03160 
04728 | 1.03975 
05828 | 1.04878 


00046 
-13071 | 1.10695 


00183 
00409 
00719 | 
O1113 
01588 | 
.02141 
0277 

03479 
04259 

09235 
15927 
24293 
34347 
46154 
59820 
-75490 


.00192 

.00336 

.00517 

.00732 | 
.00981 | 
01262 | 
.01573 
.01913 
.04012 
.06690 
.09862 
. 13469 
.17472 
.21843 | 
. 26568 
.31637 
. 37049 
42804 
.48907 
.55365 
.62190 
. 69392 


ONO vwkwne 


23330 | 1.18656 
36855 | 1.28786 
54021 | 1.41187 
75335 | 1.56023 | 
01453 | 1.73527 
33199 | 1.93994 
71592 | 2.177 93348 | 
17882 | 2.45343 | 2.13615 
73588 | 2.77175 | 2.36550 | 
40558 | 3.13888 | 2.62455 
21031 | 3.56184 | 2.91676 
17712 | 4.04881 | 3.24607 
_33875 | 4.60927 | 3.61698 | 


.07042 | 2.72708 

.82782 | 3.28984 

.81629 | 3.99685 

.10581 | 4.88079 

.78856 | 5.98934 | 

.98605 | 7.37822 
12.8585 | 9.11894 
116.6172 [11.3019 
| 
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TABLE III 
e °/1.5 yaa; 1.5; 2) 

aru | 

fe ; 0.001 0.01 0.05 0.10 0.20 | 0.25 0.30 

0.1 | 1.00000 1.00001 1.00007 | 1.00013 | 1.00023 | 1.00028 | 1.00032 
0.2 | 1.00001 | 1.00006 1.00027 | 1.00051 | 1.00095 1.00114 | 1.00131 
0.3 | 1.00001 | 1.00013 1.00061 | 1.00118 | 1.00217 | 1.00260 1.00298 
0.4 | 1.00002 | 1.00023 | 1.00111 | 1.00213 | 1.00392 | 1.00469 | 1.00537 
0.5 | 1.00004 | 1.00037 1.00177 | 1.00339 | 1.00622 1.00743 | 1.00850 
0.6 | 1.00005 | 1.00054 | 1.00259 1.00497 | 1.00909 1.01085 | 1.01240 
0:7 1.00008 | 1.00075 | 1.00360 1.00689 | 1.01258 1.01498 | 1.01710 
0.8 | 1.00010 1.00099 1.00480 1.00917 | 1.01669 1.01986 | 1.02265 
0.9 | 1.00013 1.00129 | 1.00620 | 1.01183 | 1.02148 1.02552 | 1.02906 
1.0 | 1.00016 | 1.00162 | 1.00781 | 1.01489 | 1.02696 | 1.03200 | 1.03639 
1.5 | 1.00041 1.00409 | 1.01958 | 1.03709 1.06630 | 1.07817 | 1.08834 
2.0 | 1.00083 | 1.00821 1.03912 | 1.07364 | 1.12991 | 1.15217) 1.17084 
2.5 | 1.00148 | 1.01463 | 1.06937 1.12971 | 1.22578 | 1.26269 | 1.29292 
3.0 1.00246 | 1.02428 | 1.11454 | 4.21269 1.36513 | 1.42186] 1.46711 
3.5 | 1.00391 1.03853 | 1.18070 | 1.33319 1.56382 1.64669 | 1.71083 
4.0 1.00603 | 1.05935 1.27674 | 1.50651 | 1.84429 | 1.96103 2.04826 
4.5 | 1.00912 1.08969 | 1.41563 1.75482 | 2.23847 | 2.39845 2°51310 
5.0 | 1.01364 1.13389 | 1.61648 | 2.11047 | 2.79188 | 3.00631 | 3.15239 
5.5 | 1.02024 1.19846 | 1.90756 | 2.62076 3.56970 3.85166 | 4.03193 
6.0 | 1.02995 | 1.29313 | 2.33090 | 3.35520 | 4.66549 | 5.02961 | 5.24399 
6.5 | 1.04427 | 1.43257 2.94918 | 4.41633 | 6.21403 6.67552 6.91825 
7.0 1.06552 | 1.63899 3.85635 5.95593 | 8.40997 | 8.98248 9.23742 
1.5 | 1.09718 | 1.94606 5.19364 8.19940 | 11.5352 12.2266 12.4594 

8.0 1.14463 | 2.40514 | 7.17412 | 11.4825 | 15.9992 | 16.8035 16.9491 

TABLE IV 
e °2/2-0 \7(q: 2.0; 2) 
a 

e 0.001 0.01 0.05 0.10 0.20 0.25 0.30 

0.1 | 1.00000 | 1.00001 1.00004 | 1.00008 1.00015 | 1.00018 | 1.00022 
0.2 | 1.00000 | 1.00003 | 1.00017 1.00033 | 1.00062 1.00075 1.00087 
0.3 | 1.00001 | 1.00008 | 1.00038 | 1.00075 1.00141 | 1.00171 1.00198 
0.4 1.00001 | 1.00014 1.00069 | 1.00135 1.00254 1.00307 1.00357 
0.5 1.00002 | 1.00023 1.00110 | 1.00214 | 1.00402 1.00487 1.00565 
0.6 1.00003 | 1.00033 1.00162 | 1.00313 1.00588 1.00710 | 1.00824 
0.7 | 1.00005 1.00046 1.00224 1.00433 1.00812 1.00981 1.01137 
0.8 1.00006 1.00061 1.00297 1.00576 | 1.01076 1.01299 | 1.01505 
0.9 1.00008 | 1.00079 | 1.00383 | 1.00741 | 1.01384 1.01669 | 1.01931 
LG 1.00010 | 1.00099 1.00482 | 1.00931 | 1.01735 1.02091 | 1.02418 
1.5 1.00025 | 1.00246 | 1.01192 1.02294 | 1.04241 1.05091 | 1.05861 
2.0 1.00049 | 1.00485 | 1.02347 | 1.04500 | 1.08251 1.09863 1.11309 
2.5 1.00086 | 1.00850 1.04095 | 1.07820 | 1.14220 1.16929 1.19329 
3.0 1.00139 1.01382 1.06643 1.12632 1.22780 1.27003 | 1.30699 
3.5 1.00217 1.02146 1.10280 | 1.19467 | 1.34806 | 1.41077 1.46492 
4.0 | 1.00326 1.03230 | 1.15419 | 1.29072 | 1.51520 | 1.60525 | 1.68187 
4.5 1.00481 | 1.04762 1.22647 | 1.42510 | 1.74637 1.87266 1.97836 
5.0 1.00701 | 1.06923 1.32810 | 1.61296 2.06576 2.23983 2.38292 
5.5 1.01011 | 1.09983 1.47125 | 1.87605 2.50758 2.74452 2.93534 
6.0 1.01453 1.14329 1.67366 | 2.24575 3.12050 | 3.43996 3.69137 
6.5 1.02084 | 1.20535 | 1.96120 2.76755 3.97396 4.40159 4.72933 
7.0 1.02991 | 1.29444 | 2.37183 | 3.50767 5.16750 5.73660 | 6.15961 
‘0 1.04303 1.42309 2.96149 4.56291 6.84423 | 7.59785 | 8.13828 
8.0 1.06210 | 1.60997 | 3.813801 | 6.07539 9.21068 10.2040 | 10.8867 





| 
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TABLE Ill 


e 1-5 (a; 1.5; 2) 


0.60 70 | 0.75 | 0.80 





.00044 
.00176 
00393 
00695 
.01080 
.01548 
.02097 
.02727 
.03438 
.04228 
.09368 
. 16503 
. 25700 
.37095 
.50881 
.67316 
86727 
.09510 
.36143 
.67194 
.03438 
.45359 
.94188 
50911 


00050 | 1.00050 
00200 | 1.00199 | 
00451 | 1.00447 | 
00802 | 1.00794 | 
01254 | 1.01241 
01809 | 1.01787 | 
02467 | 1.02434 | 
03229 | 1.03182 
04096 | 1.04031 | 
05070 | 1.04983 
11607 | 1.11334 
21140 | 1.20506 
34074 | 1.32826 


00048 | 1.00050 
00194 | 1.00200 
00438 | 1.00450 
00782 | 1.00802 
01227 | 1.01257 
01777 | 1.01815 
02432 | 1.02478 
03196 | 1.03247 


.04070 | 1.04125 | 


00048 
00191 | 
00428 
00758 
O1I81 | 
01697 | 
02306 

03006 | 
03799 
04684 

10516 

18773 
29640 
43391 
60397 
81133 
06197 
36323 
.72409 
15543 
67041 

28486 
01789 
89246 


00159 
00362 

00650 
01026 

01493 
02054 
02712 

03472 
04336 
10390 
. 19829 
33537 
52730 

79064 
14804 

63035 
27970 
. 15369 
33113 
5.92009 
.06900 

-9820 
9404 


01151 
.01671 
02293 
03020 
03856 
04803 
11361 
21395 
35694 
55326 
‘81734 
16843 
63233 
24341 
04758 
10621 
50126 
34246 
7769 


0016 


.05057 | 1.05112 


Es 
ry 
1. 
1. 
E, 
a: 
1; 
r. 
L. 
Es 
.11806 | 1.11780 | 1. 
.21939 | 1.21600 | 1. 
.36100 | 1.35056 | 1. 
.55162 | 1.52806 | 1.50960 | 1.48744 | 
.80284 | 1.75715 | 1.72527 | 1.68859 
.12990 | 2.04909 | 1.99713 | 1.93938 
2. 
2. 
3. 
3. 
4. 
5. 
6. 
8. 


55278 | 2.41831 | 2.33712 | 2.24950 
(09755 | 2.88327 | 2.76037 | 2.63112 | 
.79821 | 3.46747 28592 | 3.09937 | 
69907 | 4.20080 | 3.93769 | 3.67302 
"85784 | 5.12122 | 4.74564 | 4.37534 
"34980 | 6.27693 | 5.74739 | 5.23510 
27316 | 7.72917 | 6.99005 | 6.28789 


.7562 | 9.55576 | 8.53265 | 7.57769 | 
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TABLE IV 
e *4/2-9 (a: 2.0; 
0.50 | 0.60 | 0.70 | 0.75 0.80 | 0.90 1.00 


ONO mote 


00041 
. 00166 
00374 | 
00666 | 
01043 

01507 
.02058 
02697 | 
.03426 

.04247 
.09790 | 
.17946 
.29099 
.43769 


.00042 
.00167 
.00375 
.00668 
.01045 
.01507 
.02054 
.02688 
.03409 
.04219 
.09642 
. 17520 
. 28154 
- 41952 
59452 
.81343 
.08497 
. 42008 
. 83249 
. 33929 
.96179 
. 72646 
.66634 


| 1.00040 

1.00161 
1.00364 
1.00650 

1.01020 
1.01476 

1.02019 
1.02650 
1.03373 
1.04188 
1.09736 
1. 18003 
1.29455 
1.44719 


00035 | 1.00038 
00142 | 1.00153 
“00322 | 1.00347 | 
.00576 | 1.00620 | 
00907 | 1.00975 
01317 | 1.01413 
01808 | 1.01936 | 
02382 | 1.02546 
03041 | 1.03246 | 
03789 | 1.04037 | 
08958 | 1.09464 
"16852 | 1.17652 
-28066 | 1.29140 
43402 | 1.44651 
63932 
91088 | 1.91897 
26775 | 2.26564 
.73531 | 2.71334 
34742 | 3.29075 
14921 | 4.03549 
20102 | 4.99696 
.64158 | 6.74651 | 6.58362 | 6.24014 | 6.02259 | 5.78402 | 5.26680 | 
69761 | 8.73093 | 8.40528 | 7.85055 | 7.51901 | 7.16480 | 6.41983 
.5070 [11.3973 |10.8114 | 9.94091 | 9.44437 | 8.92566 | 7.86413 


00032 
00127 | 
.00289 
00518 
00817 | 
01188 
01633 
02155 
.02756 
03440 
08201 
15560 
26142 
.40793 
.60656 
.87276 | 
.22734 
69847 
.32504 
15642 
.13163 | 5.26526 


00027 
00109 
00248 
00445 
00703 
01024 
01410 
01864 
02387 | 
02984 
(07175 
13728 
23265 
36628 
54970 
-79863 
(13457 
58702 
19652 
01897 


m-cocooco$c] | 
s 


.91345 | 1.90238 
.25139 | 2.22966 
.68462 | 2.64612 
.23913 | 3.17508 
.94873 | 3.84661 
.85743 | 4.69955 


.86581 
. 16729 
.54516 | 
.01764 
. 60788 
34521 


LOC OWN 


CUOMONOMOMNOMSONS: 


=> 
OO RWS Re Re ee eee eee eee | 








00 DOR WDD DDR RR Re RR RRR | 
> OU 69 9 DD DNDN ee 





1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
3. 1 
1.65148 | 1. | 1.64621 | 1.62637 

1 1 
2 2 
2 2 
3 3 
3 3 
a + 
6 5 
* 6 
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THE MICROWAVE ABSORPTION SPECTRUM OF 
METHYL ALCOHOL! 


By BarBARA M. GIRDWOOD? 


Abstract 


The microwave absorption spectrum of methyl alcohol was studied in the 
region 23,400 to 25,200 Mc., using a double modulation system in which the usual 
radio-frequency sine wave was replaced by a sharp, pip-shaped pulse with a 
repetition rate of 55 ke. Twenty-five new lines were found and these, along with 
those previously reported outside the region examined, were grouped into four 
series. When the methyl alcohol was mixed with a large percentage of water 
vapor the spectrum was found to be modified considerably. 


Introduction 


Lines in the microwave absorption spectrum of methyl alcohol have been 
reported previously by Hershberger and Turkevitch (5), B. P. Dailey (2), and 
by B. Bak, E. S. Knudsen, and E. Madsen (1), using either the Hughes and 
Wilson (6) Stark modulation or the Gordy and Kessler (4) double modulation 
techniques. The author used a double modulation method but replaced the 
r.f. sine wave with a sharp pip with a repetition rate of 55 kc. With this ar- 
rangement methyl alcohol and methyl alcohol in water vapor were studied in 
the 1.25 cm. region and a number of new lines were discovered. It was found 
that the presence of water vapor in large percentages had a marked effect on 
the position of the lines. 

Apparatus 


The spectroscope which was used is shown schematically in Fig. 1 and the 
pip modulator in Fig. 2. The receiver was a RCN type SCR-5A and was tuned 
to 110 kc. Frequencies were measured by a cavity wavemeter which had been 
calibrated previously from the known ammonia lines. 

As shown in Fig. 3, the irregularity of the signal from the detecting crystal 
for a sine wave and a very narrow absorption line would make it impossible to 
amplify, but the output from the pip modulator and a narrow line, as shown 
in Fig. 4, is regularly spaced and easily amplified. The signal from a line whose 
width is greater than the amplitude of the modulating wave is shown in Fig. 5. 
It is very likely that the pip modulation adds a false width to the line and thus 
distorts the line shape. 


Part of the gain in sensitivity which was obtained was due to a reduction in 
noise background. The sawtooth oscillator and pip generator were built to 
operate without a trigger signal and there was very little 55 kc. signal which 
could be picked up by the tuned receiver because the Fourier frequencies 
involved in the pip were much higher. Undoubtedly a further gain could be 
achieved by building a receiver which could be tuned to the original modu- 
lation repetition rate. 


1 Manuscript received November 8, 1949. 
Contribution from the Department of Physics, McGill University, Montreal, Que. 
2 Present address; National Research Laboratories, Ottawa, Canada. 
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Schematic diagram of microwave spectroscope. 
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Fic. 3. Klystron frequency vs. time for narrow line and sine wave modulator. 
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Fic. 4. Klystron frequency vs. time for wide line and pip modulator. 
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Fic. 5. Signal from wide line using sine wave modulator. 


The Absorption Spectrum of Methyl Alcohol 


The lines which have been reported are in a range 20,898 to 26,562 Mc. 
while the present work was done in the range 23,400 to 25,200 Mc. The fre- 
quencies of the new lines which were found, along with the ones which were 
verified, are shown in Table I. All the lines appeared as doublets. An accuracy 
of + 5 Mc. was obtained. 

TABLE I 
ABSORPTION SPECTRUM OF METHYL ALCOHOL 


w—weak, b—broad, bb—very broad, s—sharp, ss—very sharp 


Py | l 


Frequency, 


Frequency, 














Me. Reported at Comments Mc. Reported at Comments 
25,182 b 24,280 s 
25,135 25,132.0 | b 24,235 ss 
25,056 25,050 b | 24,180 ss 
25,015 25,017 bb 24,132 | ss 
24,966 24,954 | bb 24,068-76 24,080 | quadruplet 

|2 Mc separation 
24,928 24,928.8 | bb, width— || 24,017 ss 
| 24,928-913 || 
24,863 bb 23,955 ss 
24,806 b 23,908 ss 
24,752 b 23,858 23,861 ss 
24,697 b, w 23,812 s, W 
24,640 b, w || 23,759 | s, W 
24,588 b, w 23,702 Ss, W 
24,530 | | b, w 23,650 s,W 
24,478 b, w || 23,594 | s 
24,426 bb 23,538 s 
24,372 bb i 23,493 | s 


24317 | 24317 | bb _ 23,445 | 23450 | __ Triplet 
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Fic. 6. Methyl alcohol spectrum in the region 23,400 to 25,200 Mc. 
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Methyl alcohol spectrum in the region 21,000 to 27,000 Mc. 
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The very nearly equal spacings between the lines were at once apparent and, 
when the lines were empirically numbered 1, 2, 3,—in order of decreasing 
frequency and a plot made of frequency against line number (Fig. 6), the lines 
were seen to fall into three separate series. On the assumption that the lines 
which had been reported but which were outside the region:examined were 
part of a similar arrangement, the values were studied and were finally plotted 
as shown in Fig. 7. In all there were four separate series. 


The accuracy of this plot may be seen from Tables II, III, 1V, and V, where 
the first and second differences have been tabulated. The lines were grouped 
in series as indicated by the graph, and the integer dividers for the intervals 
were also determined from the graph. 





TABLE II 
SERIES | 
sis © agi First difference Second difference 
26,562 
26,310 on 
| 252 _ 50 
vo 
25,898 | 412 _ 51 si 
8 
25,796 a = 51 0 
25,546 250 _. 50 +1 
o 
See 161 ; 
25,38 | —— == -§3 —3 
; 3 
25,135 250 .. 50 +3 
0 
24,929 206 _ 5) inf 
4 
649 _ , 
24,280 — = 50 +1 
13 
24,235 45 +5 
‘24,180 | 55 | ~10 
24,132 48 +7 
24,076 | 56 —8 
—6 


Bak, Knudsen, and Madsen (1) have measured lines in the region of 25,200 
Mc. and these values, along with the others available, were plotted in Fig. 8 
to show the detail in this part of the spectrum. The above groupings include 
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TABLE III 


Frequency, 
Mc. 


25,015 


First difference Second difference 





25,966 
24,863 
24,806 
24,752 
24,697 
24,640 
24,588 
24,530 
24,478 
24,426 
24,372 
24,317 


TABLE IV 
Series III 


Frequency, 


Mc First difference Second difference 


24,017 


23,955 62 
23,908 47 
23,858 
23,812 
23,759 
23,702 
23,650 
23,594 
23,538 
23,493 
23,445 








186 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 














TABLE V 
SERIES IV 
Frec ci | ‘ - a 2 
Me a First difference Second difference 
23,415 
23,121 | 294 _ 49 
| 4 
23,033 Ss ith | ms 
22,095 938 _ 44 0 
21 
21,551 o44 45 —1 
12 
20,989 562 _ 43 2 
13 * 
20,898 91. 45 wi 
2 
Tae ke. > es —3 
Note: The identification of these as a series is uncertain owing to the size of 


the intervals. 


all the published lines except one at 24,040 Mc. which was not found by the 
author and could not be fitted into the picture. 


Alcohol in Water Vapor 


The frequencies obtained when the methyl! alcohol was in a mixture con- 
taining roughly 70% water vapor are listed in Table VI and are plotted in 
Fig. 9. The shift of the lines from their former positions—the shift being 
greater the further the line was from the region of 24,312 Mc.—and the pre- 
sence of two single lines in the same region were the most noteworthy features 
of the spectrum. 


When the mixture of alcohol and water vapor was first admitted to the silver 
absorption cell the lines were displaced but they returned to their original posi- 
tion within an hour. The second time the mixture was admitted the lines were 
displaced and remained displaced for the three hours during which they were 
watched. 
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Fic. 8. Methyl alcohol spectrum in the region 24,600 to 25,900 Mc. 


TABLE VI 
ABSORPTION SPECTRUM OF METHYL ALCOHOL IN WATER VAPOR 


w—weak, b—broad, bb—very broad, s—sharp, ss—very sharp 


Frequency | : 7 : 
{te wid Comments Frequency Comments 


24,884 ,w 24,076 ww 
24,776 24,046 s 

24,721 23,991 

24,666 23,938 

24,617 23,890 

24,570 23,835 

24,511 23,781 

24,463 23,723 

24,410 23,662 

24,364 23,606 

24,312 b, single 23,546 

24,265 b, single 23,504 

24,209 ss 23,457 

24,156 ss 23,440 

24,100 ss 
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Fic. 9. Absorption spectrum of methyl alcohol in water vapor. 


The Interpretation of the Spectra 


The methyl alcohol spectrum has been attributed by Dennison* to the 
hindered internal rotation of the OH radical and he has calculated the fre- 
quencies of eight of the lines on this basis. However this gives a value of the 
barrier height lower than that estimated from the infrared spectrum and does 
not explain the observed Stark effect.There is now the added difficulty of 
explaining the effect of the water vapor. The addition of a H atom or a com- 
plete water molecule would be expected to have a more serious effect on the 
spectrum. Another explanation might be that the spectrum is due to the 
rotation of two molecules joined together by the OH radical. Such a double 
molecule would have a rotation line at approximately 26,000 Mc. and, owing 
to vibrations of the low energy bond joining the two molecules, would have a 
number of equally spaced lines on either side of the pure rotation line. The 
number of observed lines seems to be rather large unless the energies of the 


*Unpublished results referred to by Gordy (3). 
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vibrations involved are extremely small. Further information is required 
before the methyl alcohol spectrum can be completely understood. 
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ENERGY DISTRIBUTION OF FISSION FRAGMENTS 
FROM U* AND U2! 


By D. C. BRUNTON AND G. C. HANNA 


Abstract 


New measurements have been made of-the energy distribution of fission frag- 
ments. A coincident pulse technique has been used which measures simultan- 
eously the energies of both fragments of a fission pair. Previous work in this 
field was limited in statistical accuracy owing to low fission rates. With the 
Chalk River high flux pile and a 30-channel pulse analyzer a rapid accumulation 
of data is made possible. The improved statistical accuracy has led to new 
results on the detailed association of fragment pairs and on the energy variation 
occurring in each fission mode. <A careful comparison of U®%* and U2* has been 
made. 


I. Introduction 


The first measurements of the ionization pulses from fission fragments were 
made by Frisch (13), who recorded fragment energies of the order of 100 Mev. 
in an ionization chamber in 1939. Shortly afterwards Jentschke and Prankl 
(16) resolved the fragment energy spectrum into two groups at about 60 and 
100 Mev. Booth, Dunning, and Slack (2) found a continuous distribution of 
fragment energies with two pronounced maxima. Subsequently many measure- 
ments have been made of the energy spectrum from fission fragments (8, 10, 
11, 12, 17, 18, 19, 22). 


The experiments of Jentschke (18), of Flammersfeld, Jensen, and Gentner 
(10) and of Deutsch and Ramsey (8) are of particular interest because, in these 
experiments, the energies of both fragments of a fission pair were measured 
simultaneously. ‘The source was placed at the center of a double “back-to- 
back” ionization chamber. When fission occurred the two fragments traveled 
in opposite directions into the two chambers. The ionization pulses were 
amplified and displayed on a double beam oscilloscope. A continuous photo- 
graphic record was made, and the film was examined for coincident pulses. A 
rather tedious examination of the film is required, and the results of these 
experiments show poor statistical accuracy. One of the purposes of the present 
experiment was to decrease the statistical errors and obtain reasonably accurate 
data on the spread of kinetic energies associated with each splitting ratio of the 
uranium nucleus. This objective was feasible at the Chalk River Laboratory 
because of the exceptional facilities available, namely, a high flux pile and a 
30-channel pulse analyzer. 


The work of Deutsch and Ramsey (8) included a comparison of the energy 
distributions for fission in U** and Pu**. The present experiment has extended 
this comparison to fission in U**, 

1 Manuscript received November 14, 1949. 
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As in previous work, the present experiment employed a double ionization 
chamber with a thin source mounted on the common cathode. The general 
plan of the experiment was to select in one chamber fragments in a small energy 
interval, and to measure the energy distribution of the coincident fragments 
in the other chamber. A preliminary note on some of the results obtained has 
already been published in the Physital Review (4). 


II. Experimental Equipment 


The double ionization chamber used is shown in Fig. 1. The chamber was 
constructed almost entirely of aluminum to reduce the background activity 
when irradiated with a high flux of neutrons. The insulation consisted of 
porcelain and of glass coated with ceresin wax. The electrode connections were 
introduced through Kovar seals, and the main vacuum seals were made with 
neoprene gaskets. 


The chamber employed electron collection to obtain fast rising pulses. For 
this purpose Frisch grids were used to shield the collecting electrodes from 
charges induced by the slowly moving positive ions. The operation of grid 
chambers has been carefully investigated by Bunemann, Cranshaw, and Harvey 
(5). The grids used were constructed of copper wire, 0.013 cm. diameter, 
spaced 0.188 cm. apart, and separated from the collecting electrodes by 1 cm. 
This design gave a screening efficiency greater than 98%. 


ae 
NA 


CEA 


Fic. 1. The double ionization chamber. 
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The common cathode carried the collimator. The uranium source was con- 
tained in a collodion film mounted on this collimator. Preparation of the source 
is described in Section III. The collimator was a dural plate 0.010 in. thick, 
with approximately 2000 closely spaced holes, 0.010 in. in diameter. Thus the 
collimator allows only fragments within 45° from the normal to the cathode to 
enter one chamber. Collimation is required to limit ionization losses, Losses 
associated with fragments at large angles to the normal are twofold; firstly, 
self-absorption in the source becomes objectionable with long oblique paths; 
secondly, electrons may escape to the case if a portion of the track is allowed 
to enter the fringing field near the edges of the electrodes. The cathode shields 
shown in Fig. 1 were added to reduce further this electron loss. 


The chamber was filled with pure argon plus about 3% carbon dioxide to 
a pressure of 65 cm. of mercury. The pressure was so chosen that the fragments 
of longest range traveled only about 60% of the distance from the cathode to 
the grid. Under these conditions the natural a-particles from the sources were 
also contained well within the cathode-grid region. The carbon dioxide was 
added to the argon to increase the electron mobility and permit fast collection. 


As the source was mounted on one side of the collimator only, fragments 
entering one of the chambers, ‘‘the collimated chamber’’, were collimated 
directly. Collimation in the other chamber was secured by counting only 
coincident pulses. Fragments in a small energy interval were selected in the 
uncollimated chamber by an electronic ‘gate’. Pulses from the collimated 
chamber in coincidence with these were passed to a 30-channel pulse analyzer 
(26), which measured their energy distribution. A detailed description of the 
electronic apparatus is given in Section IV(0). ’ 


The neutron source was the NRX heavy water pile at Chalk River. Most 
of the measurements were made with the chamber inside the thermal column. 
On a few occasions, where a lower flux of neutrons was sufficient, an external 
beam was used. 


III. Preparation of the Source 


The source was prepared by the evaporation of a uranium—collodion solution 
on water. This source solution was prepared by mixing four parts of a solution 
of uranium acetate (2 mgm. per ml.) in amyl acetate with one part of a solution 
of collodion in ether (87 mgm. of collodion per ml.). A small quantity of the 
final solution (~ 0.03 ml.) was dropped onto a clean surface of distilled water. 
The solution spreads out quickly and uniformly to cover an area of about 100 
sq. cm. and is allowed to dry for about 10 sec. The drying time is important 
since the uranium acetate is slowly soluble in water. With a quickly soluble 
salt (e.g., uranium nitrate) the uranium in the film may dissolve in the water 
before the film can be removed from the water surface. In preliminary trials 
with uranium nitrate and slowly drying films, up to 99% of the uranium was 
lost. Many variations of the materials used offer possibilities, but the com- 
bination of uranium acetate, amyl acetate, collodion, and water was found to 
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be the best of a wide variety tested and even with this prescription approxi- 
mately 20% of the uranium was lost to the water. 


The film was lifted from the water by means of a wire frame. A rectangular 
loop of wire was brought up underneath the center of the source. As the frame 
is withdrawn the film folds over the loop, forming a double layer. Removing 
a double film in this manner was found easier than handling a single film of 
twice the thickness. 


The film was allowed to dry for one minute and then laid carefully over the 
collimator plate. This drying time was so chosen that the film was slightly 
moist, to provide proper adhesion to the plate. Excess surface moisture gives 
rise to fillets of water in the collimator holes, which dissolves uranium out of 
the film. As the film dries, a significant fraction of the fissile material is left 
deposited on the side walls of the collimator holes. Since such a deposit can 
“see” only into one chamber, it cannot give rise to truly coincident pulses and 
contributes only an objectionable background. 


In preliminary experiments (when. the sources were prepared by different 
techniques) four-fifths of the source material collected in such deposits in the 
collimator. This gave rise to a false-to-true counting ratio of 4: 1 in the col- 
limated chamber and a significant per cent of chance coincidences. In the final 
experiments the false-to-true counting ratio was about unity and the chance 
coincident counting rate was reduced toa negligible proportion (Section IV(c)). 


IV. Experimental Conditions 


(a) OPERATION OF THE CHAMBER 


The chamber was tested under a wide range of voltage conditions. Cathode- 
grid fields from 5 to 400 v. per cm. were used. If this ‘‘collecting’’ field is too 
low, a portion of the ionization is lost through recombination, and the resulting 
pulses are not “‘saturated’’. Pulse amplitudes were measured for a series of 
field values until saturation was indicated by no further increase in pulse size 
with increase in the electric field. Saturation occurred for a-pulses in pure 
argon at about 10 v. per cm. whereas fission fragments in an argon plus carbon 
dioxide mixture required cathode-grid fields of ~ 200 v. per cm. 


lf the grid-collector field is not maintained sufficiently large relative to the 
cathode-grid field, the pulse size is decreased by loss of electrons to the grid. 
lor the design of grid used, the theoretical minimum ratio required for these 
fields is about 3: 1 according to Bunemann et al. (5). Our experimental value 
for this ratio of fields was 3.5: 1, which is in reasonable agreement with the 
theory, considering nonuniformities in the fields and inaccuracies in grid con- 
struction. To ensure complete freedom from such losses a ratio of 4: 1 was 
used in the experiment. 


The type of gas used in the chamber is important. Negative ion formation 
by electron attachment must be avoided. The relatively immobile negative 
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ions give a slowly rising pulse which is not recorded by the fast amplifiers used 
here to obtain a good signal-to-noise ratio and a low chance coincidence rate. 
Even small impurities of electronegative gases must be eliminated. Argon, 
99.9% pure, with less than five parts per million oxygen, was used. Argon, 
alone, however, suffers from the disadvantage of low electron mobility, and 
hence a long pulse rise time results. 

The pulse collection time for a track of ionization is, the time between 
the arrival of the first electrons in the grid-collector region and the arrival of 
the last electrons on the collector. Thus, if /; is the projection of the length of 
the track in the direction of the field and /, is the grid-collector spacing (/.= 


1 cm.), the collection time is 
lL, le 


. ee 


where W,; and W, are the drift velocities of electrons in the cathode-grid and 
grid-collector fields respectively. 
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Fic. 2. Block diagram of electronic arrangement for coincident pulse measurement. 


With a filling of pure argon and at field strengths of about 100 v. per cm., 
Wi= 4 X 10° cm. per sec. and W.= 8 X 10° cm. per sec. (Healey and Reed 
(14)), and consequently the maximum collection time was about 10 usec. 


The addition of a few per cent carbon dioxide to the argon increases the 
electron mobilities so that W;= 3 X 10° cm. per sec., W2= 6 X 10° cm. per 
sec., and the maximum collection time is reduced to about 1.5 ysec. This allows 
much faster amplifier response and a corresponding reduction in the coin- 
cidence resolving time. » The consequent improvement in noise background and 
chance coincident rate is discussed in Section IV(c). In addition, it was ob- 
served that with the argon plus carbon dioxide mixture less stringent conditions 
are imposed on gas purity. This may arise from the much lower agitational 
velocities of electrons in this mixture and the consequent reduction in the 
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probability of electron attachment to impurities such as water vapor and 
hydrogen sulphide. As mentioned above, however, the argon plus carbon 
dioxide mixture requires much higher fields to produce saturation. 







(6) THE ELECTRONIC EQUIPMENT 

A block diagram of the general layout of the electronic equipment is shown 
in Fig. 2. Output pulses from each chamber are fed through cathode-followers 
to a pair of linear amplifiers. The amplified output from the uncollimated 
chamber is then fed through a system of discriminators, trigger circuits and 
gating stages (D,;, Dy, T1, T2, T3, Ts, Gi, G2) by means of which pulses lying 
within a selected energy interval trigger the coincidence mixer G;. The pulses 
from the collimated chamber are converted into square pulses of the same 
amplitude by the “peak voltmeter”’ circuit (7s controls the duration of these 
shaped pulses). Of these G; passes only those that are in coincidence with the 
pulses of selected energy from the uncollimated chamber. These coincident 
pulses are then passed through a biased amplifier and their energy distribution 
is measured on a 30-channel pulse analyzer. The components of this system 
will now be described in detail. 


(1) The Cathode Followers (Fig. 3) 

As the chamber itself was in a high neutron flux, the bulk of the electronic 
equipment had to be operated some distance away. The output of the chamber 
was fed through cathode followers to avoid the loss in pulse amplitude that 
would have resulted if direct feeding of long cables had been employed. A point 
of interest is the low value of input resistance used. It is usual in this type of 
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Fic. 3. Input cathode follower and main amplifier: circuit diagram. All valves are 12SH7’s, 
filaments in series fed from a 250 v. stabilised supply. Gain control: wire wound resistors to give 
4db steps. Control R is a fine gain control for balancing the two amplifiers. 


Band width control capacities 
(condenser banks ganged) 







1 = 0.001 pf. 7 = 0.001 pf. 

2 = 0.0005 pf. 8 = 0.0005 pf. 
3 = 0.00025 pf. 9 = 0.00025 pf. 
4 = 100 wypf. 10 = 50 ppf. 

5 = 50 ups. 11 = 20 pf. 
6= 20 upf. 12 = Opp. 
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work to use resistances of about 10° ohms in the interests of good signal-to- 
noise ratio. However, first circuit noise is in this case quite negligible compared 
with the noise resulting from 8- and y-ray ionization in the chamber. On the 
other hand, the steady current from this cause flowing through a resistance as 
high as 10° ohms was found to alter appreciably the potential of the collecting 
electrode (see Section IV(d)). This slight but somewhat inconvenient depen- 
dence of chamber operating conditions on ionization was eliminated by re- 
ducing this resistance to 107 ohms. 


(2) The Main Amplifier (Fig. 3) 

The design of the main amplifier is virtually identical with that of Cranshaw 
and Harvey (6). Two high stability negative feed-back rings are separated by 
a condenser-resistance network which determines the pulse shape. 


It is important to avoid a dependence of output amplitude on pulse rise time. 
Although the 45° collimation used restricts the spread of rise times consider- 
ably, the amplifier time constants have been chosen long enough to tolerate a 
variation from zero to the maximum (1.5 usec.) in the rise time. The ‘‘smooth- 
ing’ and “‘differentiating’’ time constants used were | and 5 usec. respectively 
with a resultant pulse rise time of 2.5 usec. and decay time (from maximum to 
20% of maximum)~ 10 usec. The expected independence was confirmed 
experimentally by measuring the ratio of a-particle pulses to signal generator 
pulses (of effectively zero rise time) as a function of amplifier time constants. 


Negative pulses of amplitude ~ 50 v. were obtained at the amplifier outputs. 
The balancing of the amplifiers was achieved by the continuous gain control 
R in Fig. 3. 


(3) Coincidence Discriminator (Fig. 4) 

So far the treatment of pulses from both chambers has been identical. The 
two outputs from the amplifier now enter the coincidence discriminator where 
those pulses from the uncollimated chamber lying in a chosen energy interval 
are selected. 

Consider the energy ‘“‘gating’’ of pulses from the uncollimated chamber 
(“‘uncollimated pulses’). These pulses are first amplified by a factor of four 
and passed to two similar discriminator and trigger circuits DpT, and DyT», 
the detailed operation of which is described more fully in Appendix I. 


If a pulse amplitude exceeds the level set by D,, 7; will fire; if, further, it 
exceeds the level of Dy, T2 will fire. T: is made to cancel T; so that a gating 
pulse is passed through G; only if the actuating pulse lies between the thresholds 
of D; and Dy. This cancelling involves the use of an intermediate trigger 
circuit 7. Referring to Fig. 5, 7; is adjusted to give an output pulse of length 
just greater than the rise time, T72 (and 7) are of length about twice 73. Thus 
the back of the square wave from T; always lies within the pulse from 7, at 
whatever point on the front of the input pulse 72 is fired. The wave 7; is 
differentiated with a | usec. time constant and applied to the grid of a normally 
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cut off 6AS6 valve. The negative output from 7» is applied to the suppressor. 
Consequently there appears on the anode of this gating stage (Gi) a negative 
pulse if, and only if, D, 7, fires and DyT>2 does not. 


It was originally intended to derive from G; a square wave which would 
open the pulse analyzer to a pulse from the collimated side of the chamber (via 
amplifier B) since there was within the pulse analyzer a convenient way of 
gating pulses without introducing uncertainties into their amplitudes. How- 


High level discriminator voltage, Vy . initial Pulse from 
Amplifier A 


10 psec Pulse 
from T,- 


10 psec Pulse 
from To: 


3 psec Pulse 
from Tz 
Tz Pulse 
Differentiated. 


Output Pulse from Energy Gate, G, 


(no pulse through gate for initial pulse shown above) 


Coincident Pulse 


from 
8 Amplifier B. 


8 psec Pulse from T4 to trigger 
Collimator Gate, Go- 


10 psec Pulse from Ts 
to trigger Coincidence Gate,G3.- 


Pulse from Tg triggered by pulse from 
Amplifier B (opens peak voltmeter) 


-_---- Peak Voltmeter. 
Charging Pulse 


B (input to coincidence 
gate,G,)- 


Output Pulse from 
Coincidence Gate 


Ve ( fed through biased amplifier 
| to pulse analyser) 


10 20 30 40 psec. 


Fic. 5. Coincidence discriminator: wave-form diagram. 
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ever this would have involved a long (80 usec.) coincidence resolving time and 
resulted in a high chance coincidence rate. Consequently a coincidence stage 
of shorter resolving time was inserted between G; and the pulse analyzer. 
This—called the ‘‘Collimator Gate’’ (G2)—passes a pulse to 7; only if there 
is a coincident pulse from the collimated side, and thus assumes the onus of 
coincidence discrimination. For various reasons it was subsequently decided 
to operate the pulse analyzer at a large distance (~ 1500 ft.) from the rest of 
the equipment, and to avoid a separate lead for the gating pulses it was neces- 
sary to gate the collimated pulses on the spot using the gate G; (see below). 
This scheme was not subject to the limitations (in resolving time) of the pulse 
analyzer and the reasons for the gate G2 eliminated. However G2 was retained 
and the gate G; operated with a long resolving time. 

The gating action of Gz is accomplished as follows. Any collimated pulse 
produces a positive square wave in 74, which is applied to the suppressor grid 
of a 6AS6 valve, the control grid of which is connected to an inverted version 
of the output of G;. Both grid and suppressor are normally biased beyond 
cutoff. Consequently a pulse appears on the anode only if both grid and sup- 
pressor are pulsed simultaneously. Clearly the resolving time which deter- 
mines the chance coincidence rate is the sum of the durations of the wave 
forms from G; and 74. The former can be made very short, but there is a lower 
limit imposed on 7. The wave G; may be delayed from pulse onset time (at 
which point 7, fires) up to a time equal to the sum of the rise time and the 
duration of T;, that is, approximately twice the input pulse rise time. Clearly 
the duration of 7, must be greater than this. Actually with the rise time of 
2.5 usec. obtained from the amplifiers and 7; 3 yusec., G; was 2 usec. and 7, 
8 usec., giving a total resolving time of 10 usec. The duration of 7; (equal to 
T2) was chosen to be 10 usec. 


The trigger circuit 7; produces a 10 usec. pulse to operate the coincidence 


gate G3. 


(4) Biased Amplifier Unit (Fig. 6) 

The biased amplifier unit which contains the peak voltmeter and trigger 
circuit 7 is similar to the equipment described by Cranshaw and Harvey (6) 
with the small addition G; to enable the gating to be carried out. (For a more 
detailed description of the circuit the paper of Cranshaw and Harvey should 
be consulted. ) 

Briefly, signals entering this unit are converted into square waves of adjus- 
table length and of the same amplitude as the incoming pulses. This is accom- 
plished by the switching trigger circuit V;, V2 (Fig. 6) and the peak voltmeter 
circuit D; and C. The square-wave output is available at the cathode of V5. 
In the original design the pulses then passed directly to the biasing or “‘cut- 
ting’’ stage V, and the subsequent amplifier V;, to which has been added a 
reversal stage Vy, and an output cathode follower Vio. The advantage of using 
square waves is that the output pulse duration is independent of the amount 
of cutting, and thus subsequent amplifiers need not be of very wide band width. 
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The gating stage G; was inserted between V; and V, and operates as follows. 
Unless the circuit 7; is operated the gating input line is at +30 v. approxi- 
mately. In this condition V,, is conducting; its cathode potential is deter- 
mined by R2, so long as it is more negative than the cathode of V,4 (with this 
limitation R» acts as a ‘‘zero-set’’). If T; is triggered, by coincident pulses, the 
gating line will rise to +300 v. approximately. V,, now cuts off and V,,4 opens, 
its anode rising to the same potential as the cathode of V5, i.e., to the height 
of the pulse from the collimated chamber. 


In this experiment the biased amplifier was used merely to bias the pulses 
so that the range of fission energies extended over the whole of the pulse 
analyzer range instead of just the upper half. This of course is a simple matter 
not requiring the elaboration of the circuit used, which was designed to provide 
high resolution. However the equipment was available and was quite suitable 
for this role. 


The output cathode follower of the biased amplifier provided a convenient 
local monitoring point. The cathode circuit was completed by a 1500 ft. 50 ohm. 
transmission line, properly terminated at the far end, feeding the step-up 
amplifier and pulse analyzer. This amplifier comprised a ring-of-three circuit, 
similar to one of the stages in the main amplifier, followed by an output stage 
capable of providing the +100 v. swing necessary to drive the pulse analyzer. 
This instrument has been described by Westcott and Hanna (26). 


(c) CHANCE COINCIDENCES AND THE PULSE SPECTRUM 


Any system with finite pulse lengths and resolving times will produce results 
which contain some spurious information. In this section we shall examine the 
quality and quantity of such erroneous data in the present experiment. 


Spurious pulses arise from two causes: firstly, overlapping of pulses when 
two occur within a time defined as the pulse length and, secondly, purely 
chance coincidences. These effects become the more serious as the counting 
rate increases, and it is important to consider what limit on counting rate is 
dictated by these spurious effects. 


The phenomenon of pulse overlapping or “build-up” occurs when the mea- 
sured amplitude of a pulse is increased by the tail of a previous one. It is not 
specific to this problem but arises whenever pulse heights are being measured. 
Obviously, it results in a smearing of an energy distribution in the direction 
of higher energy; but, for a continuous distribution, the calculation of the effect 
in detail is very tedious. We limit the discussion here to a semiquantitative 
appraisal of a much simplified situation so as to demonstrate that the smearing 
effect, whatever its detailed nature, is of very small importance under the 
conditions of this experiment. 

Let N counts per second be the counting rate in the uncollimated chamber, 


CN counts per second the rate of genuine coincidences in both chambers, and 
CNR counts per second the counting rate in the collimated chamber. 
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The factor R arises from the material deposited in the collimator holes which 
cannot give rise to true coincidences as pointed out in Section III. 


A fraction, f, of the N counts per second, lying in a chosen energy interval 
is selected (f depends on the width and position of the energy gate). The genuine 
coincidence rate, CfN will be compared with the build-up and chance coin- 


cidence rate. 


1. Build-ups 


It should be pointed out that we are concerned here with the effect of pulses 
neighboring (in time) on the two pulses of a genuine coincidence only. Chance 
coincidence pairs which are affected by build-ups are automatically included 
(as a very small minority) in the chance coincidence rate (part 2). 


(i) In the Collimated Chamber 

Since the total counting rate in this chamber is CNR, the chance of a pulse 
occurring f¢; sec. before or fz sec. after any given pulse is CNR (t:+ fo). The 
genuine coincidence rate is CfN; therefore the build-up rate is CfN CNR 
(t1+ te). The first term involving ¢; expresses the chance of the true pulse 
sitting on the tail of a preceding one. The value of ¢; depends on the pulse 
shape and the amount of distortion to be considered. It will be denoted by 
ar, where 7 is the time constant of the pulse fall and a is a factor depending on 
the distortion allowed. The second term involving tf, arises from a pulse fol- 
lowing the true pulse whose amplitude is being recorded. The peak voltmeter 
circuit measures the maximum amplitude reached during the time fp+ ts; + és, 
where fp, ts, and ¢; are the rise time, and the pulse lengths in trigger circuits 
T;, and 7; respectively (see Fig. 5). If, then, a pulse occurs within this 
time after the true one and is of larger amplitude (usually by virtue of the fact 
that it is sitting on the tail of the true one) a false count will be recorded. The 
time ¢2 must include a factor determined by how much larger the distorted 
amplitude of the second pulse is, and will be denoted by p(tr+ ts+ és). 

Let us illustrate the magnitude of these effects by some numerical examples. 
For this purpose we require the decay time 7 defined as the time required for 
the pulse to decay to 20% of its maximum amplitude. The time + was mea- 
sured to be 12.5 usec. (see also Appendix I1). 

If a = 1, then build-up pulses of the first type appear as distortions of 20% 
or more in amplitude. The number of such pulses is given by CfN CNR 7 or 
as a per cent of true coincidence, CNR + X 100. Typical values are 

N = 200 counts per sec. 
C 1/40 
R = 2. 
.. per cent pulses distorted by > 20% in amplitude 
= 1/40 X 200 X 2 X 12 XK 10-*X 100 = 0.012%. 


For a = 5 it is shown in Appendix II that the pulse has decreased to 0.002% 


ll 
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of its maximum value in the time 5 7. Hence 0.06% of the pulses are build-ups 
suffering a distortion of > 0.002%. It has been assumed in these calculations 
that all pulses are of equal size. The variation in pulse size throughout the 
fission spectrum results in a larger per cent distortion on the low energy pulses 
but cannot substantially alter these results. 








For build-ups of the second kind we may use similar arguments and if we 
retain the fiction of pulses of equal amplitude the calculation is quite simple. 
In the time (fg+ tst+ ts) = (2.6 + 8 + 10) = 20 usec., we find that a pulse 
decays to 3% of its maximum amplitude. Hence if we choose p = 1 (its max- 
imum value), CNR(20 X 10~*) X 100 = 0.02% of the pulses are distorted 
by > 3%. This 0.02% includes the cases of pulses of widely different initial 
amplitude from a true pulse being recorded. Furthermore, with a continuous 
distribution in pulse size, the operation of the peak voltmeter circuit permits a 
distorted measurement only when the second pulse is of larger amplitude than 












the true pulse. 






(ii) In the Uncollimated Chamber 

The counting rate of pulses in this chamber which are coincident with col- 
limated pulses is CN. With a total counting rate N, the counting rate for 
build-ups is CN Nb 7, where 0 is a factor, similar to a above, giving a measure 
of the per cent distortion considered when the true pulse falls on the tail of a 
preceding one. The factor b, however, is more restrictive than a since the 
amplitude distortions of significance in this case are only those which introduce 
or remove pulses from the selected energy band. There is a contribution from 
build-ups of the second kind, that is when a pulse sitting on the tail of the true 
pulse opens the energy gate. This pulse must clearly lie within the time 7, (the 
coincidence resolving time, ~ 10 usec.) after the true one. Apart from this 
additional restriction the effect is much the same as that of build-ups of the 
first kind. Consequently only the latter will be illustrated numerically. 

















The energy gate selects a fraction /’ of the coincident pulses. f’ is not exactly 
equal to f but since the build-up rate decreases rapidly for large amplitude 
shifts and since in any case the number introduced into the selected band is 
approximately equal to the number removed, it is reasonable to replace f’ by f 
in practice. This means that the energy gate is opened cfN Nb 7 times per 
second by pulses (in true coincidence with collimated pulses) which have been 
distorted by build-ups. This quantity gives, therefore, the number of col- 
limated pulses which are false because they are actually in coincidence with 
pulses lying outside the selected energy range. This effect will obviously depend 
on the position and width of the energy gate. For simplicity we will again 
assume equal amplitude pulses of, say, 70 Mev. For errors greater than 5% 
these 70 Mev. pulses must trigger a gate whose lower level is 73.5 Mev. The 
effect will be computed for a 5 Mev. gate width. This gate will be opened if 
the triggering pulse occurs while the amplitude of a preceding pulse lies between 
8.5 and 3.5 Mev. This time, which is b7, above, is given in Appeiidix II as 
(18.1 — 13.7) = 4.4 psec. 
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The per cent rate of these build-ups is therefore 

Nbr X 100 = 200 X 4.4 X 10°*XK 100 = 0.09%. 
The effect of build-up in the uncollimated chamber produces a higher counting 
rate of false information than do build-ups in the collimated chamber but the 
range of error is more restricted. Thus in the case cited 1 in 1000 of the pulses 
opening the energy gate will be more than 5%, but none more than 12% in 
error. For a continuous spectrum of pulses the range of error is not thus 
restricted, but the number of pulses with a given error decreases rapidly as this 
error increases. 


2. Chance Coincidences 

For an ideal system with R = 1 there are no chance coincidences. Each 
count in the collimated chamber is accompanied by a true coincidence count 
in the uncollimated chamber and since the coincidence resolving time, 7,, is 
less than the pulse length, 7, any second pulse occurring within r, falls into the 
classification of a build-up. Thus it is only the source deposits within the col- 
limator holes that give rise to chance coincidences and the counting rate for 


these spurious pulses in the collimated chamber is C(R — 1)N. 


The counting rate from the uncollimated chamber at Gz isfN. Consequently, 
the chance coincidence rate is fN CN(R — 1)r-. 

As a per cent of the true coincidence rate the chance rate, N(R — 1)7., 

= 200 x 10 X 10°*X 100 = 0.2%, 
for r-= 10 usec. and R = 2 (see Section III). 

In contrast to build-up distortions, the chance coincidence pulses are dis- 
tributed in a normal fission spectrum. As mentioned in Section III, the source 
conditions were much worse in preliminary experiments where R was 5.4. In 
addition, pure argon was used as a chamber filling for these runs and conse- 
quently a resolving time 7, of 32 usec. was required. Furthermore a counting 
rate of N = 500 counts per second was used. The resulting per cent of chance 
coincidences, 500 X 4.4 X 32 K 107° 100 = 7%, produced a considerable 
distortion in the distributions. 


(d) THERMAL NOISE AND 8-y BACKGROUND 

There are two causes of the ‘‘noise’’ apparent at the amplifier output: firstly, 
thermal noise in the cathode follower and the first circuit of the main amplifier 
and, secondly, pulses from the 8-y background in the chamber. The first 
effect is much smaller than the second. The equivalent value in kev. for the 
standard deviation of a signal generator pulse distribution obtained in the 
absence of 6-y background was about 65 kev. This is about twice the 
optimum value attainable, but is reasonable considering that the circuit em- 
ployed involved two “‘first circuit’”’ noise sources and a low input impedance. 


The second effect, under the conditions of irradiation employed in the 
experiment, gave rise to an equivalent standard deviation of signal generator 
pulses of 450 kev. The effect may be understood semiquantitatively by the 
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following argument. Suppose (for various reasons)  §-particles traverse the 
chamber per second giving an average pulse size corresponding to an energy E. 
When x is large compared to 1/7, where 7 is a measure of the pulse duration, 
a steady current flows which is subject to the usual statistical fluctuations of 
random phenomena. The standard deviation would then be of the order of 
(nr) in the number of 8 pulses per amplifier pulse duration; that is, a Gaussian 
distribution in amplitude of pulses (of a frequency ~ 1/7) of standard devi- 
ation og= E(nr)'. 


In argon a fast 8 particle produces about 90 ion pairs per cm. The average 
track length in the chamber is about 10 cm. Thus on an energy scale corres- 
ponding to 28 ev. per ion pair, E = 90 X 10 X 28 = 25 kev. +r = 107 sec. 
The quantity 7 is a little difficult to determine but one may make an estimate 
from this knowledge that there is ~ ly per thermal neutron in the pile flux 
and this may well swamp any local 8-y process (in an aluminum chamber). 
If we assume an efficiency for the production of 6’s by y quanta of 107 anda 
thermal neutron flux of pv = 10° neutrons per sq. cm. per sec. over an area of 


10 sq. cm. 


n = 108X 10?X 10-?= 108 sec. 
nr = 10° 
og ~ 750 kev. 


(The mean current = » X 90 X 10 X 1.6 X 107% = 1.5 & 1078 amp.) 


(€) CALIBRATION AND DISCUSSION OF ERRORS IN THE 
ABSOLUTE ENERGY MEASUREMENTS 


Measurements of energy were made effectively at the chamber collecting 
electrodes by feeding pulses to these points from a signal generator through 
two small identical condensers (Fig. 3). If a signal generator pulse of amplitude 
V volts is applied through a coupling condenser of C farads it will produce a 
pulse at the input of the amplifier equal to the collection of a charge due to 
ionization of VC coulombs, independent of the total input capacity provided 
C is always connected. Clearly if | and C are both known this provides an 
absolute method of calibration. Pulses from the signal generator were thus 
treated in exactly the same way as ionization pulses, and any nonlinearity in 
the amplifier system was readily avoided. 


Attenuators in the signal generator allowed the pulse amplitude to be varied 
over a wide range in accurately known small steps. Such pulses fed into the 
collimated side were recorded on the pulse analyzer, thus calibrating the pulse 
analyzer channel against signal generator millivolts, Fig. 7(a@). These signal 
generator pulses were likewise fed into the uncollimated chamber. For each 
pulse amplitude, the firing level of the low level discriminator D; was increased 
until 7, just failed to fire. In this way D, was calibrated in terms of signal 
generator millivolts, Fig. 7(b). Calibration of Dy followed automatically from 
its voltage relation to D;. It is worth mentioning here that, since the equipment 






















206 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


associated with each chamber was calibrated directly from the signal generator, 
there was no necessity for accurate balancing of the two main amplifiers. 


Absolute calibration was obtained by measuring the natural a-particle ener- 
gies from the sources in terms of signal generator millivolts and using the known 
energies of these particles (U** — 4.823 Mev., U?*% — 4.396 Mev.) to obtain 
the Mev. equivalent of one millivolt. 


Such a method assumes, of course, that the number of electron volts per ion 
pair, w, in argon is the same for fission fragments as for a-particles. As yet 
no accurate experiments have been performed to test this assumption. 
Henderson (15) measured fission energies by a calorimetric method and ob- 


FIG. 7 (a) 
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Fic. 7. (a) Calibration curve for pulse analyzer. 
(b) Calibration curve for discriminator. 


tained 177 + 1.8 Mev. This measurement included part of the 6 decay energy 
of the fragments but very little of the neutron, y or high energy £6 radiation. 
Henderson estimated that the kinetic energy was 12 + 6 Mev. less or 165 + 
8 Mev. lonization chamber measurements give a somewhat lower value. The 
mean of the values obtained by Deutsch and Ramsey (8) and by the present 
experimenters (see Table I) is 152.8 Mev. These results appear to be the most 
reliable since in each case the source was very thin and the emergent angle 
of the fragment track was restricted by a collimator. In our case there is a 
correction of about 2% (see below) to be made for source thickness and col- 
limator loss: a similar correction would seem to be necessary in Deutsch and 
Ramsey’s experiment. The corrected mean value then becomes 155.8 Mev. 
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Theoretically one would expect w to increase for fission fragments (relative 
to a-particles) since the heavy fragment has a velocity only 0.6 times that of a 
4.5 Mev. a-particle. If we assume that w is the same for fission fragments and 
a-particles of the same velocity, w for the heavy fragment will be 2% greater 
than for a 4.5 Mev. a-particle, but for the light fragment will be almost the 
same. This variation of w with a-particle velocity is discussed in the paper of 
Cranshaw and Harvey (6). The value of the total kinetic energy thus corrected 
is 157 Mev. We estimate the limits of error here as about + 3 Mev., and there 
would seem consequently to be no serious discrepancy between the ionization 












and calorimetric measurements. 







TABLE I 
NEUTRON FISSION IN [235 






SLow 






















Jentschke | Flammersfeld, | Deutsch| Fowler | Present 
and Jensen, and | Jentschke| and and | experi- 
Prankl Gentner | | Ramsey | Rosen | ment 
Most probable energy of 
light group, Mev. 91 92 92.5 94 92.5 92.7 
Most probable energy of| | 
heavy group, Mev. a7 59 65 60 61.2 59.0 
Ratio of most probable! | 
energies 1.60 1.56 1.47 } 57 1.51 1.57 
Most probable mass ratio | 
of coincident pairs (on 
an equal ratio-interval 
curve) - - - 1.49 - 1.485 
Width at half-maximum of 
high energy peak, Mev. 17 15 13 12 15 12 
Width at half-maximum of 
low energy peak, Mev. 22 19 | 20 19 24.5 20 
Ratio of peak heights 1.3 1.22 1.49 1.57 1.46 1.37 
Foil thickness, mm. air 
equivalent 0.05 0.1 
(backing 0.34)| 0.15} 0.1 0.2 0. 

















The 2% correction to the energy measurements mentioned above arises from 


two sources: (1) loss of energy by self-absorption in the source, and (2) loss of 






ions in the collimator. 






The energy absorption in the source may be evaluated from the value 6.5 
Mev. per mm. given by West (25) for the initial energy loss of an average 
fission fragment in standard nitrogen. A source (85% collodion, 15% uranium) 
of thickness 14 ugm. per sq. cm. is the equivalent of ~ 0.1 mm. of air and 
hence gives rise to an energy loss of 0.6 to 0.9 Mev. depending on the path of 
the particle through the source. The mean loss from this cause is, therefore, 
~™ 0.75 Mev. or ~ 1%. The a-particle loss from the same cause is ~ 0.2%. 
Since the fission energies are determined by comparison with a-energies, the 
error in fission energies from this cause is about 0.8%. 









The collimator losses arise from the recombination of ions produced in the 
weak field region inside the collimator holes. For the 0.010 in. collimator and 
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65 cm. gas pressure used the energy loss (assuming complete recombination 
but neglecting obliquity of track) is at most 0.010 X 25.4 X 6.5 K 0.9 = 1.25 
Mev. or 1.7%. The corresponding a-particle loss is 0.4% and hence the net 
error is 1.3%. The combined error from both causes is therefore ~ 2%. 


V. Experimental Results 


(a) RESULTS witH U*%* 
The energy distribution was measured, first, with the energy gate wide open 
so that the whole range of fragment energies was recorded. The result is the 
usual double-humped curve shown in Fig. 8. The details of the distribution 
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8 


| 


aR 


tt 


1 
i | i i 
4 SO 60 70 80 930 0 IO 
ENERGY ( Mev.) 


Fic. 8. Spectrum of fragment energies with U?*, 
Distribution measured with energy gate wide open. 


were then examined by measuring the distributions associated with an energy 
gate 5 Mev. wide. The gate was set on the low energy end of the spectrum and 
moved up in discrete steps of 5 Mev. In each case the coincident distribution 
was recorded. The results are shown in Fig. 9 (a—m). In these figures the 
normal distribution curve of Fig. 8 is shown by a broken curve on a reduced 
scale, and the gating position relative to this curve is indicated. The curves 
are normalized to equal counting times: the standard error (of counting) com- 
puted from the unnormalized curves, is shown for the maxima only. The curves 
are uncorrected for the source and collimator losses discussed in Section IV(e). 
The corrected values are: 

Most probable energy of light fragment = 94.5 Mev. 

Most probable energy of heavy fragment = 60.2 Mev. 


(b) REsuLts witH U2 


The experimental conditions using the U*** source were identical with those 
for U** except that the higher isotopic purity of the source allowed the use of 
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FIG.9 (m) uss 
ENERGY GATING BAND 
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Fic. 9(m). Spectra of U?* (energy gate 5 Mev. wide). 


9 


a lower neutron flux for a reasonable counting rate. This was a practical 
advantage since the chamber background ionization, which was proportional 
to the neutron flux, gave considerable ‘‘sparking’’ difficulties at high fields 
and high fluxes. The total distribution for a wide open gate is given in Fig. 10, 
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Fic. 10. Spectrum of fragment energies with U?, 
Distribution measured with energy gate wide open. 


and the partial distributions for a sequence of runs with a 5 Mev. wide gate 
are shown in Fig. 11 (a — k). Again the curves are uncorrected for source and 
collimator losses. The corrected values are: 

Most probable energy of light fragment = 93.0 Mev. 

Most probable energy of heavy fragment = 56.6 Mev. 


VI. Analysis of Experimental Results 


(a) DEPENDENCE OF GATED DISTRIBUTIONS ON GATING ENERGY 


In Table I the results on U?® derived from the normal distribution curves 
are compared with those of previous workers. In this table the fragment 
energy values have been left uncorrected for source and collimator loss to form 
a better comparison with previous work which was likewise uncorrected. It 
will be noted that the agreement is very satisfactory and that in particular the 
widths of the distributions indicate adequate resolution in the present ex- 
periment. 
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The results of the gated runs are rather interesting. When the energy gate 
is set on the heavy fragment group, the energy distribution of the corresponding 
light group is almost independent of the position of the gate, and conversely. 
This result is shown in Tables II and III for the two fissionable materials. 


TABLE II 
ANALYSIS OF GATED RUNS WITH U?% 


Me: 7 Coincident distribution 
Mean energy 
of gate, — eer eens a 
Mev. Most probable Width, Mev. 

energy, Mev. (at half-maximum) 


Light fragment group 


Hed 
ie 


Ungated 


47.0 
52.0 
56.0 
60.2 
64.6 
68.9 
73.3 


AE = 26.3* 


Sowa w¢ 


Heavy fragment group 


Ungated 


io. 


83.4 


AE = 4.3* | 


* AE is the maximum spread for the values in each column. 


In each case the partial distribution covers nearly the whole range of the 
complete spectrum of one group. The distributions are not identical, however, 
and the shift that does occur is such that as the gating energy is increased, the 
corresponding distribution maximum also increases. This relationship between 
the gating energy and the energy of the corresponding maximum is shown in 
Fig. 12. 


Flammersfeld ef al. (10) presented their results in a manner similar to Figs. 
9 and 11, but their statistical accuracy was too low to define accurately the 
position of the distribution peaks, and no curve similar to Fig. 12 is given. 
Neither Jentschke (18) nor Deutsch and Ramsey (8) have attempted to show 
such distribution peaks as the statistical accuracy in both cases is, again, too 
low. The situation may be summarized by the statement that a heavy frag- 
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ment.of lower than average energy for the heavy group will be paired, on the 
average, with a light fragment of /ower than average energy. The corresponding 
situation with mass distribution curve is quite different. Since the sum of the 
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TABLE III 
ANALYSIS OF GATED RUNS WITH U3 













Coincident distribution 





Mean energy 
of gate, ee eames 
Mev. Most probable | Width, Mev. 
energy, Mev. | (at half-maximum) 










Light fragment group 








Ungated 










46.4 Oo 9. 

50.9 89.6 11.2 
55.7 90.7 12.4 
60.5 90.7 10.9 
65.5 90.3 11.6 
70.1 93.1 12.6 
76.2 Ve ‘ 














AE=298 | AE=5.3 | AE=3.0 











Heavy fragment group 










Ungated 









76.2 2 4 
82.1 54.7 15.3 
86.6 55.4 17.7 
91.5 56.6 16.8 
95.8 5 1.5 






AE = 19.6 
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masses is constant, a heavy fragment of higher than average mass for the heavy 
group must be paired with a light fragment of Jower than average mass. Thus 
the double-humped energy distribution curve is not even approximately a 
simple inversion of the double-peaked mass distribution curve, as was first 
pointed out by Jentschke and Prankl (17). 


Further indication of the relationship between fragment energies may be 
found by examination of Figs. 9(g) and 11(g). Here the gating band is set on 
the minimum between the two peaks of the total distribution. If the energy 
curves were an inverse picture of the mass curves, this experiment would have 
resulted in a single maximum at the same energy as the gate or perhaps two 
small maxima very close to this energy. The experimental result, however, is 
a curve with two maxima close to the energy peaks of the total distribution. 


From Tables II and III it may be observed that the width of the gating band 
(~ 5 Mev.) can have very little effect on the measurements. Since a shift of 
mean gating energy of 5 Mev. produces only a small change in the distribution, 
the spread of energies about the mean through a 5 Mev. band will have con- 
siderably less effect. Some runs made with a 2.5 Mev. gate width verified this 
conclusion. 


(6) Mass Ratio CURVES 


Fig. 13 shows the distribution of mass (or energy) ratio for U?* and U**, 
The mass ratio is the reciprocal of the energy ratio, since by conservation of 


momentum E,;M;= E,My, where E and M refer to mass and energy and the 
subscripts ZL and // refer to light and heavy fragments respectively. The curves 
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Fic. 13. Mass ratio curves for U* and U8, 


were computed by determining directly from Figs. 9 and 11 the number of 
counts in equal energy intervals for a given energy (or mass) ratio and con- 
verting these figures to equal mass ratio intervals. The most probable ratios 
are 1.485 and 1.60 for U** and U** respectively. This ratio for U may be 
compared with 1.49 found by Deutsch and Ramsey (8). 
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(c) Mass DISTRIBUTION 






Figs. 14 and 15 show the mass number curves for U** and U** respectively. 
These curves were computed from the mass ratio curves of Fig. 13 by assuming 
total masses of 236 and 234 units in the two cases and converting the frequency 
of occurrence to equal mass intervals. The mass curves from isotopic identi- 
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Fic. 14. Mass curve for U?* from present experiment and comparison with chemical curve. 
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fication of fragment abundance (the ‘chemical method’’) are taken from 
Plutonium Project report (21) and from the paper by Steinberg, Seiler, 
Goldstein, and Dudley (23). The agreement is not very satisfactory, and 
especially in regard to the comparison of the peak positions in the U*** curves 
it is definitely poor. In this case the light fragment peaks differ by approx- 
imately 4 and the heavy fragment peaks by 9 mass units. A coincident fission 
fragment experiment yields, of course, the mass distribution before neutron 
emission whereas the chemical method applies to the fragments after neutron 
emission. Thus a shift towards higher mass must be expected in a coincident 
pulse experiment. The observed shifts, however, indicate that neutron emission 
occurs predominantly from the heavy fragment. This conclusion is contrary to 
the evidence obtained by DeBenedetti, Francis, Preston, and Bonner (7) in 
their experiment on the angular correlation of fission neutrons. The results of 
Deutsch and Ramsey (8) give no clear indication of sucha shift but the accuracy 
in their results is hardly better than + 5 mass units. 


Examination of the curves of Steinberg ef a/. indicates an accuracy of + 2 
mass units in drawing the curves from the experimental points. In addition 
the sum of the most probable masses is 93 + 136 = 229 which indicates a shift 
of 5 mass units from the coincident pulse experiment where a total mass of 234 
units is assumed. The accuracy of the peak positions from the present experi- 
ment is estimated as + 2 mass units. It has been found that the mass curves 
plotted from preliminary data taken with considerably different experimental 
conditions (different source, collimator of same geometry but six times as thick, 
chamber filled with pure argon only, resolving time 30 usec.) are in agreement 
to this accuracy. 


If the heavy fragment pulses in the chamber were slightly unsaturated while 
the light fragment pulses were saturated, the correction to be applied would 
be in the direction of removing the discrepancy. This condition, however, 
would have been apparent during the course of the experiment from the energy 
ratio of the fragments in the most probable case. This ratio was checked with 
widely varying collecting fields, and even when all pulses were well below 
saturation (70% of their maximum value) the ratio of the two peaks was 
unchanged. However, the measured energy ratio may be in error from another 
cause. As pointed out in Section IV(e), there is good reason to believe that 
the heavy fragment produces less ionization per Mev. than the light fragment. 
This would reduce the energy (or mass) ratio. A difference in w of 2% was 
estimated in Section IV(e), but it would have to be as high as about 5% to 
remove the discrepancy between the chemical data and the coincident pulse 
experiment. 

(2d) ContouR DIAGRAMS 


Figs. 16 and 17 are contour diagrams of the modes of fission. The energies 
of the fragments are given by the rectangular co-ordinates and the total kinetic 
energy and mass ratio are given by the oblique co-ordinates. The contours 
represent relative frequency of occurrence. These curves are plotted directly 
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from the data of Figs. 9 and 11 and are uncorrected for source and collimator 
errors. Neglecting differences in w for the two fission fragments the most pro- 
bable mode of fission in U** is one of total kinetic energy 156 Mev. (corrected 
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5 


Fic. 18. Variation of the most probable total kinetic energy with mass ratio. 


value) and mass ratio 1.57. Corresponding values for U*** are 150 Mev. and 
1.64. The two diagrams are on the whole very similar, although in U** the 
constant frequency curves are more nearly circular. 


(e) VARIATION OF KINETIC ENERGY WITH Mass Ratio 


Fig. 18 shows the most probable total kinetic energy as a function of mass 
ratio. The energy increases as the ratio decreases to a value ~ 1.25. From 
1.25 to symmetrical fission, however, the energy decreases. Calculations of 
total fission energy (kinetic + excitation energies of the fragments) from 
nuclear masses (e.g., Bohr and Wheeler (1)) show that the total energy increases 
continuously to unity mass ratio. Consequently, in the region between ratios 
of 1.25 and unity the excitation energy is increasing rapidly. This decrease in 
kinetic energy confirms the work of Katcoff, Miskel, and Stanley (20), who 
require such a condition to explain their range measurements of plutonium 
fission fragments. 


VII. Discussion of Experimenta! Results 
(a) RANGE OF KINETIC ENERGIES IN FISSION 


The wide range of kinetic energies released in different modes of fission is 
somewhat surprising. From the contour diagrams (Figs. 16 and 17) the vari- 
ation in kinetic energy over all modes is about 60 Mev., while the maximum 
variation associated with a fixed mass ratio is nearly as great (~ 50 Mev.). 
For a fixed mass ratio, the total energy release in fission is determined by the 
difference in mass energy of the uranium nucleus plus a neutron and the mass 
energy of the final stable nuclei. This variation in total kinetic energy for a 
fixed mass ratio must arise either from variations in the number of neutrons 
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“ 


emitted or from variations in the primary charge division. The first alternative 
has been discussed by Way and Wigner (24), the latter has been considered by 
one of the present authors (3). 


(b) COMPARISON OF SLOW NEUTRON FISSION IN DIFFERENT ELEMENTS 


Table IV shows a comparison of energy distribution measurements in 
different materials fissile by slow neutrons. The figures for Pu are taken 
from the work of Deutsch and Ramsey (8). Again the uncorrected values are 
used for comparison. There is an indication here that fission becomes more 
symmetrical as the mass of the fissile nucleus increases. At the same time the 
total kinetic energy increases, which might be expected since the total available 
energy also increases as symmetry is approached. 


TABLE IV 


COMPARISON OF FRAGMENT ENERGY MEASUREMENTS 
FOR THREE FISSILE MATERIALS 


Present Deutsch 
experiment and Ramsey 


Pu? 


Most probable energy of the light fragment, Mev. ¢ | ‘O2. 93 
Most probable energy of the heavy fragment, Mev. 65.5 | & 65 
Ratio of most probable energies 1.43 
Most probable mass ratio (from an equal ratio interval 

curve) 1 
Width at half-maximum of high energy peak, Mev. | 12 
Width at half-maximum of low energy peak, Mev. 
{ Total energy 


Most probable fission mode 
\ Mass ratio 
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APPENDIX I 






The discriminator proper consists of V; and the associated bias setting 
potentiometer, Fig. 19. Consider first the operation with C2 omitted. If the 
input pulse is of greater amplitude than the positive bias applied to the cathode 
of V;, the current flowing through V34 via R» will be diverted through V. 
Thus if the cathode of V; falls v» volts below earth, a current of (300/R.2+ 
vo/ Rez) will flow into C;. If this happens for a time 7, ” times per second there 
will be a net current through V; of m7(300 + vo)/Re. The circuit is completed 


4 aoe : 5 7 nrtR, 
via R,, and consequently a shift in the bias level will occur given by — 


9 










(300 + v0). This is small here since 7 is only 150,:and in the worst case rt = 
20 usec. and vp= 200 v., which gives 0.75 v. in 200 v. 
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Fic. 19. Discriminator and trigger circuit. This is a detail from the main coincidence dis- 
criminator circuit, Fig. 4, and shows Dx and T, (or Du and T:). 
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If substantially higher counting rates were to be employed it would be 
desirable to eliminate this effect by inserting a cathode follower between CiR; 


and V3. 


It should be noted that no d-c. restoration of pulses to the bias line is 
attempted. Such d-c. restoration is superfluous because of the low counting 
rate used here. (If counting rates are kept low enough for the probability 
of pulses sitting on each others’ tails to be small d-c. restoration is always 


unnecessary). 


R2 is made fairly small, that is, the standing current through V3, = 7 is large, 
so that steeply falling pulses can be dealt with. If the pulse falls at a rate 
dvo/dt, the self-capacity of V; will give rise to a current Cy,dvo/dt. If this is 
less than 7 only a small change in the triggering level will result. If it exceeds 7 
the grid of V2 will be driven negative (and the trigger circuit fired) while V1 
is still nonconducting. The value i = 3 ma. enables pulses falling at the rate 
of 200 v. in 0.2 usec. to be tolerated (Cy,= 3yuyf.). 


However this small change in current will cause a very small potential change 
across V;, sufficient to trip the trigger circuit since this is extremely sensitive. 
Rather than decrease this sensitivity it was decided to maintain the facility of 
dealing with fast pulses by adding a neutralizing condenser C2. This is con- 
nected to a positive ‘‘antiphase’”’ version of the pulse which is already available 
earlier in the complete circuit. In view of the behavior described above it is 
clear that fast pulses are acceptable even if C2 is not exactly correct in value, 
and with antiphase pulses of slightly different shape (in practice an almost 
inevitable difficulty). 

Vo, Vs form a very sensitive trigger circuit which responds to any drop 
greater than ~ 1/10 v. in the potential of the grid of V2. 

The main point of interest is the self-stabilization of the anode potential of 
V2, which is achieved by d-c. negative feedback (R; and Rs) to the suppressor 
grid, which in this valve has good control of the anode current. The grid of V4 
may then be run very close to the point at which conduction starts in V4 
without the usual danger of instability. It has previously been usual in circuits 
of this type to avoid this danger either by heavy biasing of V, or by running V2 
so that the anode potential in the quiescent state lies below the ‘‘knee” of the 
pentode i,/V, characteristic, with a consequent loss of gain in this stage. 
Either of these solutions means a large, and consequently somewhat uncertain 
triggering threshold. 

In this circuit the potential of the suppressor grid of V2 is —1 + 1 v. de- 
pending on the health of the valve; the anode potential is 148 + 2 v. if R; 
and R, are exactly equal, which is of course not necessary. 

By means of R; and Rs the grid potential of V4 is set 7.5 v. below the sup- 
pressor potential of V2, that is, V4 is run with a bias of 8.5 + 1 v. 


The gain between grid and anode of V2 is about 100 for fast pulses (C, 
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“holding’”’ the suppressor) and thus the uncertainty of 1 v. in the level of V4 
grid is reflected as a 1/50 v. uncertainty in the triggering level at the grid of 
V2. Since it requires a positive pulse of amplitude about 4 v. on the grid of V4 
to cause conduction the triggering level is about 1/12 v. 


The grid of V4 will respond to fast pulses by virtue of C; which is made 
approximately equal to the input capacity of V4, thus giving an essentially 
aperiodic coupling network. 

The time constant C,R, determines the length of trigger pulse produced. 
If, however, the operating pulse holds V2 in the cutoff state for a longer period 
than this, it prevents the circuit recovering. Consequently in this case a pulse 
is produced of duration equal to the time V2 is cut off. This behavior is fre- 
quently an advantage as it avoids double triggering on any spurious long 
pulse. 

If, as in this case, the full gain of V2 is to be exploited, that is, the valve run 
right up to its rating, it is particularly important to consider the screen dis- 
sipation. A wide range of valve performances was encountered in practice. 
Consequently a circuit designed to work satisfactorily with the poorest valve 
will, if a good valve is inserted, divert a large current to the screen. This pre- 
cludes running the screen from a constant voltage source. Rather R; must be 
made large enough (and a consequently higher supply line voltage used) so 
that the screen dissipation is automatically limited. This means that at high 
counting rates the screen potential will fall and the triggering level will change, 
but this effect is exceedingly small in the present application. With this excep- 
tion there will be no change in any operating conditions with repetition fre- 
quency. 

The reason for splitting the anode load of V4 is merely that a relatively low 
output was needed and the capacitative loading on V4 could thus be reduced. 


APPENDIX II 
Pulse Duration 
The pulse rise and pulse length are determined by the main amplifier time- 
constants—the equivalent circuit is shown in Fig. 20, where we neglect the 
loading of C2R2 on CiR. 


i= 9m&g 


Fic. 20. Equivalent circuit for main amplifier band width control. 


The pulse duration will depend slightly on the shape of the pulse rise from 
the electron collection in the chamber, but when the amplifier rise time con- 





226 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


stant is longer than the collection time a close approximation to the true pulse 
duration is obtained using a step function input pulse. 

A simple analysis of the transient response of this circuit is achieved using 
the Laplace Transform method (9). 

For a unit step function input voltage e, the output voltage eo is given by 


T2 oe tien os) ‘ 
a (e 1/8» ven). 
EE = FE 


where G = gmRi, 7:= CiRi and r2= C2Rz. 


To calculate the time at which the pulse reaches its maximum amplitude 
we set deo/dt = 0. That is, 


eo= e,G 


‘. ae 
— etm’? 4. —_ etm’: = (0, where t = tm, for €9= €0 maz 
T2 T1 


T2 
*, ta = —— Ind, ‘where 
A- 1 
2 r/i => 
oo €o mat — Ge rd > 


= K(e’" —e-/") where K 


€0 max 


= 2Ke~™ sinh Bt where a 


Case (a): time f; for pulse to decay to 20% of its maximum amplitude. 


‘ €o 
1.e., — = 0.2, 


eo max 
ty 
e - 
—— = 10K. 
sinh Bf, 
Values of 7 and 72 in the amplifier are somewhat uncertain owing to stray 
capacities. The values 7;= 1 ysec., t2= 5 usec. fit the measured pulse shape 
reasonably closely. With these values \ = 5, K = 1.87, a = 0.6, B = 0.4, 


then, fy = 11.25 usec. 
Case (b): t:= 57 = 56.25 usec. 


€o aie cee os = 
— = 3.1/4 e *-% sinh 22.5 


€0 max 
2X 1075. 

Case (c): t1= 20 usec. 

es 3.74 e~® sinh 8 

€o mar 

0.03. 

Case (d): time br (see Section IV(c), 1(ii)) for a 70 Mev. pulse to decrease 
from 8.5 Mev. to 3.5 Mev. 
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8. 3.5 


Let ¢; and tz be the times for the pulse to reach oe = 0.12 and 70 = 0.05 
7 7 


of the maximum respectively, 

then f; is given by —_—_—_——_—- = 
sinh 0.4 fy 

from which 


similarly ———— = 75, 
sinh 0.4 te 


from which to = 18.15 usec. 


br = floe—t = 4.4 psec. 
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